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I.  Introduction 

Plasma  Enhzutced  Chemical  Vapor  Deposition  (PECVD)  is  a  process  whereby 
thin  films  of  amorphous,  microcrystalline  and  crystalline  films  my  be 
deposited.  As  the  activation  energy  for  the  production  of  radicals  in  PECVD  is 
provided  by  electron  Impact  collisions  in  the  plasma,  the  substrate 
temperature  may  be  significantly  lower  in  PECVD  than  that  used  in  thermally 
activated  chemical  vapor  deposition  (CVD) .  Hence.  PECVD  is  attractive  for 
applications  which  are  thermally  sensitive. 

A  typical  PECVD  reactor  and  oi>erating  conditions  are  shown  in  Fig.  1. 
PECVD  is  usually  performed  at  low  pressures  {<  0.5  Torr)  in  a  parallel  plate 
capacitlvely  coupled  radio  frequency  (rf)  discharge.  The  power  deposition  is  a 
few  to  100 's  of  mW-cm'^  at  frequencies  of  100 's  of  kHz  to  10 's  of  MHz.  Feed 
stock  gases  are  flowed  through  the  discharge  with  residence  times  of  10 ’s  of 
ms  to  1  s.  resulting  in  an  energy  investment  of  1  -  lO's  of  eV  per  molecule. 
Electron  impact  dissociation  and  ionization  collisions  create  neutral  amd 
charged  radicals  which  diffuse  or  drift  out  of  the  plasma  to  the  substrate. 
While  approaching  the  substrate,  the  radicals  may  undergo  further  chemical 
reactions.  At  the  substrate,  the  radicals  adsorb,  bond,  sputter,  and  otherwise 
react  with  the  surface  to  form  the  film.  Saturated  molecules  may  be  desorbed 
from  the  surface  and  enter  the  plasma. 

The  PECVD  of  amorphous  hydrogenated  silicon  (a-Si:H)  has  received  a  great 
deal  of  attention  during  the  past  few  years  as  this  material  is  used  in  the 
fabrication  of  thin  film  transistors,  photo-voltaic  devices  (eg.,  solar  cells) 
and  photo-receptors,  as  used  in  copiers.  In  this  process,  feedstock  gases 
consisting  of  mixtures  of  silane  {SiH4),  disilane  (SiaH^).  emd  buffers  (eg., 
Hj,  Ar)  are  used  to  generate  radicals  such  as  SiHs  and  SiH2  which  are  the 
major  contributors  to  film  growth.  The  a-Si:H  film  consists  of  a  network  of 
Si-Si  and  Si-H  bonds  having  an  atomic  hydrogen  fraction  of  .5-20%.  Device 
quality  material  may  be  deposited  at  growth  rates  of  lO's  -  lOO’s  of  A/mln. 

The  processes  leading  to  the  deposition  of  a-Si:H.  described  below,  are 
complex.  As  a  result,  there  has  in  the  past  been  no  a  priori  method  to  predict 
film  properties  (ie..  growth  rate,  hydrogen  fraction,  hydr i de/d 1 hydride  bond 
ratios,  surface  roughness)  given  physical  parameters  such  as  gas  mix,  power 
deposition  and  reactor  geometry.  The  technology  associated  with  thin  film 
growth  has  therefore  been  an  empirical  discipline,  and  the  optimization  of 
film  properties  has  been  slow  to  occur. 

The  PECVD  process  is  conceptually  divisible  into  three  interlinked 
subprocesses : 

•  Gas  phase  electron  kinetics 

•  Heavy  particle  plasma  chemistry  and  transport 

•  Surface  deposition  kinetics. 

Gas  phase  electron  kinetics  includes  the  processes  in  which  power  is  coupled 
through  the  free  electrons  and  deposited  in  the  gas  by  electron  impact 


•  Pressure:  10’s  -  100’s  mTorr 

•  Power:  10’s  -  100  mW/cm® 

•  Electrode  Separation:  1-10  cm 

•  Electron  Density;  10®-  10^°/cm^ 

•  Bulk  Electron  Temperature:  1-2  eV 

•  Gas  Residence  Time:  10’s  ms  -  1  s 

•  Fractional  Dissociation;  1-10’s  % 

Figure  1  -  Typical  plasma  reactor  conditions  for 
low  pressure  PAMP  (eg.,  etching,  deposition) 
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collisions.  To  describe  these  processes,  one  must  be  able  to  specify  the 
electron  distribution  function  for  energy,  position  in  the  reactor  and  time 
during  the  rf  cycle:  f(e.x,t).  The  heavy  particle  plasma  chemistry  amd 
transport  includes  the  processes  during  which  radicals  amd  ions  are  generated 
by  electron  impact  amd  further  reactions  occur  between  the  radicals  amd  with 
the  feed  stock  gases.  These  processes  also  include  the  transport  of  radicals 
and  charged  species  in  the  plasma  by  either  diffusion,  convection,  or  drift  of 
charged  species  in  the  electric  field.  The  surface  deposition  kinetics  include 
those  processes  occuring  on  and  near  the  substrate  which  result  in  the 
formation  of  the  film. 

Although  conceptually  divisible,  these  three  subprocesses  are  strongly 
coupled.  The  electron  Impact  dissociation  of,  say,  silame  creates  radicals 
which  by  further  reactions  chemge  the  electron  transport  coefficients  amd 
hence  the  electron  distribution  function.  This  chamge  in  the  electron 
distribution  function  causes  a  chamge  In  the  rate  of  radical  generation. 
Radicals  and  ions  impacting  on  the  substrate  form  the  aunorphous  silicon  film. 
While  doing  so,  these  processes  generate  species  which  are  released  back  into 
the  plasma  which  affect  the  plasma  and  electron  kinetics. 

To  obtain  a  predictive  capability  for  the  PECVD  of  a-Si-H,  one  must 
address  the  three  sub-processes  described  above  in  a  self  consistent  fashion. 
In  this  research  project,  a  computer  model  for  the  PECVD  of  a-Sl:H  has  been 
developed  consisting  of  interlinked  sub-models  for  the  electron  kinetics, 
plasma  chemistry,  and  surface  deposition  kinetics.  A  Monte  Carlo  simulation 
is  used  to  generate  the  electron  distribution  function  and  electron  impact 
rate  coefficients  as  a  function  of  position  and  time  in  a  parallel  plate 
capacitlvely  coupled  rf  discharge.  These  rate  coefficients  are  used  as  input 
to  a  1-1/2  dimensional,  time  dependent  plasma  chemistry  model.  The  plasma 
chemistry  model  produces  the  densities  and  fluxes  of  radicals  and  ions.  These 
fluxes  are  then  used  as  input  to  a  model  for  the  surface  deposition  kinetics 
which  simulates  bulk  and  surface  properties  of  the  film.  This  work  represents 
the  first  iteration  of  a  totally  self  consistent  model  for  the  PECVD  of 
a-Si :H. 
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II.  Research  Results 


The  results  of  this  research  project  are  sunmarized  in  publications 
included  in  the  Appendix.  The  abstracts  from  those  publications  are  reproduced 
below. 


On  the  Balance  Between  Silylene  and  Silyl  Radicals  in  RF  Glow  Discharges  in 
Silane:  The  Effect  on  Deposition  Rates  of  a-Sl:H  (M.  J.  Kushner,  J.  Appl. 
Phys.  62.  2803  (1987)) 


Abstract 

Plasma  Enhanced  Chemical  Vapor  Deposition  (PECVD)  of  amorphous 
hydrogenated  silicon  {a-Sl:H)  is  typically  performed  in  low  pressure 
{i  0..5  Torr)  radio  frequency  (RF)  discharges  in  gas  mixtures 
containing  silane  (SiH.4).  TTie  initiating  step  is  electron  impact 
dissociation  of  silane  whose  products  are  primarily  the  silylene 
(SiH2)  and  silyl  (SIH3)  radicals.  The  ratio  of  SiH2  to  SIHq  in  the 
plasma  is  Important  because  these  radicals  Incorporate  differently 
into  the  a-Si:H  film;  and  therfore  the  characteristics  of  the  film 
are  a  function  of  the  ratio  [SiH3]/[SiH2] •  The  initial  step  is 
followed  by  a  series  of  hydrogen  abstraction,  silylene  insertion, 
and  silyl  association  reactions  which  subsequenty  alter  this  ratio. 
The  branching  ratios  for  electron  impact  dissociation  of  silane  as 
well  as  the  rate  constants  for  the  subsequent  reactions  have  not 
been  measured  or  are  uncertain.  Using  a  model  for  am  RF  discharge  in 
silane  gas  mixtures,  the  effects  of  branching  ratios  for  silame 
dissocation  and  rate  constants  for  key  reactions  are  investigated. 

We  find  that  in  order  for  SIH3  to  be  the  dominant  radical,  as 
currently  thought,  its  bramching  ratio  from  electron  impact 
dissociation  of  SIH4  must  exceed  0.7.5.  and  the  yield  for  H  atoms  in 
the  bramch  for  SiHa  must  also  exceed  0.75.  The  deposition  rate  of 
a-Si:H  is  then  controlled  by  the  yield  of  H  atoms  cmd  the  subsequent 
generation  of  radicals  by  hydrogen  abstraction  from  SlH^.  These 
results  require  that  the  silyl  association  reaction  must  have  a  rate 
constant  <,  10'“cm^s"‘,  amd  that  the  rate  constamt  for  insertion  of 
SiH2  into  silane  must  exeed  10*“cm^s"‘. 


A  Phenomenological  Model  for  Surface  Deposition  Kinetics  During  Plasma  and 
Sputter  Deposition  of  Amorphous  Hydrogenated  Silicon  (M.  J.  Kushner,  J.  Appl. 
Phys.  62,  4763  (1987)) 


Abstract 

The  surface  processes  during  the  plasma  enhanced  chemical  vapmr 
deposition  and  reactive  sputter  deposition  of  amorphous  hydrogenated 
silicon  (a-Si:H)  are  investigated  by  use  of  a  phenomenological 
model.  The  model  consists  of  an  accounting,  in  rate  equation  form, 
of  adsorption  of  radicals  from  the  plasma  onto  the  surface,  surface 
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diffusion,  incorporation  into  the  lattice,  interconnection  of  bonds 
in  the  lattice,  eind  burial  of  species  on  the  surface,  thereby 
constituting  film  growth.  By  accounting  for  the  coordination 
p>artners  of  Si  atoms  in  the  film,  the  atomic  fraction  of  hydrogen  in 
the  film  is  computed  for  both  the  lattice  and  for  hydrogen  in 
polymeric  or  isolated  configurations.  Results  from  the  model  are 
discussed  while  parameterizing  the  probability  for  hydrogen 
elimination  during  incorporation,  and  the  probability  for  saturation 
of  dangling  bonds  by  gas  phase  species.  We  find  that  the  mode  of 
hydrogen  elimination  during  incorporation  distinguishes  films  grown 
dominantly  from  SiH2  or  SlHa  radicals.  Qiaracteristics  of  films 
grown  by  sputter  deposition  are  Investigated  as  a  function  of  the 
composition  of  the  radical  flux.  We  find  that  films  grown  from 

hydrogen  rich  fluxes  are  composed  dominantly  of  dihydride  {=Si.  ) 

configurations,  whereas  hydrogen  lean  mixtures  are  composed  of 
dominantly  hydride  (=S1-H)  configurations. 


A  Model  for  the  Discharge  Kinetics  and  Plasma  Chemistry  During  Plasma  Enhanced 
Chemical  Vapor  Deposition  of  Amorphous  Silicon  (M.  J.  Kushner,  J.  Appl.  Phys. 
63,  2532  (1988)) 


Abstract 

A  model  for  the  plasma  enhanced  chemical  vapor  deposition 
(PECVD)  of  amorphous  hydrogenated  silicon  (a-Si:H)  in  rf  and  dc 
discharges  is  presented.  The  model  deals  primarily  with  the  plasma 
chemistry  of  discharges  sustained  in  gas  mixtures  containing  silane 
(SiH^).  The  plasma  chemistry  model  uses  as  input  the  electron  impact 
rate  coefficients  generated  in  a  separate  simulation  for  the 
electron  kinetics  and  therefore  makes  no  apriori  assumptions  as  to 
the  manner  of  power  deposition.  Radical  densities  and  contributions 
to  film  growth  are  discussed  as  a  function  of  gas  mixture,  electrode 
separation,  and  locale  of  power  deposition;  and  comparisons  are  made 
to  experiment.  A  compendium  of  reactions  and  rate  consteuits  for 
silane  neutral  euid  ion  chemistry  is  also  presented. 


Phase  and  Energy  Distribution  of  Ions  on  Electrodes  in  Radio-Frequency 
Discharges  (G.  Hebner  emd  M.  J.  Kushner.  J.  Appl.  Phys.  62,  22.56  (1987)) 


Abs  t  rac  t 

A  Monte  Carlo  particle  simulation  and  parametric  models  for  the 
sheath  voltage  and  thickness  were  used  to  calculate  the  arrival 
phase  and  energy  of  ions  striking  the  electrodes  in  low-pressure 
capacitively  coupled  rf  discharges.  Ion  phase  and  energy 
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distributions  are  presented  as  a  function  of  rf  frequency,  ion  mass, 
rf  voltage,  dc  bias,  sheath  thickness,  and  gas  pressure.  When  the 
rf  frequency  is  below  the  ion  response  frequency,  the  ions  arrive  at 
the  electrode  in  phase  with  the  applied  rf  voltage.  As  the  rf 
frequency  Increases,  the  highest  ion  arrival  probability  shifts 
towards  higher  phase  until,  at  sufficiently  high  frequencies,  it  is 
nearly  uniform.  The  transition  from  a  highly  peadced  ion  phase 
distribution  at  low  frequencies  to  a  uniform  distribution  at  high 
frequencies  requires  at  least  an  order  of  magnitude  change  in  rf 
frequency.  The  implication  of  these  calculations  on  the  electron 
energy  distribution  is  discussed. 


Important  Considerations  for  Optimizing  Production  Rates  in  RF  Discharge 
Chemistry  (M.  J.  Kushner  and  A.  Garscadden  in  L.  G.  Christophorou  and  D. 
Bouldln,  Gaseous  Dielectrics  V  (Pergammon,  New  York.  1987)  pp.  334-341) 

Abs  t  rac  t 

Important  issues  which  should  be  considered  when  optimizing  RF 
discharges  for  use  in  plasma  processing  are  discussed  using  results 
from  electron  kinetics  and  plasma  chemistry  models.  Mixtures  of 
Ar/SiH^  and  H2/SiH4  are  discussed  as  exemplary  systems  for  the 
plasma  deposition  of  amorphous  silicon.  The  effect  of  dissociation 
of  the  feed  stock  gas  on  subsequent  rates  of  dissociation  is 
discussed,  as  is  the  effect  of  diluting  silane  discharges  with 
argon. 


Simulation  of  the  Bulk  and  Surface  Properties  of  Amorphous  Hydrogenated 
Silicon  Deposited  from  Silane  Plasmas  (M.  J.  McCaughey  and  M.  J.  Kushner, 
submitted  to  J.  Appl.  Phys.) 


Abstract 

A  Monte  Carlo  simulation  for  the  growth  of  amorphous 
hydrogenated  silicon  (a-Si'H)  thin  films  by  plasma  enhanced  chemical 
vapor  deposition  (PECVD)  is  presented.  The  focus  of  the  model  is  to 
predict  the  bulk  properties  of  films  having  thicknesses  of  <  2000A 
(eg.,  hydrogen  content,  buried  hydride/di hydride  ratios,  porosity, 
and  surface  roughness),  the  effects  on  the  film  properties  of  the 
composition  of  the  radical  flux  incident  on  the  surface  are 
extunined.  Film  properties  were  found  to  be  critically  dependent  on 
the  ratio  of  SiH3/SiH2  in  the  radical  flux,  with  high  values  for 
this  ratio  causing  depositon  resembling  chemical  vapor  deposition 
and  low  values  resembing  physical  vapor  deposition.  Rough  films  (rms 
roughness  of  lO’s  of  A)  result  from  radical  fluxes  having  high  SiH2 
fractions.  We  find  that  surface  roughness  and  hydrogen  fraction 
increase  with  increasing  growth  rate  and  increasing  film  thickness. 
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though  thin  films  have  large  hydrogen  fractions  due  to  there  being  a 
hydrogen  rich  surface  layer  with  a  thickness  approximately  equal  to 
the  surface  roughness.  We  also  find  correlation  between  porosity 
(subsurface  voids)  and  hydrogen  fraction,  implying  that  hydrogen  is 
trapped  in  voids  as  a  result  of  their  being  lined  with  ^i-H 
configurations. 
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III.  Publications  and  Technical  Presentations 

The  following  publications  and  techninal  presentations  were  prepared 
under  this  contract. 


Book  Chapters 

L.  E.  Kline  sind  M.  J.  Kushner.  "Computer  Simulations  of  Materials  Processing 
Plasma  Discharges",  to  be  published  in  CRC  Critical  Reviews  in  Solid  State  and 
Materials  Science. 


Journal  Articles 

M.  J.  Kushner,  "Mechanisms  for  Power  Deposition  in  Ar/SiH^  Capacitively 
Coupled  RF  Discharges,"  IEEE  Trans.  Plasma  Scl.  PS- 14.  188  (1985). 

G.  Hebner  and  M.  J.  Kushner.  "The  Phase  and  Energy  Distribution  of  Ions 
Incident  on  Electrodes  in  Radio  Frequency  Discharges."  J.  Appl.  Phys.  62,  2256 
( 1987) 

M.  J.  Kushner,  "On  the  Balance  Between  Silyle  and  Silylene  Radicals  in  RF  Glow 
Discharges  in  Silane:  The  Effect  on  Deposition  Rates  of  a-Si:H,"  J.  Appl. 
Phys.  62,  280.3  (1987) 

M.  J.  Kushner,  "A  Phenomenological  Model  for  the  Surface  Deposition  Kinetics 
During  Plasma  and  Sputter  Deposition  of  Amorphous  Silicon,"  J.  Appl.  Phys.  62. 
4763,  (1987) 

M.  J.  Kushner,  "A  Model  for  the  Discharge  Kinetics  and  Plasma  Chemistry  During 
Plasma  Enhanced  Chemical  Vapor  Deposition  of  Amorphous  Silicon",  J.  Appl. 

Phys.  63  .  25.32  (1988) 

M.  J.  McCaughey  and  M.  J.  Kusliner,  "A  Simulation  of  the  Bulk  and  Surface 
Properties  of  Amorhous  Hydrogenated  Silicon  Deposited  From  Silane  Plasmas", 
submitted  to  J.  .Appl.  Phys. 


Invited  Conference  Presentations  with  Proceedings 


M.  J.  Kushner,  "A  Plasma  Chemistry  and  Surface  Model  for  the  Deposition  of 
a-Si:H  from  RF  Glow  Discharges:  A  Study  of  Hydrogen  Content,"  Plasma 
Proceedings.  Symposia  Proceedings,  vol.  68,  J.  W.  Coburn.  R.  A.  Gottscho  and 
D.  W.  Hess,  Editors.  Material  Research  Society.  Pittsburgh,  pp.  293-307,  1986. 

M.  J.  Kushner  and  A.  Garscadden  "Important  Considerations  for  Optimizing 
Production  Rates  in  RF  Discharge  Chemistry",  Gaseous  Dielectrics  V. 

Proceedings  of  the  Fifth  International  Symposium  on  Gaseous  Dielectrics, 
Knoxville,  Tennessee  1987,  L.  G.  Chr Istophorou  and  D.  W.  Bouldin,  eds. 
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A  Monte  Carlo  particle  simulation  and  parametric  models  for  the  sheath  voltage  and  thickness 
were  used  to  calculate  the  arrival  phase  and  energy  of  ions  striking  the  electrodes  in  low- 
pressure  capacitively  coupled  rf  discharges.  Ion  phase  and  energy  distributions  are  presented  as 
a  function  of  rf  frequency,  ion  mass,  rf  voltage,  dc  bias,  sheath  thickness,  and  gas  pressure. 

When  the  rf  frequency  is  below  the  ion  response  frequency,  the  ions  arrive  at  the  electrode  in ' 
phase  with  the  applied  rf  voltage.  As  the  rf  frequency  increases,  the  highest  ion  arrival 
probability  shifts  towards  higher  phase  until,  at  sufficiently  high  frequencies,  it  is  nearly 
uniform.  The  transition  from  a  highly  peaked  ion  phase  distribution  at  low  frequencies  to  a 
uniform  distribution  at  high  frequencies  requires  at  least  an  order  of  magnitude  change  in  rf 
frequency.  The  implication  of  these  calculations  on  the  electron  energy  distribution  is 
discussed. 


I.  INTRODUCTION 

Low-pressure  radio-frequency  ( if)  discharges  are  now 
commonly  used  for  plasma  processing  of  semiconductor  ma¬ 
terials.''^  In  order  to  predict  the  characteristics  of  the  plas¬ 
ma  reactor  and  the  materials  produced,  studies  have  been 
performed  of  the  electron  kinetics  in  such  systems.^'^  The 
results  of  theoretical^’  and  experimental^'^  investigations 
have  indicated  that  the  electron  energy  distribution  is  a  func¬ 
tion  of  space  and  time.  In  general,  the  majority  or  “bulk”  of 
the  electron  energy  distribution  is  composed  of  low-ene»'gy 
electrons  with  a  small  “tail”  of  high-energy  electrons.  The 
high-energy  tail  of  the  distribution  results  from  secondary 
electrons  produced  by  ion  bombardment  which  are  acceler¬ 
ated  by  the  sheath  potential,  and  from  the  acceleration  of 
bulk  electrons  by  the  oscillating  sheath  boundary.  Since  a 
portion  of  the  high-energy  tail  of  the  electron  energy  distri¬ 
bution  is  the  result  of  ion  bombardment,  knowing  the  time 
and  energy  distributions  of  the  ions  striking  the  electrode  is 
important  to  accurately  describe  the  temporal  and  spatial 
behavior  of  the  tail  of  the  electron  energy  distribution.  In  this 
paper  we  report  on  a  theoretical  investigation  of  the  time  and 
energy  dependence  of  ions  bombarding  the  electrodes  in  a 
low-pressure  rf  discharge.  Previous  studies  of  ion  bombard¬ 
ment  of  the  electrodes  in  rf  discharges  have  only  examined 
the  distribution  of  energies  or  angles  of  the  ions  that  strike 
the  electrode.^'  ’  As  a  result,  the  main  emphasis  of  this  study 
is  the  investigation  of  the  arrival  time,  or  phase  of  the  ions 
with  respect  to  the  applied  rf  voltage,  and  the  energy  distri¬ 
bution  of  the  ions  that  strike  the  electrode  after  traversing 
the  sheath. 

The  model  is  described  in  Sec.  II,  followed  by  a  discus¬ 
sion  of  the  results  for  frequency,  ion  mass,  voltage,  and 
sheath  thickness  in  Sec.  III.  Concluding  remarks  are  in 
Sec.  IV. 

II.  DESCRIPTION  OF  THE  MODEL 

The  study  was  performed  using  a  Monte  Carlo  simula¬ 
tion  similar  to  that  previously  described.'’  Therefore,  the 


model  will  only  be  discussed  briefly.  The  voltage  drop  across 
the  sheath  and  the  sheath  thickness  are  assumed  to  have  the 
form 

K(f)  =  max(  sin  eo^t.V^ )  ( 1 ) 

and 

/(/)  =  max(J',fc,  +  sin  )  ,  (2) 

where  Vl,t)  is  the  time  dependent  voltage  across  the  sheath, 
/(r)  is  the  sheath  thickness,  and  is  the  rf  frequency.  The 
subscripts  on  V and  X  denote  to  the  time-independent  (dc) 
and  the  time-dependent  (rf)  components  of  the  applied  vol¬ 
tage  and  sheath  thickness,  respectively,  max  {a,b)  indicates 
the  maximum  value  of  a  or  6.  As  our  concern  is  with  electro¬ 
positive  plasmas  where  the  plasma  potential  is  always  posi¬ 
tive,  the  function  max  is  used  to  guarantee  a  positive  sheath 
potential.  andAT,..  as  equal  to  0.001.  The  spatial  distri¬ 
bution  of  the  electric  held  within  the  sheath  is  chosen  to  be 
linear  and  has  the  form 

E{x,t)  =  EoVit)  - x]  ,  (3) 

where 

fo  =  2(F(f)//(f)]  .  (4) 

Using  the  above  parametric  models,  the  trajectories  of 
the  ions  in  the  sheath  are  computed  using  standard  Monte 
Carlo  techniques.'*  '*  The  ions  are  injected  into  the  simula¬ 
tion  with  a  random  velocity  chosen  from  a  Maxwellian  dis¬ 
tribution.  Their  initial  spatial  positions  are  greater  than  the 
maximum  sheath  thickness  (X^c  +  X^).  As  previously  de¬ 
scribed,'*  elastic  and  charge  exchange  collisions  are  includ¬ 
ed  in  the  simulation  by  specifying  a  mean  free  path  A,  and 
randomly  choosing  a  collision  distance  7^,  such  that 

4  =  -/I  ln(r)  ,  (5) 

where  r  is  a  random  number  on  (0,1).  When  the  ions  strike 
the  electrode,  their  energy  and  arrival  time  with  respect  to 
the  applied  rf  voltage  are  recorded. 

Simulations  were  performed  in  helium  and  argon  for 
frequencies  in  the  range  of  0. 1-40.0  MHz,  rf  voltages  of 400- 
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800  V,  dc  biases  of  0-300  V,  gas  densities  of  1.0  X  10'*-S.0 
X  10'*  cm  '  ^  and  sheath  thicknesses  of  0.5-1. 5  cm. 

III.  RESULTS  AND  DISCUSSION 

Unless  noted  otherwise,  the  discharge  conditions  for  the 
simulation  of  helium  and  argon  ions  striking  the  electrodes 
were  a  peak  rf  voltage  (.V^)  of  400  V,  maximum  sheath 
thickness  (X^)  of  0.5  cm,  and  gas  density  of  l.Ox  10'* 
cm  ~ .  The  dc  bias  ( )  and  the  time-independent  sheath 
thickness  )  were  set  to  zero  for  the  first  cases  discussed 
below.  For  the  calculations  in  helium,  an  energy-indepen¬ 
dent  charge  exchange  collision  cross  section  of  1.5  X  10~ 
cm^  and  an  elastic  collision  cross  section  of  l.Ox  10~'*  cm^ 
were  used.  We  found  that  the  results  were  insensitive  to  the 
choice  of  the  elastic  collision  cross  section  as  long  as 
The  calculations  with  argon  ions 
were  performed  with  energy-independent  charge  exchange 
and  elastic  collision  cross  sections  of  2.5  X 10  ~  "  cm^ .  For 
both  gases,  the  initial  ion  velocity  was  chosen  to  be  close  to 
the  Bohm  criterion  for  an  electron  temperature  of  10  000  K. 
Trials  showed  that  variation  of  the  initial  ion  velocity  had  a 
very  small  effect  on  the  final  distributions.  Good  statistics 
were  obtained  with  runs  of  10  000  particles. 

A.  rf  frequency  and  ion  mass 

The  dependence  of  the  ion  energy  distribution  and  the 
arrival  phase  at  the  electrode  of  light  and  heavy  ions  is  shown 
in  Figs.  1-4.  Figures  1  and  2  show  the  distribution  of  arrival 
times  for  He  ions  and  Ar  ions  as  a  function  of  frequency  and 
phase  with  respect  to  the  rf  cycle.  A  phase  of  O-g  represents 
the  cathodic  part  of  the  cycle,  while  a  phase  oiir-ltr  denotes 
the  anodic  part  of  the  cycle.  We  define  a  parameter,  the  ion 
response  frequency,  <U/  =:  1/r,  where  r  is  the  time  required 
for  the  ion  to  cross  the  sheath  while  being  accelerated  by  the 
average  sheath  potential.  The  parameter  a,  for  He  and  Ar 
ions  are  approximately  3.6  and  0.76  MHz,  respectivly.  For 
low  rf  frequencies,  that  is  frequencies  that  are  significantly 


FIG.  I.  Relative  toul  arrival  probability  of  He  ions  striking  the  electrodes 
as  a  function  of  phase  and  rf  frequency  for  a  peak  rf  voltage  of  400  V,  nuxi- 
mum  sheath  thickness  ofO.S  cm,  and  a  pressure  of  0.3  Torr.  A  phase  of  0  to 
<r  represents  the  cathodic  part  of  the  cycle. 
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PHASE  (It) 

FIG.  2.  Relative  total  arrival  probability  of  Ar  ions  striking  the  electrodes  as 
a  function  of  phase  and  rf  frequency  for  the  conditions  listed  in  Fig.  I . 


below  <U/,  the  ions  follow"  the  applied  rf  voltage  with  the 
majority  of  the  ions  arriving  in  phase  with  the  voltage.  The 
probability  of  ions  arriving  at  the  electrode  during  the  anodic 
part  of  the  cycle  is  small  since  the  sheath  region  is  swept  free 
of  ions  during  the  cathodic  part  of  the  cycle.  Ions  must  then 
diffuse  into  the  sheath  region  from  the  plasma. 

As  the  frequency  increases,  the  phase  of  the  highest  ion 
arrival  probability  shifts  towards  higher  phase  since  the  ions 
no  longer  have  sufficient  time,  during  a  single  rf  cycle,  to  be 
accelerated  and  traverse  the  sheath  in  phase  with  the  voltage. 
At  sufficiently  high  rf  frequency,  ions,  due  to  their  finite 
mass,  can  no  longer  respond  to  the  time  varying  field  but 
instead  respond  to  a  time-averaged  field.  As  a  result,  the  ions 
tend  towards  a  uniform  distribution.  However,  even  at  fre¬ 
quencies  much  higher  than  u, ,  the  ion  distribution  is  not 
uniform  but  is  skewed  towards  the  cathodic  portion  of  the 
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FIG.  3.  Relative  arrival  probability  of  He  ions  striking  the  electrodes  ax  a 
function  of  phase,  energy  range,  and  rf  frequency  for  the  conditions  listed  in 
Fig.  I.  Within  each  small  block  the  horizontal  axis  is  phase  from  0  to  Iv, 
while  the  vertical  axis  is  relative  arrival  probability.  The  higher  energy 
ranges  have  been  multiplied  by  a  display  factor  of  3  or  6. 
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FIG.  4.  Relative  arrival  probability  of  Ar  ions  striking  the  electrodes  as  a 
function  of  phase,  energy  range,  and  rf  frequency  for  the  conditions  listed  in 
Fig.  I.  Within  each  small  block  the  horizontal  axis  is  phase  from  0  to  27, 
while  the  vertical  axis  is  relative  arrival  probability.  The  higher  energy 
ranges  have  been  multiplied  by  a  display  factor  of  3  or  6. 
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cycle.  The  transition  from  a  peaked  behavior  at  low  frequen* 
cies  to  a  quasiuniform  distribution  at  high  frequencies  does 
not  happen  rapidly,  as  illustrated  in  Fig.  S.  The  transition 
requires  at  least  an  order  of  magnitude  in  rf  frequency. 

As  noted,  the  high-energy  tail  of  the  electron  energy 
distribution  is  influenced  by  the  ion  bombardment  of  the 
electrode  because  of  the  production  of  secondary  electrons 
which  are  then  accelerated  across  the  sheath  potential.  As 
the  peak  of  the  ion  phase  distribution  moves  towards  higher 
phase,  the  production  of  high-energy  secondary  electrons 
will  also  move  towards  higher  phase.  Since  the  tail  of  the 
distribution  is  responsible  for  the  optical  excitation  in  the  He 
rf  discharge,*  the  optical  emission  from  the  discharge  should 
show  some  shift  with  increasing  frequency.  This  predicted 
phase  shift  of  a  portion  of  the  high-energy  tail  of  the  electron 
energy  distribution  is  suggested  by  recent  experimental  re- 


FIG.  S.  .Measured  full  width  at  half  maximum  of  the  relative  total  arrival 
probabilities  shown  in  Figs.  I  and  2. 


FIG.  6.  Optical  emisaion  ( SS7S 
A,r  •  14  na)  1  cmbdowtherf 
driven  electrode  in  a  0.5-Torr 
He  tf  diachatfe  at  two  frequen¬ 
cies.  The  rf  voltage  on  the  driv¬ 
en  electrode  ia  aa  ahown.  The 
negative  dc  adf hiaa  waa  aO  V 
for  0.4  MHz  and  a  1/3  the 
zero  to  peak  voltage  for  2.3 
MHz.  IncTcaaing  optical  emia- 
aton  ia  plotted  downward.  The 
electro^  acperation  ia  13  cm. 
The  rf  circuit  ia  described  in 
Ref.  6.  At  0.4  MHz  the  optical 
emisaion  peaks  approximately 
in  phase  with  the  cathodic  por¬ 
tion  of  the  applied  rf  voltage. 
At  2.3  MHz,  the  peak  emission 
is  shifted  to  higher  phase.  The 
optical  emiaaion  at  2.3  MHz  is 
longer  due  to  the  negative  dc 
self-bias. 

suits  using  a  high  speed  framing  camera.*  The  time-depen¬ 
dent  optical  emission  from  a  rf  discharge  has  also  been  ob¬ 
served  using  an  apertured  photomultiplier  tube  (PMT)  and 
an  interference  filter  (Fig.  6).  These  results  show  that  at  low 
(0.4  MHz)  frequencies,  the  optical  emission  peaks  during 
the  cathodic  half  cycle  and  is  approximately  in  phase  with 
the  applied  rf  voltage.  During  the  anodic  half  cycle  the  opti¬ 
cal  emission  is  small.  At  higher  frequencies  (2.S  MHz),  the 
peak  optical  emission  has  moved  towards  higher  phase. 
These  results  indicate  that  the  arrival  time  of  the  ions  on  the 
electrode  and  the  subsequent  production  of  high-energy 
electrons  is  phase-  and  frequency-dependent  and  that  this 
can  have  an  influence  on  the  discharge  kinetics. 

The  energy  distributions  of  the  helium  and  argon  ions 
striking  the  electrodes  for  the  conditions  of  Fig.  1  are  shown 
in  Figs.  3  and  4.  Each  small  plot  in  Figs.  3  and  4  represents  a 
different  energy  range  (stacked  vertically)  and  rf  frequency 
(horizontal).  Within  each  small  plot  the  horizontal  axis  is 
phase  from  0-27r  and  the  vertical  axis  is  the  probability  of  an 
ion  arriving  at  that  phase  within  that  energy  range. 

As  previously  discussed,  for  low  frequencies,  the  ions 
track  the  rf  voltage.  Those  ions  that  arrive  during  the  ca¬ 
thodic  part  of  the  cycle  have  a  peak  in  the  ion  distribution  at 
a  phase  of  approximately  v/l,  independent  of  their  energy 
bin.  In  addition,  those  ions  which  diffuse  in  from  the  plasma 
or  are  left  outside  the  sheath  region  when  the  sheath  retreats 
to  the  electrode  are  collected  during  the  anodic  part  of  the 
cycle  with  low  energy,  since  the  electric  field  is  small.  As  the 
rf  frequency  increases,  both  the  high-  and  low-energy  ion 
distributions  tend  to  shift  toward  higher  phase.  However, 
the  higher  energy  ions  shift  towards  rr  and  subsequently  de¬ 
grade  to  lower  energy  since  the  ions  can  no  longer  follow  the 
field  and  acquire  the  maximum  sheath  potential.  Those  ions 
with  high  velocities  during  the  anodic  half  cycle  rapidly  lose 
their  energy  through  charge  exchange  collisions  and  popu¬ 
late  the  lower  energy  bins.  A  decrease  in  the  peak  ion  energy 
above  oi/  has  been  observed  in  chlorine  rf  discharges.  The 
ion  distribution  in  the  lowest  energy  bin  tends  to  spread  out 
and  become  more  uniform  over  the  entire  rf  cycle  as  a  conse¬ 
quence  of  the  ions  now  responding  to  a  quasitime  averaged 
sheath  voltage.  The  end  result  is  that  the  average  ion  energy 
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FIG.  7.  Average  ion  energy  as  a  function  of  rf  frequency  for  He  and  Ar  ions 
striking  the  electrodes.  Conditions  are  as  described  in  Fig.  I. 


slowly  decreases  with  increasing  frequency  until  a  frequency 
approximately  equal  to  a  few  when  the  ions  can  no  longer 
track  the  sheath  held.  At  that  point  the  average  energy  falls 
more  rapidly  (Fig.  7). 

The  effect  of  the  ion  mass  (or  o, )  can  be  observed  by 
comparing  Figs.  1  and  3  with  Figs.  2  and  4.  The  frequency 
when  the  ions  no  longer  follow  the  rf  frequency  is  much 
lower  in  argon  when  compared  to  helium.  In  both  gases  the 
ion  phase  distributions  follow  the  rf  cycle  at  low  frequency, 
with  the  most  probable  arrival  time  moving  towards  higher 
phase  and  spreading  out  with  increasing  rf  frequency.  How¬ 
ever,  while  the  frequencies  where  the  average  energy  drops 
off  and  the  ion  phase  distribution  becomes  more  uniform  in 
time  are  very  different  (Fig.  7),  they  nevertheless  scale  as 

CO,. 

B.  rf  voltage  and  dc  bias 

Increasing  the  rf  voltage  across  the  sheath  from  400  to 
8(X)  V  resulted  in  a  linear  increase  of  the  average  ion  energy 


NORMALIZED  ENERGY 

F10.  g.  Normalized  He  ion  energy  ai  a  function  of  rf  voltage  normalized  to 
400  V.  The  ion  energy  is  normalized  to  the  S8.4  eV  average  energy  obtained 
with  a  sheath  potential  of  400  V.  The  rf  frequency  is  2.0  MHz. 


Striking  the  electrode  (Fig.  8).  The  probability  distribution 
of  ion  arrival  phase  changed  only  slightly  with  increasing  rf 
voltage.  In  general,  there  is  a  small  increase  in  the  height  of 
the  ion  distribution  peak  (20%)  at  the  expense  of  the  ions 
that  would  normally  arrive  before  the  peak.  In  addition,  the 
absolute  position  of  the  peak  shifted  slightly  (less  than 
tr/lO)  towards  the  beginning  of  the  cycle.  Over  this  range  of 
rf  voltage,  only  the  average  energy  appears  to  significantly 
differ.  Thus,  the  linear  relation  between  average  energy  and 
sheath  voltage  appears  to  be  a  valid  and  useful  scaling  law. 

The  effect  of  imposing  a  dc  voltage  and  a  dc  sheath  on 
the  ion  phase  distribution  and  the  average  energy  was  also 
investigated.  For  these  calculations,  the  dc  voltage  and 
sheath  thickness  were  assumed  to  scale  linearly  with  the  rf 
voltage  and  sheath  thickness.  For  example:  with  =  400 
V  and  JTrf  =  0.5  cm.  if  =  200  V.  then  *  0.25  cm. 
Using  this  assumption,  the  ion  phase  distributions  were  cal¬ 
culated  for  He  ions,  as  shown  in  Fig.  9. 

When  the  dc  bias  and  dc  sheath  are  increased,  the  ion 
distribution  begins  to  take  on  the  characteristics  of  both  a  rf 
system  with  its  distinctly  peaked  behavior  and  a  dc  system, 
where  the  ion  flux  is  constant  in  time.  With  larger  bias,  the 
peak  becomes  smaller  and  more  ions  are  collected  during 
what  was  the  anodic  part  of  the  cycle.  However,  even  though 
the  distribution  becomes  more  uniform  with  increasing  dc 
bias,  the  average  energy  only  increases  by  approximately 
10%.  This  is  likely  a  result  of  the  competing  processes  of  an 
increased  sheath  voltage  raising  the  average  energy  while  an 
increased  sheath  thickness  decreases  the  average  energy  due 
to  an  increased  probability  for  charge  exchange  collisions 
(to  be  discussed  in  Sec.  Ill  C).  The  same  behavior  was  also 
observed  for  Ar  ions. 

C.  Sheath  thickneM  and  praaaura 

Increasing  the  rf  sheath  thickness  and  the  pressure  did 
not  have  a  large  effect  on  the  ion  phase  distribution.  With 
increasing  sheath  thickness  (0.5- 1.5  cm)  and  pressure 
( l.Ox  10'*-5.0x  10‘*  cm  "  ’ ),  the  peak  height  decreased  by 


no.  9.  Relative  total  arrival  probability  of  He  kme  ttriUni  the  electrode 
with  an  increaeing  dc  bias  and  dc  sheath  thicknew.  The  dc  sheath  thickneu 
is  assumed  to  increase  linearly  with  increasing  dc  bias.  The  rf  frequency  it 
2.0  MHz. 
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FIG.  10.  Normalized  He  ion  energy  as  a  function  of  the  reciprocal  of  the 
normalized  pressure  and  sheath  thickness.  The  ion  energy  is  normalized  to 
the  88.4  eV  average  ion  energy  in  a  2.0  MHz  discharge. 


15%  and  10%,  respectively,  with  a  resulting  increase  in  the 
number  of  ions  arriving  before  the  peak  of  the  distribution. 
In  both  cases  the  absolute  position  of  the  phase  peak  re¬ 
mained  in  approximately  the  same  position.  However,  as 
stated  previously,  the  average  ion  energy  increased  approxi¬ 
mately  linearly  with  respect  to  the  reciprocal  of  the  sheath 
thickness  and  pressure  (Fig.  10).  The  high-energy  tail 
dropped  off  quickly  with  increasing  sheath  thickness  and 
pressure  due  to  rn  increased  number  of  charge  exchange 
collisions. 

IV.  CONCLUSIONS 

Parametric  models  for  the  sheath  voltage  and  thickness 
and  a  Monte  Carlo  particle  simulation  were  used  to  compute 
the  arrival  phase  and  energy  of  ions  striking  the  electrodes  in 
low  pressure  rf  discharges.  When  the  rf  frequency  is  below 
the  ion  response  frequency,  the  ions  arrive  at  the  electrode  in 
phase  with  the  applied  rf  voltage.  As  the  rf  frequency  in¬ 
creases,  the  phase  of  the  highest  ion  arrival  probability  shifts 
towards  higher  phase  and  the  ion  energy  decreases.  The  opti¬ 
cal  emission  from  the  plasma  was  also  observed  to  peak  ap¬ 
proximately  in  phase  with  the  applied  rf  voltage  at  low  fre¬ 
quencies.  At  higher  frequencies,  the  peak  optical  emission 


moves  to  higher  phase.  The  transition  from  a  highly  peaked 
ion  distribution  at  low  frequencies  to  a  uniform  distribution 
at  higher  frequencies  required  at  least  an  order  of  magnitude 
change  in  rf  frequency.  Increasing  the  ion  mass  or  decreasing 
Oil  decreased  the  ion  energy  and  reduced  the  frequency  at 
which  the  phase  distribution  tends  towards  uniform.  In¬ 
creasing  the  rf  voltage  linearly  increased  the  average  ion  en¬ 
ergy,  while  the  sheath  thickness  and  pressure  are  approxi¬ 
mately  linearly  related  to  the  reciprocal  of  the  average 
energy.  Imposing  a  dc  bias  on  the  electrode  produced  a  su¬ 
perposition  of  a  peaked  ( if)  distribution  and  a  more  uniform 
(dc)  ion  distribution  with  little  change  in  average  ion  ener¬ 
gy.  Changing  the  rf  voltage,  sheath  thickness  and  pressure 
had  only  a  small  effect  on  the  shape  of  the  ion  phase  distribu¬ 
tion. 
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A  phenomenological  model  for  surface  deposition  kinetics  during  plasma 
and  sputter  deposition  of  amorphous  hydrogenated  silicon 
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The  surface  processes  during  the  plasma-enhanced  chemical  vapor  deposition  and  reactive 
sputter  deposition  of  amorphous  hydrogenated  silicon  (a-Si:H)  are  investigated  by  use  of  a 
phenomenological  model.  The  model  consists  of  an  accounting,  in  rate  equation  form,  of 
adsorption  of  radicals  from  the  plasma  onto  the  surface,  surface  diffusion,  incorporation  into 
the  lattice,  interconnection  of  bonds  in  the  lattice,  and  burial  of  species  on  the  surface,  thereby 
constituting  film  growth.  By  accounting  for  the  coordination  partners  of  Si  atoms  in  the  film, 
the  atomic  fraction  of  hydrogen  in  the  film  is  computed  for  the  lattice  and  for  hydrogen  in 
polymeric  or  isolated  configurations.  Results  from  the  model  are  discussed  while 
parametrizing  the  probability  for  hydrogen  elimination  during  incorporation  and  the 
probability  for  saturation  of  dangling  bonds  by  gas  phase  species.  We  find  that  the  mode  of 
hydrogen  elimination  during  incorporation  distinguishes  films  grown  dominantly  from  SiHj  or 
SiH]  radicals.  Characteristics  of  films  grown  by  sputter  deposition  are  investigated  as  a 
function  of  the  composition  of  the  radical  flux.  We  find  that  films  grown  from  hydrogen-rich 
fluxes  are  composed  dominantly  of  dihydride  (=Si  <{] )  configurations,  whereas  hydrogen- 
lean  mixtures  are  composed  of  dominantly  hydride  (sSi-H)  configurations. 


I.  INTRODUCTION 

Thin  films  of  amorphous  hydrogenated  silicon  (a-Si;H) 
are  produced  by  a  variety  of  means,  including  thermal  chem¬ 
ical  vapor  deposition  (CVD),'  homogeneous  chemical  va¬ 
por  deposition  ( HOMOCVD )  plasma-enhanced  chemical 
vapor  deposition  (PECVD),^  and  reactive  sputter  depo¬ 
sition  (RSD).^  The  films  thus  produced  typically  have 
atomic  hydrogen  fractions  of  5%-20%,  and  thicknesses  of 
3:0.5 /im  for  use  in  electronic  devices  and  10-15^m  for  use 
as  photoreceptor  material.’  Of  the  methods  of  fabrication 
listed  above,  the  former  two  rely  upon  thermal  pyrolysis  of 
gaseous  silane  compounds  (e.g.,  SiH4,  Si^H^)  to  generate  the 
radicals  from  which  the  film  is  grown;  the  latter  two  meth¬ 
ods  generate  the  radicals  by  electron  or  ion  impact  processes. 
The  purpose  of  this  paper  is  to  theoretically  investigate  the 
property  of  a-Si:H  films  produced  by  the  latter  two  methods. 

The  apparatus  typically  used  for  PECVD  is  a  parallel- 
plate  capacitively  coupled  radio-frequency  discharge  sus¬ 
tained  in  mixtures  of  silane  (or  disilane },  hydrogen,  and  no¬ 
ble  gases  at  pressures  <  0.5  Torr.*-’  The  substrate 
temperature  is  typically  400-700  K.  Silicon-containing  radi¬ 
cals  (e.g.,  SiHj,  SiHj)  are  generated  by  electron  impact  dis¬ 
sociation  of  the  feed  stock  at  discharge  power  levels  of  tens  to 
hundreds  of  mW  cm~’.  The  radicals  then  diffuse  out  of  the 
plasma  to  the  substrate  where  they  are  incorporated  into  the 
growing  film.  The  apparatus  typically  used  for  RSD  is  func¬ 
tionally  similar  to  that  for  PECVD,  except  that  one  of  the 
electrodes  (the  target)  is  covered  by  a  silicon  plate,  and  the 
gas  pressure  is  lower  (  <0. 1  Torr)."  Noble  gas  ions  from  the 
discharge  impact  into  the  target,  thereby  sputtering  Si 
atoms,  which  then  drift  to,  and  deposit  on,  the  substrate, 
typically  the  opposite  electrode.  By  flowing  controlled 
amounts  of  hydrogen  through  the  sputter  chamber,  the  frac¬ 


tion  of  hydrogen  incorporated  into  the  film  can  be  con¬ 
trolled.’  Additional  control  of  the  amount  of  incorporated 
hydrogen  in  the  film  is  obtained  in  both  methods  by  regulat¬ 
ing  the  substrate  temperature;  the  higher  the  substrate  tem¬ 
perature,  the  lower  the  hydrogen  fraction. 

The  processes  through  which  radicals  generated  by  the 
plasma  during  PECVD  and  RSD  are  adsorbed  onto  the  sur¬ 
face  and  incorporated  into  the  film  are  complex,  and  the 
study  of  those  processes  are  current  topics  of  research. 

The  processes  can  be  divided  into  the  following  general  cate¬ 
gories:  adsorption,  incorporation,  hydrogen  elimination,  in¬ 
terconnection,  etching,  sputtering,  and  burial.  Adsorption  is 
the  process  whereby  radicals  from  the  plasma  stick  to  the 
surface  and  become  available  for  bonding  to  the  film.  Incor¬ 
poration  is  the  process  whereby  a  radical  on  the  surface  co¬ 
valently  bonds  to  another  silicon  atom  in  the  film.  We  call 
this  bonding  to  the  lattice,  although  the  use  of  the  term  lat¬ 
tice  does  not  claim  any  long-term  order.  Incorporation  can 
be  direct,  in  which  a  radical  attaches  to  the  lattice  at  the  site 
of  a  dangling  bond;  or  indirect,  in  which  the  radical  displaces 
hydrogen  from  saturated  bonds  (i.e.,  Si-H),  thereby  elimi¬ 
nating  hydrogen  from  the  surface.  Etching  is  the  process 
whereby  radicals  remove  atoms  from  the  lattice  or  remove 
adsorbed  species  by  colliding  and  bonding  to  those  species 
but  not  sticking  on  the  surface.  Sputtering  is  the  analogous 
process  performed  by  energetic  ions  colliding  with  the  sur¬ 
face.  Burial  is  the  process  during  which  interconnection  and 
incorporation  cover  over  atoms  previously  on  the  surface 
and  e.xposed  to  the  plasma,  and  confine  them  to  the  lattice. 

The  fraction  of  hydrogen  in  a-Si:H  films,  5%-20%,  is 
small  compared  to  the  ratio  of  hydrogen  to  silicon  atoms  in 
the  flux  of  radicals  which  adsorb  or  collide  with  the  growing 
film.  The  manner  in  which  the  excess  hydrogen  is  eliminated 
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is  therefore  an  important  consideration.  Heuristic  models 
for  growth  of  a-Si:H  films  by  PECVD  have  been  suggested 
by  Kampas,’  Scott,  Reimer,  and  Longeway,"  Longeway,'"* 
Gallagher, and  others.  These  models  differ  qualitatively  in 
the  manner  in  which  hydrogen  is  eliminated  from  the  film. 
Kampas*  suggests  that  cross  linking  between  adjacent  Si-H 
bonds  on  the  surface  resulting  in  elimination  and  desorption 
of  is  the  dominant  mechanism.  Scott"  also  suggests  that 
hydrogen  atoms  are  removed  from  Si-H  bonds  on  the  sur¬ 
face  via  etching  reactions  by  radicals  from  the  plasma.  An¬ 
other  suggested  method  is  that  radicals  (e.g.,  SiHj)  bond  to 
the  lattice  as  activated  silyl  groups,  which  spontaneously 
eliminate  Hj  as  a  method  of  relaxation.  Gleason  et  al.  “  have 
recently  published  a  Monto  Carlo  simulation  for  thin-film 
growth  of  a-Si:H.  In  their  model  the  initiating  step  (called 
the  addition  reaction)  proceeds  by  incorporation  of  SiH4 
directly  into  the  film  with  the  simultaneous  elimination  of 
H^.  This  step  resembles  heterogeneous  pyrolysis,  and  there¬ 
fore  the  model  is  most  directly  applicable  to  CVD  (as  op¬ 
posed  to  PECVD).  In  their  model,  though,  hydrogen  is  eli¬ 
minated  from  the  film  by  cross  linking  between  adjacent 
^i-H  bonds,  in  a  manner  similar  to  that  suggested  by 
Kampas.’ 

In  this  paper  a  phenomenological  model  for  the  surface 
deposition  kinetics  during  PECVD  and  RSD  of  a-Si;H  is 
presented,  and  results  from  the  model  are  discussed  for  a 
variety  of  plasma  conditions.  The  model  consists  of  a  set  of 
rate  equations  describing  the  major  processes  of  adsorption, 
surface  diffusion,  incorporation,  interconnection,  hydrogen 
elimination,  and  burial.  From  the  results  of  the  model,  one 
obtains  the  fraction  of  hydrogen  in  the  film,  and  the  distribu¬ 
tion  of  hydride  (^i-H)  and  dihydride  (=Si  <{} )  bonds  in 
the  film.  The  model  is  described  in  Sec.  II  after  which  rate 
constants  used  in  the  model  are  discussed  in  Sec.  III.  Results 
from  the  model  are  presented  in  Sec.  IV  in  the  subtopics  of 
hydrogen  elimination  during  incorporation,  the  burial  ef¬ 
fect,  saturation  and  etching,  temperature  effects,  reactive 
sputter  deposition,  and  deposition  in  the  absence  of  SiH, 
insertion.  Concluding  remarks  are  in  Sec.  V. 

II.  DESCRIPTION  OF  THE  MODEL 

A  model  for  gas  phase  chemical  kinetics  can  be  formu¬ 
lated  by  defining  a  set  of  species  and  rates,  and  expressing 
them  as  a  set  of  rate  equations.  The  rate  equations  have  the 
general  form 

^  =  (I) 

at  rS  l.m 

where  k  is  the  rate  constant  for  formation  of  species  N,  by 
collisions  between  species  N,  and  N„.  The  terms  on  the 
right-hand  side  of  Eq.  ( 1 )  simply  represent  the  sources  and 
sinks  for  species  N,.  The  surface  deposition  model  (SDM) 
described  in  this  paper  consists  of  a  set  of  rate  equations 
similar  to  that  which  one  would  write  for  gas  phase  chemis¬ 
try.  There  are  six  classes  of  species  in  the  model.  Allowing  X 
to  denote  some  specific  species,  the  six  classes  are  the  follow¬ 
ing:  X(F),  the  flux  of  radicals  incident  onto  the  surface; 
X(A),  adsorbed  nonmobile  radicals;  X{M),  adsorbed  mo¬ 
bile  radicals;  X(G),  desorbed  gas  phase  molecules;  X(.L),  a 


silicon  atom  bounded  to  the  lattice  in  the  film  but  residing  on 
the  surface;  and  X(B),i  silicon  atom  bonded  to  the  lattice  in 
the  film  and  buried  beneath  the  surface. 

The  individual  species  in  each  class  described  above  are 
listed  in  Table  I.  The  species  X(L)  and  X(B)  represent  sili¬ 
con  atoms  on  and  below  the  surface.  The  relative  densities  of 
species  in  class  X(.B)  contain  the  information  from  which 
the  fractions  of  hydrogen  and  silicon  in  the  film  are  obtained 
and  the  types  of  bonds  in  the  film  are  calculated.  This  infor¬ 
mation  is  obtained  by  tagging  the  silicon  atoms  X(L)  and 
X(B)  with  their  coordination  partners.  For  example, 
SSSH(R)  represents  a  silicon  atom  that  is  part  of  the  lattice 
and  buried  beneath  the  surface.  It  is  bonded  to  three  other 
silicon  atoms  (SSS)  and  a  hydrogen  atom  (H).  SSHD(Z.) 
represents  a  silicon  atom  on  the  surface  bonded  to  two  sili¬ 
con  atoms  (SS),  one  hydrogen  atom  (H)  and  having  one 
dangling  bond  (D).  The  atomic  fractions  of  hydrogen  and 
silicon  in  the  film  are  obtained  by  summing  the  silicon  and 
hydrogen  bonds  of  the  buried  silicon  atoms  and  taking  the 
appropriate  ratios. 

The  following  conceptual  sequence  of  events  is  simulat¬ 
ed  by  the  SDM.  Radicals  are  generated  in  the  gas  phase  by 
electron  impact  and  molecular  reactions,  and  diffuse 
towards  the  surface.  The  composition  of  the  radical  flux  is 
either  specified  as  part  of  a  parametric  study  or  obtained 
from  the  results  of  a  companion  plasma  chemistry  model  for 
rf  discharges  in  silane  gas  mixtures."  The  radicals  then  ad¬ 
sorb  on  the  surface  as  mobile  adsorbates  [X{.M)  ] ,  or  etch  or 
saturate  the  surface  (see  below).  As  energetically  allowed, 
sputtering  by  ions  directed  towards  the  surface  may  also  oc¬ 
cur.  Adsorbed  mobile  radicals  diffuse  on  the  surface.  These 
radicals  thermally  desorb  from  an  immobile  state  and  read¬ 
sorb  to  the  immobile  state  desorbing  many  times  before  col¬ 
liding  with  a  bonding  site  in  the  lattice  on  the  surface  [e.g., 
SSSD(L) ].  At  that  time  it  is  incorporated  into  the  lattice, 
and  hydrogen  may  be  eliminated  from  the  surface  in  the 
manner  described  below.  Mobile  radicals  which  collide  and 
react  with  other  mobile  radicals  may  desorb  from  the 
surface  as  saturated  molecules  [e.g.,  SiH](^)  SiHjlAf) 
—  SijHgl  (7)  ].  Heterogeneous  production  of  saturated  mole¬ 
cules  in  this  fashion  is  believed  to  be  an  important  contribu¬ 
tion  to  the  density  of  gas  phase  species.  Once  a  member 
of  the  lattice,  further  reactions,  such  as  interconnection,  may 
take  place  between  adjacent  silicon  atoms.  The  silicon  atoms 
on  the  surface  are  gradually  buried  [i.e.,  SSSH(L) 
— SSSH(R)  ]  at  a  rate  determined  by  the  magnitude  of  the 


TABLE  I.  Species  included  in  the  surface  deposition  model. 


Surface  and  lattice 
species  IX{L).  X{B)  ] 

Radical  flux,  adsorbed, 
and  buried  isolated  species 
{XiF),XlM),XlA),X(B)] 

SSSS 

H  Si  SiH  SiH,  SiH, 

SSSD  SSSH 

Si,H,  Si,H,  Si,H,  Si,H, 

SSDD  SSHD  SSHH 

SDDD  SHOD  SHHD  SHHH 

Saturated  gas  phase  species  [XiO  ] 

H,  SiH,  Si.H, 
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flux  of  radicals  incident  onto  the  surface  and  by  the  require¬ 
ment  that  the  total  density  of  surface  species  remains  a  con¬ 
stant.  Burial  is  distinct  from  interconnection.  Radicals  on 
the  surface  may  interconnect;  however,  they  are  not  buried 
beneath  the  surface  unless  there  is  a  flux  of  radicals  to  cover 
them. 

There  are  two  types  of  surface  sites.  The  first  site  is  the 
location  at  which  a  silicon  atom  bonds  to  the  lattice;  the 
second  site  is  a  location  at  which  a  radical  is  immobily  ad¬ 
sorbed.  The  density  of  bonding  sites,  is  a  ~  where  a  is  the 
average  bond  length,  taken  to  be  2.7S  A.  The  density  of  ad¬ 
sorption  sites  was  assumed  to  be  O.Spa-  The  distinction 
between  mobile  and  immobile  adsorbed  radicals  is  that  mo¬ 
bile  radicals  are  free  to  diffuse  along  the  surface.  Immobile 
adsorbed  radicals  are  confined  to  their  adsorptions  sites,  but 
they  may  thermally  desorb  and  become  mobile.  Mobile  and 
immobile  adsorbed  radicals  were  invoked  to  account  for  the 
observation  that  hydrogen  in  u-Si:H  films  produced  by 
PECVD  is  contained  in  both  the  lattice  and  in  isolated  poly¬ 
meric  configurations,  (SiH,),,  We  attribute  the  isolated 
configurations  to  clusters  of  adsorbed  SiH„  molecules  which 
the  film  has  grown  over  and  buried,  but  which  are  not  incor¬ 
porated  (i.e.,  bonded)  to  the  lattice. 

Rate  equations  were  formulated  for  the  processes  de¬ 
scribed  above  using  the  coordination  partner  notation  to  ac¬ 
count  for  bond  densities.  For  example,  the  interconnection 
of  dangling  bonds  of  two  adjacent  silicon  atoms  on  the  sur¬ 
face  is  represented  schematically  as 

SSHDd)  -f-  SSSDfZ.)  -SSSHd)  -f  SSSS(Z,).  (2) 

The  contribution  of  this  process  to  the  rate  equations  is  sim¬ 
ply  (SSHD(£)l(SSSD(L)lk,  where  (A"]  is  the  surface  den- 
sity  of  the  species  and  k  is  the  reaction  rate  constant.  A  set  of 
sample  reactions  appear  in  Table  II.  All  energetically  al¬ 
lowed  combinations  of  reactants  from  Table  I  for  the  sample 
reactions  were  included  in  the  model. 

III.  REACTION  RATE  CONSTANTS 

Reaction  rate  constants  are  required  for  all  the  pro¬ 
cesses  described  above;  adsorption,  desorption,  surface  dif¬ 
fusion,  lattice  building,  interconnection,  and  burial.  Al¬ 
though  some  of  these  rates  are  known  for  crystalline  silicon, 
such  as  the  adsorption  probability,'"  they  are  generally  not 


known  for  amorphous  silicon,  particularly  when  dangling 
bonds  on  the  surface  may  be  passivated  by  hydrogen.  There¬ 
fore,  the  rate  constants  must  be  estimated  in  an  internally 
consistent  fashion,  as  described  below. 

A.  Adsorption,  passivation,  and  etching 

This  category  of  reactions  includes  processes  involving 
the  flux  of  radicals  or  molecules  from  the  plasma.  These 
reactions  are  represented  schematically,  or  by  example,  by 

X(F)—X{.M)  (adsorption), 

(3) 

-I- SSSD(Z,)  —  SSSH(Z,)  -\-X{F)  (passivation), 

(4) 

X{F)  -|-SSSH(I,)-.SSSD(Z.)  -l-;irH(/0  (etching). 

(5) 

Since  the  values  of  X(F)  in  the  SDM  are  specified  param¬ 
eters  and  are  nearly  independent  of  the  surface  processes 
(see  below),  their  densities  are  constants  in  the  expressions 
above.  The  time  rate  of  change  of  the  density  of  a  surface 
species,  X(L),  expressed  as  cm  ~  ^  s  ~  ‘ ,  for  process  j  involving 
A’(F)  is  therefore  [X{P)]p^[X(L)  \/po,  where Py  is  the  prob¬ 
ability  of  the  process  and  [X{L)]  /py  is  the  fractional  surface 
coverage  of  the  reactant  surface  species.  For  adsorption,  the 
fractional  surface  coverage  is  set  to  unity  so  that  the  rate  of 
adsorption  depends  only  on  the  probability  Pj .  We  treat p^  as 
a  parameter  in  the  model. 

B.  Mobile  and  immobile  adsorbed  radicals 

Mobile  radicals  X(M)  are  immobilized  by  “colliding" 
with  an  active  adsorption  site.  Since  the  adsorption  site  can 
be  considered  equivalent  to  a  lattice  site,  the  rate  of  immobile 
adsorption  is  as  given  in  the  next  section.  The  rate  of  desorp¬ 
tion  for  X{.A )  -‘X(Af)  is  v^e  ~  where  is  the  de¬ 

sorption  frequency,  T,  is  the  substrate  temperature,  and  Af „ 
is  the  adsorption  energy.  By  parametrizing  the  results  of  the 
model  for  the  hydrogen  fraction  in  isolated  configurations  as 
a  function  of  surface  temperature  (see  below),  and  compar¬ 
ing  to  experiment,’  we  estimated  that  Vj  =  10’  s~  ‘  and  Af^ 
=  0.2  eV. 


TABLE  II.  Examplary  reactions  included  in  the  model  (see  text  Tor  rales).* 


SiH,(f)-SiH,(/Vr) 

SiH,(.W)=:SiH,(/<) 

SiH,(/^  H-SSSH(L)-SiH,(n  -)-SSSD(L) 

M  *(n  -f-SSSH(/.)-.W(f)  +SSSD(L)  +  H{C) 
SiH,(n  +SSSD(Z.)-SiH,(n  -t-SSSHCD 
SiH,(.Vf)  +  SiH,(W)-Si,H„(C) 

SiH.(.ir)  +  SSSD(L)  -SSSS(L)  +  SSHH(L) 

SiH.C.tT)  +SSSH(L)-SSSS(L)  +SHHD(t)  +  8,(0 
SSSD(Z.)  +  SSSD(i)-SSSS(L)  +SSSS(/.) 

SSSH(L)  +SSSH(L)-SSSS(Z.)  +SSSS(£)  +  H,(C) 
SSSH(Z.)-SSSH(«) 

SiH2(^)-SiH,(fi) 


Adsorption  to  mobile  slate 

Mobile  adsorption^Immobile  desorption 

Etching 

Sputtering'’ 

Saturation 

Surface  clustering  and  desorption 
Incorporation 

Incorporation  with  H,  elimination 
Interconnection 

Interconnection  with  H-  elimination 
Burial  (lattice  species) 

Burial  (isolated  configurations) 


*  All  allowed  permutations  of  exemplary  reactions  for  the  species  listed  in  Table  I  are  included  in  the  model. 
"  denotes  any  energetic  ion  incident  onto  the  surface. 


4765 


J.  Appl.  Phys.,  Vol.  62,  No  12. 1 5  December  1987 


Mark  J.  Kushner 


4765 


C.  Surface  clustering,  lattice  building,  and 
interconnection 

The  reactions  in  this  category  have  the  general  or  exam- 
plary  form  of 

X(M)  +  Y{M)  -‘XY(M)  (surface  clustering), (6) 

X(M) Y(L)—XY(L)  (lattice building),  (7) 

SSSD(L)  +SSHD(I,) 

—  SSSS(Z.)  +  SSSH(Z.)  (interconnection).  (8) 

The  time  rate  of  change  of  a  surface  species  resulting  from 
process  j  in  this  category  is  [X][Y]kjPje~  where 

[X]  and  ( Y]  are  the  densities  of  the  reactants  (cm”*),  kj  is 
the  rate  constant  (cm*  s”'),  Afy  is  an  activation  energy,  and 
Pj  is  a  relative  probability  of  the  process.  For  processes  in¬ 
volving  reactants  that  are  both  lattice  species,  Pj  =  \.  For 
processing  where  one  of  the  reactions  is  a  mobile  species,  and 
therefore  diffusing  on  the  surface,  pj  =  {{T,/5QO  K) 
X{M /{2p.)]y'-,  where  M  is  the  molecular  weight  of 
SiH](iVf)  and  p.  is  the  reduced  mass  of  the  reactants.  This 
scaling  parameter  merely  accounts  for  lighter  adsorbates 
having  a  higher  thermal  velocity  and  therefore  a  higher 
probability  for  reaction  per  unit  time.  Unless  otherwise  not¬ 
ed,  for  reactions  between  adsorbed  species,  or  adsorbed  spe¬ 
cies  and  the  lattice,  Aiy  (T,  =  500  K)  =  1.5 X  10"'*cm*s”' 
and  Af  =  0  (see  below).  For  interconnection  reactions 
between  saturated  lattice  species  [i.e.,  SSSH(L)],  kj 
=  3.0x  10”“  cm*  s" ',  and  Af  is  either  1.5  or  10  kcal/mol 
(see  below). 

For  unity  coverage  of  reactants,  the  reaction  frequency 
is  approximately  2000  s”'.  The  average  time  a  radical  is 
immobily  adsorbed  (T,  =  500  K)  is  0.75  ms.  Therefore, 
radicals  spend  little  time  in  the  mobile  state,  but  rather  move 
from  immobile  site  to  immobile  site  or  to  an  incorporation 
site.  The  results  of  the  model  are  relatively  insensitive  to  the 
choice  of  kj  provided  that  the  reaction  frequency  is  compara¬ 
ble  to  or  larger  than  the  desorption  frequency.  The  values  for 
rate  constants  were  obtained  by  parametrically  comparing 
the  results  of  this  model  to  the  experimental  results  of  Ross 
and  Jaklik.*  The  values  derived  are  not  unique  and  are  mere¬ 
ly  intended  to  represent  general  reaction  probabilities. 

The  activation  energies  Af  were  estimated  for  each  reac¬ 
tion  by  assigning  bond  energies  and  computing  the  change  in 
total  energy  going  from  reactants  to  products.  Exothermic 
reactions  were  assigned  zero  activation  energies;  endother¬ 
mic  reactions  were  assigned  an  activation  energy  equal  to  the 
endothermicity.  The  bond  energies  for  surface  configura¬ 
tions  were  estimated  from  the  work  of  Ho  et  al.  ”  who  calcu¬ 
lated  a  correction  factor  for  Si-Si  and  Si-H  bond  energies  in 
higher  silane  molecules.  In  general,  all  clustering,  insertion, 
and  lattice  building  reactions  involving  species  having  dan¬ 
gling  bonds  are  exothermic  to  the  degree  of  approximation  of 
this  method.  The  endothermic  reactions  have  energies  on  the 
order  of  a  few  to  10  kcal/mol.  These  reactions  are  primarily 
lattice  interconnection  with  molecular  hydrogen  elimina¬ 
tion.  The  activation  energy  associated  with  these  reactions  is 
responsible,  in  part,  for  the  observed  reduction  in  the  /h 
with  increasing  substrate  temperature.*  "  *®  *'  We  assigned 
an  activation  energy  of  1.5  kcal/mol  for  interconnection  re¬ 


actions  involving  the  lattice  species  SSSH(L);  and  10  kcal/ 
mol  for  reactions  involving  SSHH (/I )  and  SHHH (Z. ) .  This 
difference  in  activation  energies  will  be  discussed  in  Sec. 
IV  D. 

The  rate  of  interconnection  between  dangling  bonds 
Ie.g.,SSSD(£)  -fSSSD(Z,)-.SSSS(I,)  -I- SSSS(L)]  is  cal¬ 
culated  during  execution  of  the  model.  The  instantaneous 
rate  is  chosen  so  that  upon  burial,  the  fractional  density  of 
dangling  bonds  is  <  1  %.  The  rate  so  derived  varies  from  case 
to  case,  but  is  generally  close  to  2-5 X  10”'*  cm*  s”'.  Al¬ 
though  the  spin  density  in  a-Si:H  is  usually  much  less  than 
0.01,  we  found  that  requiring  a  smaller  density  than  this 
upon  burial,  that  is,  having  a  higher  rate  of  interconnection, 
resulted  in  unphysically  small  hydrogen  fractions.  The  im¬ 
plication  of  this  observation  is  that  interconnection  between 
dangling  bonds  continues  in  the  film  after  burial. ' ' 

C.  Burial 

Burial  is  the  process  whereby  species  residing  on  the 
surface  are  incorporated  into  the  film.  For  lattice  species, 
this  is  schematically  shown,  for  example,  as 

SSSH(Z.)-SSSH(fl).  (9) 

The  rate  of  burial  for  surface  specie  Xj{L)  is  (Xj(L)/po) 
X  (Pl  —  for  pL  >Po>  zero  otherwise,  where  pt_ 
=  2j  Xj(L)  is  the  sum  of  all  lattice  species  and  r  is  the  burial 
equilibration  time.  This  rate  “clamps”  the  density  of  surface 
lattice  species  at  po-  model  is  not  sensitive  to  the  choice 

of  r  as  long  as  r  <po/ [  Ij  njXj  (/^  ] ,  where  Xj(F)  is  the  net 
flux  of  radical  specie  j  that  sticks  on  the  surface  and  n,  is  the 
number  of  silicon  atoms  contained  in  the  radical.  For  very 
high  deposition  rates  (  >  10(X)  A/min )  the  limiting  value  for 
r  is  =;  150  ms;  therefore,  we  used  r  =  10  ms  for  all  cases. 

D.  Initial  conditions 

The  initial  conditions  used  for  the  model  are  that  the 
surface  is  covered  with  dangling  bonds  [i  e.,SSSD(L)]  hav¬ 
ing  density  po>  ^nd  unoccupied  adsorption  sites  having  den¬ 
sity  O.Spy.  The  fluxes  are  turned  on  and  the  rate  equations  are 
integrated  until  steady-state  conditions  are  reached.  We 
found  that  for  constant  fluxes,  final  surface  conditions  were 
obtained  in  1-3  s  of  real  time,  corresponding  to  only  5-10 
“monolayers”  of  film,  and  that  the  hydrogen  content  de¬ 
creased  to  its  steady-state  value.  These  observations  agree 
well  with  the  results  of  Gleason  et  a/.,'®  who  calculated  that 
steady-state  conditions  are  obtained  after  depositing  a  simi¬ 
lar  number  of  monolayers,  and  hydrogen  content  decreases 
towards  its  steady-state  value.  We  have  also  started  with  a 
prepared  surface  [e.g.,  a  partial  covering  of  SSSH(Z,)  ];  and 
found  little  or  no  change  in  the  steady-state  results  when 
compared  to  starting  with  a  fresh  surface. 

IV.  RESULTS 

In  the  plasmas  of  interest,  the  dominant  radicals  are 
thought  to  be  SiH,  and  SiH,.**  **  The  ratio  of  hydrogen 
atoms  to  silicon  atoms  in  silicon-containing  radicals  incident 
onto  the  surface,  ( H  ]  /  [  Si  ]  =  F„ ,  is  therefore  approximate¬ 
ly  2-3.  The  eHctron  impact  dissociation  of  saturated  silanes 
which  generated  the  radicals  also  generate  a  nearly  equal 


4766 


J.  Appl.  Phys..  Vol.  62.  No.  12. 15  December  1987 


Mark  J.  Kushner 


4766 


density  of  atomic  hydrogen.  The  flux  of  hydrogen  atoms 
incident  onto  the  surface  is  therefore  larger  (6-8)  than  the 
value  given  by  the  silicon-containing  radicals  alone.  Addi¬ 
tional  hydrogen  is  deposited  on  the  surface  by  hydrogen- 
containing  molecules  from  the  gas  phase  which  saturate  dan¬ 
gling  bonds,  as  discussed  below.  The  fraction  of  hydrogen 
incorporated  into  the  film,  /„,  though,  is  usually  <0.20. 
Therefore,  hydrogen  must  be  eliminated  from  the  surface 
some  time  between  the  processes  of  adsorption  and  burial;  or 
once  buried,  hydrogen  must  be  eliminated  from  the  lattice. 
In  this  work  we  e.xamine  the  former  process.  We  assumed 
that  hydrogen,  once  buried  beneath  the  surface,  is  perma¬ 
nently  part  of  the  film. 

In  the  results  discussed  below,  we  either  specify  the  radi¬ 
cal  flux  as  part  of  a  parametric  survey  or  use  the  flux  comput¬ 
ed  with  the  plasma  chemistry  model.  By  specifying  a  flux,  we 
ignore  the  effect  of  surface  chemistry  on  the  radical  flux. 
Since  the  majority  of  the  desorbed  products  of  the  surface 
chemistry  are  saturated  molecules,  this  effect  is  small. 

A.  Hydrogen  elimination  during  incorporation 

When  adsorbed  molecules  are  incorporated  into  the  lat¬ 
tice,  they  attach  to  either  dangling  bonds  [e.g.,SSSD(I.)  ]  or 
insert  into  saturated  bonds  (e.g.,  SSSH(£.)  ].  Any  adsorbed 
radical  may  incorporate  directly  into  a  site  occupied 
by  a  dangling  (e.g.,  SiHj(A/) -|- SSSD(Z.)  — SSSS(Z.) 
+  SHHH(Z,)  ].  Incorporation  by  insertion  into  a  saturated 
bond,  though,  may  necessarily  require  elimination  of  hydro¬ 
gen.  For  purposes  of  identification,  we  label  the  adsorbed 
radical  by  its  ratio  of  hydrogen  to  silicon  atoms,  =  ( H )/ 
(Si).  For  Fh  <2,  adsorbed  radicals  may  insert  directly  into 
saturated  SSSH(Z.)  bonds  without  eliminating  hydrogen; 
for  example, 

SiHjiM)  +SSSH(£)-SHHH(L)  -(-SSSS(£).  (10) 

For  adsorbed  radicals  having  F„  >  2,  the  insertion  reaction 
requires  elimination  or  displacement  of  hydrogen.  Elimina¬ 
tion  may  proceed  by  displacing  atomic  or  molecular  hydro¬ 
gen.  This  method  of  incorporation  of  radicals  having  Fh  >  2 
is  the  analogy  of  the  gas  phase  reaction  SiHj  +  SiH4 
— SijHs  +  H,  which,  at  typical  substrate  temperatures,  pro¬ 
ceeds  with  a  rate  constant  of  10"  cm’ s"  If  radicals 
having  a  single  dangling  bond  [e.g.,  SiH](A/)  ]  are  unable  to 
insert  into  a  saturated  surface  site  [e.g.,  SSSH(L)],  then 
incorporation  can  only  proceed  by  generating  dangling 
bonds  on  the  surface  by  sputtering  or  etching  mechanisms.  ” 
In  this  section  we  investigate  surface  deposition  while  allow¬ 
ing  insertion  of  SiH3(A/)  into  saturated  surface  sites.  In  Sec. 
V  F  we  examine  some  of  the  consequences  of  requiring  a 
dangling  bond  for  incorporation  of  SiHj  (A/)  into  the  lattice. 

For  the  incorporation  of  SiHj(A/)  while  allowing  inser¬ 
tion  into  saturated  surface  sites,  we  may  have  either  of’  ” 

SiHjl^Vf)  -|-SSSH(I.) 

-SHHD(Z.) -)-SSSS(Z.)  4- H,(C)  (II) 

-SHHHd) -f  SSSS(Z,) -h  HlC).  (12) 

The  former  reaction  is  thermodynamically  favored.  If  it  is 
more  likely  to  occur,  the  elimination  of  molecu  .ir  hydrogen 
during  incorporation  of  SiHjfiV/)  into  films  grown  from 


plasmas  whose  predominant  radical  isSiH}!/^  may  result  in 
less  hydrogen  remaining  on  the  surface  than  incorporation 
of  S’H}!^/),  a  radical  having  a  smaller  value  of  Fh  . 

Before  examining  the  results  of  the  model  for  actual  dis¬ 
charge  conditions,  it  is  instructive  to  examine  the  character¬ 
istics  of  o-Si:H  films  grown  from  specific  fluxes  of  radicals. 
We  will  examine  the  fraction  of  incorporated  hydrogen,  /h  , 
as  a  function  of  the  hydrogen-to-silicon  ratio  Fh  of  the  indi¬ 
vidual  radicals  incident  on  to  the  surface.  This  is  accom¬ 
plished  by  specifying  that  the  radical  flux  is  composed  of 
only  a  single  specie  having  a  given  F'h.  These  results  are 
shown  in  Fig.  1  where  /h  is  plotted  as  a  function  of  Fh  while 
allowing  adsorbed  radicals  with  a  single  dangling  bond  to 
insert  into  a  saturated  surface  site.  The  flux  of  silicon  atoms 
incident  on  the  surface  was  kept  constant  at  2  X  10' Vcm’ s, 
typical  of  moderately  power  discharges  of  tens  to  100  mW/ 
cm’  sustained  in  a  few  hundred  mTorr  of  SiHj  and  corre¬ 
sponding  to  a  deposition  rate  of  a-Si:H  of  =:250  A/min.  As 
expected,  /h  increases  with  increasing  F„,  but  only  to  a 
ratio  of  Fh  =:  2  ( as  shown  by  the  triangles  in  the  figure ) ,  after 
which  /h  decreases.  This  trend  is  a  result  of  having  (as 
opposed  to  H)  being  eliminated  during  incorporation  of  ad¬ 
sorbed  radicals,  that  is,  using  the  insertion  step  shown  in  Eq. 
(11)  as  opposed  to  that  shown  in  Eq.  (12).  If  we  instead 
specify  that  H  atoms  are  eliminated  during  incorporation 
(dashed  line  in  Fig.  1),  /h  differs  little  from  the  cases  for 
elimination  of  for  Fh  <2,  but  is  larger  than  the  previous 
cases  for  Fh  >2.  For  these  conditions,  /„  is  equally  deter¬ 
mined  by  hydrogen  elimination  during  interconnection  of 
adjacent  Si-H(Z.)  bonds  after  incorporation,  and  hydrogen 
elimination  during  incorporation.  The  latter  process  can 
have  a  significant  effect  on  /h  for  plasmas  dominated  by 
radicals  having  Fh  >  2. 

The  sensitivity  of  /h  on  the  mode  of  hydrogen  elimina¬ 
tion  during  incorporation  is  further  illustrated  in  Fig.  2 
where  /„  is  plotted  for  films  grown  from  SiHjlF)  and 
SiH)(F)  as  a  function  of  the  probability  of  eliminating  Hj. 


FtG.  1.  Atomic  rraction  of  hydrogen  in  a-Si:H  films,  /„,  as  computed  with 
the  surface  deposition  model  as  a  function  of  the  H/Si  ratio  of  individual 
radicals  incident  onto  the  surface.  For  these  results,  adsorbed  radicals  hav¬ 
ing  a  single  dangling  bond  may  insert  into  saturated  surface  sites.  The  trian¬ 
gles  are  for  H,  elimination  during  incorporation;  the  dashed  line  is  for  H 
atom  elimination  during  incorporation. 
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FIG.  2.  Total  hydrogen  fraction  (lattice  and  isolated)  as  a  function  of  the 
probability  of  eliminating  or  H  for  incident  fluxes  of  SiH;(F)  and 
SiH,(F).  A  high  probability  for  eliminating  Hj  distinguishes  between  a- 
Si:H  films  grown  from  SiH,(F)  and  SiH,(F)  when  SiHjCW)  insert  into 
saturated  surface  sites. 


The  value  of  /„  for  films  grown  from  is  relatively 

insensitive  to  this  probability,  and  when  H  atoms  are  elimin¬ 
ated  there  is  little  distinction  between  films  grown  from 
SiH;(/^  and  SiMjC/").  As  the  probability  of  H,  elimination 
increases,  though,  /„  decreases  for  films  grown  from 
SiHjC/-). 

Computed  atomic  hydrogen  fractions  in  a-Si:H  films  for 
films  grown  in  Ar/SiH^  rf  discharges  are  shown  in  Fig.  3. 
The  discharge  conditions  were  simulated  with  the  plasma 
chemistry  model  described  in  Ref.  17.  The  plasma  condi¬ 
tions  for  these  examples  are  as  follows:  gas  pressure,  0.5 
Torr;  substrate  temperature,  500  K;  and  power  deposition, 
=:  25  mW  cm  ~  Cases  are  shown  for  varying  the  fraction  of 
SiHj  in  Ar.  To  eliminate  the  burial  effect  (see  next  section) 
the  power  for  each  case  was  varied  to  obtain  the  same  depo¬ 
sition  rate.  Since  the  total  atomic  hydrogen  fraction  in¬ 
creases  with  decreasing  fraction  of  SiH^  while  the  deposition 
rate  is  constant,  the  change  in  /„  is  a  result  of  a  change  in  the 
composition  of  the  radical  flux.  Also  plotted  in  Fig.  3  are  the 
ratios  of  the  total  flux  of  H  atoms  to  Si  atoms  incident  on  the 
surface,  [  H/  ( F)  ]  /  ( Si/  ( F)  ] ;  and  the  ratio  of  radicals  which 
have  Fh<2  to  those  having  Fh  >2,  [Si,„H,,  n/m<;2)/ 
[Si„H„,  n/m  >  2].  Both  ratios  increase  with  decreasing 
fraction  of  silane.  This  trend  is  a  consequence  of  silane  radi¬ 
cals  and  hydrogen  atoms,  once  generated  by  electron  impact 
of  silane,  having  a  smaller  probability  of  colliding  with  other 
silane  molecules  before  striking  the  substrate  as  the  silane 
fraction  decreases.  Of  the  two  ratios,  the  latter,  that  compar¬ 
ing  radicals  having  Fh  greater  than  and  less  than  2,  corre¬ 
lates  better  with  /h  •  The  increase  in  /„  may  therefore  be 
caused,  in  part,  by  the  effect  shown  in  Fig.  1;  radicals  with 
F„  >  2  eliminate  hydrogen  during  incorporation  at  a  higher 
rate  than  those  radicals  having  <  2. 

C.  The  burial  effect 

It  has  been  generally  observed  that  for  a  given  gas  mix¬ 
ture,  /h  increases  with  increasing  rate  of  deposition  of  a- 
Si:H.^  ‘’  ■'*  The  deposition  rate  is  usually  increased  by  in- 
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FIG.  3.  Computed  hydrogen  fraction  /„  and  H/Si  ratios  in  the  incident 
flux  as  a  function  of  gas  mixture.  [  H  ]  /  [  Si )  is  the  ratio  of  all  hydrogen  atoms 
to  silicon  atoms  incident  onto  the  surface;  [Si, H„,  c2]/[Si.H„, 

Fh  >  2  ]  is  the  ratio  of  radicals  that  have  Fh  =  H/Si  <2  and  Fh  >  2. 


creasing  discharge  power.  Therefore,  the  increase  in  /„  can 
result  from  either  a  change  in  the  composition  of  the  radicals 
incident  on  the  surface  or  from  what  we  call  the  burial  effect. 
The  burial  effect  is  a  physical,  as  opposed  to  a  chemical, 
effect.  The  results  of  our  model  are  consistent  with  lattice 
interconnection  with  molecular  hydrogen  elimination  being 
a  major  method  of  removing  hydrogen  from  the  surface.  To 
the  extent  that  the  rate  of  interconnection  and  hydrogen 
elimination  is  finite,  the  rate  of  hydrogen  elimination  is  pro¬ 
portional  to  the  residence  time  of  adjacent  Si-H(L)  bonds 
on  the  surface.  The  higher  the  rate  of  deposition,  the  shorter 
the  residence  time  of  Si-H(L)  bonds  on  the  surface  prior  to 
burial  and  the  smaller  the  probability  that  hydrogen  is  eli¬ 
minated  from  the  surface. 

The  burial  effect  is  illustrated  in  Fig.  4,  where  we  have 
plotted  /„  as  a  function  of  deposition  rate.  Figure  4(a) 
shows  results  for  idealized  cases  of  the  flux  of  radicals  being 
composed  of  a  single  species.  Figure  4(b)  shows  simulated 
results  for  and  the  experimental  results  of  Ross  and  Jaklik.^ 
The  simulated  fluxes  were  generated  with  the  model  de¬ 
scribed  in  Ref  17.  In  both  cases  /„  increases  with  increasing 
deposition  rate.  For  these  cases  we  allow  insertion  of 
SiHjfiV)  into  Si-H(Z.)  with  Hj  elimination.  The  lattice  /„ 
for  SiHjfF)  is  therefore  less  than  that  for  SiH,(F)  due  to  H; 
elimination  during  incorporation;  /h  in  isolated  configura¬ 
tions  is  greater  for  SiHjlF)  in  the  absence  of  this  effect.  Over 
the  indicated  range  of  power  deposition,  in  the  simulated 
results  [Fig.  4(b)  ]  the  fractional  composition  of  the  radical 
flux  changes  by  only  a  few  percent.  The  dominant  effect  re¬ 
sponsible  for  increasing  /„  is  the  burial  eff.jt  For  the  re¬ 
sults  in  Fig.  4,  gas  phase  molecules  did  not  saturate  the  sur¬ 
face.  That  effect  will  be  discussed  in  the  next  section. 

D.  Saturation  and  etching  of  lattice  bonds  by  gas  phase 
molecules 

The  reaction  of  radicals  and  molecules  from  the  plasma 
which  collide  with  but  do  not  adsorb  on  the  surface  are  im¬ 
portant  because  they  can  cither  remove  hydrogen  from 
(etch)  the  surface  by  removing  a  hydrogen  atom  from  a  Si- 
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FIG.  4.  Hydrogen  atomic  ratio  in  a-Si:H  films  as  a  function  of  discharge 
power,  demonstrating  the  burial  effect,  (a)  /„  for  individual  radical  spe¬ 
cies,  (b)  this  model's  results  and  experimental  values  (Ref.  7)  for  depo¬ 
sition  of  u-Si:H  from  an  rf  discharge  in  SiH,. 

H(Z.)  bond,  or  deposit  hydrogen  on  (saturate)  the  surface 
by  inserting  a  hydrogen  atom  into  a  Si-D(L)  bond.'^  An 
e,xample  of  these  two  processes  are  the  forward  (etching) 
and  reverse  (saturation)  reactions 

SiH,(F)  -|-SSSH(I);^iHj(F)  -t-SSSD(L).  (13) 

Since  a  Si-H(Z. )  bond  on  the  surface  is  exchanged  for  a  Si- 
H  ( F)  bond  in  the  gas  phase  radical,  the  exothermicity  of  the 
class  of  reactions  with  radicals  is  determined  by  the  differ¬ 
ence  in  bond  energies  for  SiH„H(F)  and  SSS-H(£).  Rob¬ 
ertson  and  Gallagher  assigned  a  bond  energy  of  =:  80  kcal/ 
mol  to  SSS-H ( L) . The  SiHj-H (/^  bond  energy  is  approx¬ 
imately  75  kcal/mol.  To  the  level  of  accuracy  of  the  approxi¬ 
mations,  this  class  of  reactions  involving  radicals  can  be  con¬ 
sidered  thermoneutral.  The  equilibrium  ratio  of  saturated 
bonds  to  dangling  bonds,  [SSSH(I.)]/[SSSD(I,)1,  result¬ 
ing  from  these  processes  involving  radicals  is  therefore 
(SiH,  (/■)  )/[SiH(, _  (fO ].  For  typical  plasma  condi¬ 
tions*’  this  ratio  is  >  10,  so  that  the  equilibrium  density  of 
saturated  to  dangling  bonds  resulting  from  these  radical  re¬ 
actions  is  [SSSH(I.)  1/[SSSD(£,)]  >  10.  Therefore,  surface 


bonds  are  more  likely  to  be  saturated  than  etched  by  radicals 
via  this  process.  Even  though  the  steady-state  density  of  dan¬ 
gling  bonds  made  available  by  this  process  is  likely  to  be 
small,  the  process  may  be  an  important  source  of  dangling 
bonds  for  the  case  where  adsorbed  radicals  as  SiHj  ( A/)  can¬ 
not  insert  into  saturated  surface  sites  [e.g.,  SSSH(£.)  ].'* 

Analogous  etching  and  saturation  reactions  can  also  oc¬ 
cur  by  collisions  of  saturated  molecules  with  the  lattice.  Ex¬ 
amples  of  these  reactions  are 

SiH4(F)  -f  SSSH(L) 

-SiHjCF)  -I-  H2((7)  +  SSSD(£)  (etching),  ( 14) 
-f  SSSD(£) 

— SiH3(F) -h SSSH(£)  (saturation).  (IS) 

The  etching  reaction  in  Eq.  ( 14)  is  endothermic,  since  two 
Si-H  bonds  are  exchanged  for  the  single  H-H  bond.  The 
saturation  reaction  in  Eq.  (15),  though,  is  approximately 
thermoneutral  in  the  same  manner  as  for  the  analogous  reac¬ 
tion  involving  radicals  in  Eq.  (13).  Since  SiH4  is  the  feed 
stock,  the  large  ratio  of  densities  of  [SiH4(/0]/[SiH3(/^] 
insures  that  the  equilibrium  density  of  saturated  bonds  com¬ 
pared  to  dangling  bonds  resulting  from  this  process  is  large; 
the  relative  etching  probability  is  small,  although  the  contri¬ 
bution  may  be  an  important  source  of  dangling  bonds  (see 
above) .  From  this  discussion,  one  can  conclude  that  nonad¬ 
sorbing  gas  phase  radicals  and  molecules  colliding  with  the 
growing  film  are  more  likely  to  contribute  hydrogen  to  the 
surface  than  remove  it. 

The  saturation  of  dangling  bonds  on  the  surface  of  the 
lattice  by  gas  phase,  nonadsorbing  molecules  clearly  has  the 
potential  of  being  an  important  source  of  hydrogen  to  the 
film.  We  investigated  this  process  by  assigning  a  probability 
of  saturating  Si-D(£)  bonds  with  hydrogen  resulting  from 
collisions  of  SiH4(F)  and  HjlF),  the  saturated  molecules 
having  the  largest  flux  striking  the  surface.  The  rate  of  satu¬ 
ration  of  dangling  bond  Xj(L)  was  specified  to  be  p, 
((SiH4(F)l  -f  [Xj(L)]/pa,  where  p,  is  the  satura¬ 
tion  probability.  The  effect  on  /h  of  varying  p,  is  shown  in 
Fig.  5.  The  case  shown  is  the  same  as  in  Fig.  4(b)  for  Ar/ 


SATURATION  PROBABILITY 

FIG.  S.  Computed  hydrogen  atomic  fraction  as  a  function  of  the  saturation 
probability  of  dangling  bonds  [e.g.,  SSSD(  I. ) )  by  nonadsorbing  molecules 
form  the  plasma  ( Hj,  SiH4 ) .  The  plasma  conditions  are  the  same  as  for  Fig. 
3.  The  partial  composition  of  the  radical  and  molecular  flux  to  the  surface  is 
shown  at  right. 
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SiH4  =  0.75/0.25.  The  composition  of  the  flux  is  also  shown 
in  the  figure.  The  fraction  of  incorporated  hydrogen,  /h, 
increases  significantly  when  the  probability  for  saturation  is 
>  io-\ 

The  fraction  of  hydrogen  in  isolated  configurations  also 
increases  with  increasing  lattice  saturation  probability,  al¬ 
though  to  a  lesser  degree  than  /h  of  the  lattice.  The  increase 
of  /h  in  isolated  configurations  is  due  to  the  adsorbed  mo¬ 
bile  radicals  [X{M)]  having  a  smaller  probability  of  en¬ 
countering  a  lattice  site  having  a  dangling  bond  when  the 
saturation  probability  is  high.  Since  incorporation  of  X{M) 
into  a  saturated  bond  has  a  smaller  rate  constant,  the  resi¬ 
dence  time  and  density  of  X(M)  on  the  surface  increases, 
thereby  increasing  the  probability  of  being  buried  as  an  iso¬ 
lated  configuration. 

D.  Substrate  temperature 

Computed  hydrogen  fraction  /h  as  a  function  of  sub¬ 
strate  temperature  T,  for  the  experimental  conditions  of 
Ross  and  Jaklik^  appear  in  Fig.  6.  decreases  with  increas¬ 
ing  substrate  temperature  due  to  the  activation  energy  for 
interconnection  and  hydrogen  elimination  of  Si-H(Z.) 
bonds.  The  activation  energy  derived  here,  1.5  kcal/mol,  is 
approximately  equal  to  that  empirically  derived  by  measur¬ 
ing  /h  as  a  function  of  T„  1.35  kcal/mol. '  ‘  The  decrease  of 
hydrogen  in  isolated  or  polymeric  configurations  results 
from  a  decrease  in  the  residence  time  of  SiH„  (A )  as  immo¬ 
bile  adsorbates;  a  consequence  of  the  activation  energy  for 
desorption. 

The  relative  densities  of  hydrides  and  dihydrides  is  also 
a  function  of  substrate  temperature,  as  shown  in  Fig.  7.  In 
this  figure  the  experimental  results  of  Luveosky,  Nemanich, 
and  Knights'^  and  results  from  the  model  for  [SSSH(J3) )/ 
(SSHH(5)1  are  plotted  as  a  function  of  T,  and  rf  power 
(deposition  rate).”  The  relative  hydride  density  increases 
with  increasing  T,  and  is  nearly  constant  for  discharge  pow¬ 
er  greater  than  5  W  (deposition  rate  of  150  A/min).  To 
obtain  these  results,  it  was  necessary  to  invoke  a  higher  acti¬ 
vation  energy  ( 10  kcal/mol)  for  interconnection  and  hydro¬ 
gen  elimination  between  di-  and  trihydride  [SSHH(L), 
SHHH  (£)]  configurations  than  between  hydride  [SSSH 
(£))  configurations  (1.5  kcal/mol).  If  the  method  of  hy- 


FIG.  6.  Hydrogen  atomic  fraction  as  a  function  of  substrate  temperature  as 
computed  with  the  surface  deposition  model  and  the  experimental  results  of 
Ross  and  Jaklik  (Ref.  7). 


SUBSTRATE  TEMPERATURE  (K) 


FIG.  7.  The  ratio  of  hydride  (SSSH(^)]  and  dihydride  (SSHH(£))  con¬ 
figurations  in  a-Si;H  films  as  calculated  with  the  surface  deposition  model 
compared  to  the  experimental  results  of  Lucovsky  ei  at.  ( Ref.  12).  The  two 
cases  are  for  varying  the  rf  power  (deposition  rate)  and  substrate  tempera¬ 
ture  r, .  The  relative  hydride  density  increases  with  increasing  T,.To  obtain 
agreement  with  experiment  for  the  temperature  dependence,  two  activation 
energies  were  invoked;  for  interconnection  between  hydrides,  and  for  inter¬ 
connection  between  di-  and  trihydrides. 


drogen  incorporation  is  a  weak  function  of  temperature, 
then  a  uniform  activation  energy  for  interconnection 
between  all  hydride  configurations  nets  only  a  weak  tem¬ 
perature  dependence  for  their  relative  densities.  It  is  possible 
that  the  difference  in  activation  energies  between  hydrides 
and  dihydrides  lies  not  in  the  interconnection  process  but  in 
the  incorporation  step.  We  assumed,  though,  that  the  latter 
process  had  no  activation  energy.  The  suggestion  that  acti¬ 
vation  energies  for  mono-  and  polyhydrides  are  different  is 
consistent  with  the  observation  that  above  a  substrate  tem¬ 
perature  of  575  K,  hydrogen  evolves  from  HOMOCVD 
films  with  a  higher  activation  energy  (30  kcal/mol)  than 
below  that  temperature.  ‘ ' 

E.  Reactive  sputter  deposition 

Sputter  deposition  of  a-Si.-H  differs  qualitatively  from 
plasma  deposition  of  a-Si:H  in  that  the  operator,  by  adjust¬ 
ing  target  voltage  and  hydrogen  flowrate,  has  finer  control  of 
the  composition  of  radicals  incident  onto  the  substrate.  Un¬ 
like  PECVD,  where  a  large  fraction  of  the  hydrogen  incident 
onto  the  surface  is  contained  in  silane  radicals,  in  reactive 
sputter  deposition  the  hydrogen  and  silicon  fluxes  to  the  sur¬ 
face  are  dominantly  atomic.  For  sputter  deposition,  the  scal¬ 
ing  parameter  to  characterize  the  incident  flux  is  ,  the 
ratio  of  atomic  H  ( F)  to  atomic  SHF).  The  computed  atomic 
hydrogen  content  of  sputter  deposited  a-Si:H  is  plotted  in 
Fig.  8  as  a  function  of  F^  for  three  different  substrate  tem¬ 
peratures.  The  silicon  flux  is  constant  at  2.4  x  10'*  cm  “  ’  s" ' 
corresponding  to  a  deposition  rate  of  300  A/min.  As  expect¬ 
ed,  the  fraction  of  hydrogen  in  both  the  lattice  and  in  isolated 
configurations  increases  with  increasing  F  h  ■  At  higher  sub¬ 
strate  temperatures,  though,  /h  for  the  lattice  reaches  a 
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FIG.  8.  Computed  hydrogen  atomic  fraction  and  hydride/dihydride  den¬ 
sity  for  <i-Si;H  Rims  grown  by  reactive  sputter  deposition  as  a  function  of  the 
atomic  hydrogen  flux.  The  Si  flux  is  constant  at  2.4x  10'^  cm  "' s~',  corre¬ 
sponding  to  a  deposition  rate  of  300  A/min. 


nearly  constant  value  at  whereas  for  isolated 

configurations  continues  to  increase.  The  saturation  of  /h  is 
due  in  part,  to  an  equilibration  between  adsorption  and  de¬ 
sorption  of  resulting  from  collisions  between  H{M) 
radicals;  a  process  proportional  to  [H(/^  and  which  in¬ 
creases  with  increasing  T,. 

The  relative  densities  of  hydride  and  dihydride  configu¬ 
rations  [i.e.,  SSSH(.5)and  SSHH(.fl)]  in  the  sputtered  a- 
Si:H  films  are  also  plotted  in  Fig.  8.  At  low  values  of  FJ, ,  the 
films  are  composed  of  dominantly  hydride  [SSSH(F)  j  con¬ 
figurations,  and  with  increasing  larger  fraction  of  hy¬ 

drogen  is  contained  in  the  film  as  dihydride  [SSHH(F)]. 
The  cause  for  the  increase  in  dihyride  density  with  increas¬ 
ing  Fli  (at  constant  deposition  rate)  is  analogous  to  the 
burial  effect.  Given  a  finite  residence  time  of  dangling  bonds 
on  the  surface  prior  to  burial,  a  higher  flux  of  hydrogen  pas¬ 
sivates  a  larger  fraction  of  those  bonds,  thereby  leading  to 
the  higher  dihydride  density.  If  so,  then  at  higher  deposition 
rate  (shorter  residence  time  of  dangling  bonds  on  the  sur¬ 
face)  or  higher  substrate  temperature  (larger  rate  of  inter¬ 
connection),  the  ratio  [SSSH(F)1/[SSHH(F)]  should  de¬ 
crease.  This  trend  is  confirmed  by  comparing  the  computed 
ratio  [SSSH(F)  ]/[SSHH(F)  ]  to  the  experimental  results 
of  Jeffrey,  Shanks,  and  Danielson*  in  Fig.  9.  For  these  re¬ 
sults,  the  ratio  H(F)/Si(F)  was  obtained  by  normalizing  to 
the  experiment  at  a  deposition  of  2:330  A/min  (H(F)/ 
Si(F)  =  3],  Results  as  a  function  of  deposition  rate  were 
then  obtained  both  while  holding  H(F)  constant  and  assum¬ 
ing  H(F)  was  proportional  to  the  change  in  discharge  pow¬ 
er,  which  is  directly  proportional  to  deposition  rate. 


F.  Deposition  without  insertion  of  into 

saturated  surface  sites 

The  details  of  the  method  of  incorporation  of  SiH,(iV) 
into  the  lattice  are  particularly  important,  because  SiHj(F) 
has  the  highest  flux  of  silicon  containing  radicals  incident 
onto  the  surface.  In  Sec.  V  A  we  discussed  the  implications 
of  having  SiHj(A/)  incorporate  into  the  lattice  by  inserting 
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FIG.  9.  Computed  and  experimental  results  for  the  ratio  of  hydride 
(SSSH(S)  ]  to  dihydride  (SSHH(5)  ]  density  in  a-Si:H  films  grown  by  re¬ 
active  sputter  deposition.  The  experimental  results  are  from  Jeffrey  et  al. 
(Ref.  8).  The  ratio  of  atomic  H/Si  of  the  radical  flux  for  the  model  was 
obtained  by  normalizing  to  the  experimental  results  at  a  desposition  rate  of 
330  A/min.  Two  case  are  shown,  scaling  H/Si  with  deposition  rate  (dis¬ 
charge  power)  and  holding  the  hydrogen  flux  constant. 


into  Si-H(£.)  bonds  and  eliminating  hydrogen,  as  in  reac¬ 
tions  (11)  and  ( 12).  In  this  section  we  will  discuss  the  impli¬ 
cations  of  SiH,(ti/)  incorporation  in  the  absence  of  this 
mechanism.  For  these  conditions,  SiH3(A/)  can  only  be  in¬ 
corporated  into  the  lattice  at  a  surface  site  having  a  dangling 
bond,  as  in  the  reaction 

SiHj(W)  -(-SSSD(£,)-SSSS(f.)  -bSHHH(L).  (16) 

The  dangling  bonds  at  surface  sites  are  made  available  by 
etching  [reactions  (13)  and  (14)]  or  sputtering  by  bom¬ 
bardment  of  the  surface  with  energetic  ions  ( see  Table  II ) . 

When  SiH}(M)  cannot  insert  into  Si-H(L )  bonds,  the 
rate  at  which  SiH,(Af)  incorporates  into  the  lattice  is 
smaller  since  the  availability  of  the  necessary  surface  sites 
[eg.,  SSSD(Z, )  ]  is  also  smaller.  For  a  constant  SiH)(F)  flux, 
the  surface  density  of  SiHjfM)  therefore  increases.  Surface 
clustering  and  desorption  reactions  therefore  also  become 
more  likely  (see  Table  II.)  The  two  desorption  reactions 
which  are  important  are 

SiH,(Af)  -I-  SiH3(M)--SijH»(G).  (17a) 

SiHj(A/)  -I-  H(A/)  -SiH«(C).  ( 17b) 

As  a  result  of  these  reactions,  the  effective  sticking  coeffi¬ 
cient  for  SiH}(Af)  is  less  than  unity. 

We  investigated  the  characteristics  of  a-Si:H  films  when 
SiHylAf)  incorporation  occurs  only  via  a  dangling  bond. 
These  results  are  shown  in  Fig.  10  where  the  film  deposition 
rate  and  effective  sticking  coefficient  for  SiH3(M)  are  plot¬ 
ted  as  a  function  of  the  rate  of  generation  of  dangling  bonds. 
The  flux  of  radicals  incident  onto  the  surface  have  the  ratio 
SiH3/SiH2/H  =  100/5/10  with  [SiH3(F)]  =  10'*  cm"* 
s~ These  values  are  typical  for  a  low-powered  plasma  in  a 
SiH^/Hj  =  1/1  gas  mixture  at  0.25  Torr.  The  independent 
variable  is  the  number  of  dangling  bonds  generated  on  the 
surface  by  sputtering  or  etching  as  a  fraction  of  the  SiH3  (F) 
flux  incident  onto  the  surface.  As  the  rate  of  generation  of 
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FIG.  10.  Deposition  rate  and  effective  slicking  coefficient  for  SiH,(<V)  as  a 
function  of  the  rate  of  generation  of  dangling  bonds.  For  these  results, 
SiH,(.lf)  does  not  insert  inloSi-H(£).  The  discharge  conditions  are  SiH,/ 
(F)/SiH:(n/H(F)  =  iOO/5/IOwithSiHj(F)  =  IO'*cm-^s-‘.Therale 
of  generation  of  dangling  bonds  is  expressed  as  a  fraction  of  SiHjiF). 

dangling  bonds  increases  and  the  number  of  sites  available 
for  SiHjliV/)  insertion  increases,  both  the  deposition  rate 
and  effective  sticking  coefficient  of  SiHjC/O  also  increase. 
Over  the  range  investigated,  though,  there  is  little  change  in 
hydrogen  fraction  ( 9.5%  </„  <10.0%).  This  last  observa¬ 
tion  is  explained  by  the  fact  that  if  a  dangling  bond  is  re¬ 
quired  for  incorporation  of  SiHjCA/),  then  the  two  reactions 
responsible  for  deposition  for  these  conditions, 

SiH,(.V/)  -|-SSSH(I)-SSSS(L)  +SHHH(Z,)  (I8a) 

SiH3(.W)  -f  SSSD(Z,)-SSSS(L)  -f-SHHH(L),  (18b) 

are  functionally  equivalent  with  respect  to  "hydrogen  ac¬ 
counting."  The  SiHjC.V)  which  does  not  stick  is  desorbed 
primarily  as  disilane,  as  the  desorption  products  have  the 
ratio  [Si;Hft(G) ]/[SiH4(G) ]  =6. 

VI.  CONCLUDING  REMARKS 

A  phenomenological  model  for  surface  kinetics  during 
plasma  and  sputter  CVD  of  amorphous  hydrogenerated  sili¬ 
con  has  been  described.  The  model  consists  of  an  accounting, 
in  rate  equation  form,  of  adsorption,  incorporation,  inter¬ 
connection,  and  other  pertinent  processes.  The  model  has 
successfully  compared  with  experimental  results  for  hydro¬ 
gen  fraction  and  distribution  of  bond  types.  We  found  that 
molecular  hydrogen  elimination  during  incorporation  dis¬ 
tinguishes  between  films  grown  from  plasmas  whose  pri¬ 
mary  radicals  SiH„  have  «<  or  /» >  2.  We  also  found  that  the 
observed  increase  in  hydrogen  fraction  with  increasing  de¬ 
position  rate  is  largely  a  physical,  as  opposed  to  chemical, 
process.  At  higher  deposition  rates,  the  shorter  residence 
time  of  Si-H(£)  bonds  on  the  surface  results  in  a  smaller 
probability  for  interconnection  and  hydrogen  elimination. 
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The  cause  for  the  increase  in  the  fraction  of  hydrogen  con¬ 
tained  in  hydride,  as  opposed  to  dihydride,  configurations 
with  increasing  substrate  temperature  may  be  a  result  of 
there  being  a  higher  activation  energy  of  interconnection 
between  dihydride  configurations. 
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A  model  for  the  plasma  enhanced  chemical  vapor  deposition  of  amorphous  hydrogenated 
silicon  (a-Si:H)  in  rf  and  dc  discharges  is  presented.  The  model  deals  primarily  with  the 
plasma  chemistry  of  discharges  sustained  in  gas  mixtures  containing  silane  (SiH4) .  The  plasma 
chemistry  model  uses  as  input  the  electron  impact  rate  coefficients  generated  in  a  separate 
simulation  for  the  electron  kinetics  and  therefore  makes  no  a  priori  assumptions  as  to  the 
manner  of  power  deposition.  Radical  densities  and  contributions  to  film  growth  are  discussed 
as  a  function  of  gas  mixture,  electrode  separation,  and  locale  of  power  deposition,  and 
comparisons  are  made  to  experiment.  A  compendium  of  reactions  and  rate  constants  for  silane 
neutral  and  ion  chemistry  is  also  presented. 


I.  INTRODUCTION 

Plasma  enhanced  chemical  vapor  deposition  ( PECVD) 
is  a  process  whereby  thin  films  of  amorphous  hydrogenated 
silicon  (<2-Si:H )  are  made  at  substrate  temperatures  of  <S(X) 
K,  significantly  below  than  those  required  for  conventional 
chemical  vapor  deposition  (CVD).'"*  PECVD  is  usually 
performed  in  parallel  plate,  capacitively  coupled  radio  fre¬ 
quency  ( rf)  discharges  sustained  in  gas  mixtures  containing 
silane  (SiH4)  or  disilane  (SiiH*)  (see  Fig.  1).  Deposition 
rates  of  10's-l(X)’s  A/min  can  be  obtained  from  discharges 
in  pure  silane  of  few  hundred  mTorr  with  power  deposition 
of  tens  to  hundreds  of  mW  cm  ”  \  Films  so  produced  have 
atomic  hydrogen  fractions  of  5%-20%,  and  are  used  in  elec¬ 
tronic  devices  (thickness  =:0.5  pm),  photovoltaic  cells 
(0.2-0.5  /rm),  and  as  photoreceptors  (10-15  pm).  Solar 
cells  of  single  junction  a-Si:H  fabricated  by  PECVD  have 
demonstrated  efficiencies  of  1 1.7%  in  small  areas  ( 1  cm~), 
and  7,9%  in  large  area  modules  (  >  1000  cm’).* 

The  processes  leading  to  the  PECVD  of  a-Si:H  begin 
with  electron  impact  dissociation  of  silane,  creating  neutral 
and  charged  radicals  consisting  principally  of  silyl,  SiHj. 
The  radicals  undergo  neutral  chemical  and  ion-molecule  re¬ 
actions  while  they  diffuse  or  drift  to  a  substrate,  usually  one 
of  the  parallel  plate  electrodes.  At  the  substrate,  the  radicals 
are  adsorbed  and  incorporated  into  the  lattice  of  the  film.  In 
principal,  PECVD  of  a-Si;H  is  divisible  into  three  sequential 
subprocesses:  discharge  electron  kinetics,  plasma  chemistry, 
and  surface  deposition  kinetics.  Electron  impact  collisions 
from  the  first  generate  radicals  which  subsequently  react  in 
the  second  and  which  are  finally  incorporated  into  the  film  in 
the  third.  In  reality,  the  three  subprocesses  are  connected  in 
a  more  global  sense  than  a  simple  sequence  of  events.  For 
example,  depletion  of  the  feedstock  gases  by  electron  impact 
dissociation  cause  changes  in  the  electron  distribution  func¬ 
tion,  which  results  in  a  change  in  the  rate  of  production  of 
radicals.  Surface  reactions  which  cause  saturated  molecules 
to  desorb  into  the  gas  may,  in  turn,  modify  the  chain  of 
chemical  reactions.'’ 

The  electron  kinetics,  plasma  chemistry,  and  surface  ki¬ 
netics  for  the  PECVD  of  a-Si:H  have  been  the  topic  of  many 


previous  investigations.  In  particular,  discharge  and  chemi¬ 
cal  kinetics  models  have  been  developed  with  the  goal  of 
obtaining  a  theoretical  understanding  of  the  process.  Tur¬ 
ban,  Catherine,  and  Grolleau^  developed  the  first  such  mod¬ 
el  in  1979.  In  this  model,  all  processes  were  assumed  to  be 
spatially  uniform  and  a  semi-empirical  rate  constant  for 
electron  impact  dissociation  of  SiH4  was  used.  At  the  time  of 
the  study,  the  cross  section  for  electron  impact  dissociation 
of  silane  had  not  yet  been  measured.  Ion  chemistry  and  neu¬ 
tral  chain  reactions  in  silane  discharges  were  analyzed  by 
Haller*’*’  using  rate  equation  models  in  1980  and  1983.  He 
concluded  that  reactions  of  radicals  following  electron  im¬ 
pact  dissociation  of  silane  are  equally  if  not  more  important 
in  determining  the  characteristics  of  the  deposited  a-Si:H  as 
the  direct  products  of  the  dissociation.  Chen,'”  also  in  1983, 
analyzed  PECVD  from  a  mass  transport  viewpoint,  examin¬ 
ing  various  reactor  geometries  while  not  emphasizing  the 
details  of  the  plasma  chemistry.  Longeway,  Estes,  and 
Weakliem"  analyzed  the  decomposition  of  silane  in  a  dc 
glow  discharge  in  1984  using  a  steady-state  rate  equation 
model,  and  concluded  that  the  primary  dissociation  product 
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FIG.  I .  Schematic  of  a  plasma  deposition  reactor.  A  radio  frequency  power 
supply  is  capacitively  coupled  through  an  impedance  matching  network 
(IMN)  to  a  set  of  parallel  plate  electrodes  separated  by  a  few  to  10  cm. 
Reactant  gas  mixtures  containing  silane  (SiH,)  at  pressures  of  lO's  to  lOO's 
of  mTorr  flow  between  the  electrodes  with  residence  times  of  lO's  of  ms  to  a 
second.  Silane  is  dissociated  by  electron  impact  resulting  in  deposition  of  an 
a-Si:H  film.  The  products  and  undissociated  silane  are  pumped  out  of  the 
reactor. 
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of  silane  is  SiHj.  In  1984,  Tachibana'*  presented  a  quasi- 
steady-state  plasma  chemistry  model  for  rf  discharges  in  si¬ 
lane,  including  higher  silane  species  (through  SijH,;  and 
SijH^'  ).  Although  the  neutral  chemistry  and  ion-molecule 
reactions  in  his  model  are  treated  in  a  self-consistent  man¬ 
ner,  the  electron  density  was  an  input  parameter  whose  value 
was  based  on  experimental  electrical  probe  measurements. 
Electron  impact  rate  constants  were  obtained  by  convolving 
a  semi-empirical  electron  distribution  function  with  cross 
sections.  Tachibana  concluded  that  SiHj  is  the  most  abun¬ 
dant  gas  phase  radical  and  that  SiHj  is  likely  to  directly 
insert  into  the  growing  film.  The  following  year,  Chang  el 
al.'^  modeled  positive  column  discharges  of  silane  diluted  by 
argon,  helium,  and  neon.  They  assumed  a  Maxwellian  elec¬ 
tron  distribution  and,  with  the  quasi-neutrality  condition, 
solved  for  the  electron  temperature  as  a  function  of  electron 
density.  Values  of  1-2  eV  were  obtained.  They  found  that 
silane  ions  have  a  small  density  compared  to  neutral  radi¬ 
cals,  and  excitation  transfer  to  silane  from  noble  gas  metasta¬ 
bles  is  most  important  for  He-SiH4  mixtures.  Ion  chemistry 
in  the  sheaths  of  dc  discharges  was  studied  in  198S  by 
Chatham  and  Gallagher  using  a  one-dimensional  transport 
model.  They  parametrically  calculated  the  discharge  oper¬ 
ating  characteristics  which  yield  given  fractions  of  lower  and 
higher  silane  ions. 

The  electron  kinetics  of  silane  discharges  have  also  been 
separately  theoretically  studied  by  Garscadden,  Duke,  and 
Baily,*'  Hayashi,‘*Ohmori,  Shimozuma,  and  Tagashima,’’ 
and  Kushner. These  studies  resulted  in  derived  sets  of  elec¬ 
tron  impact  cross  sections,  and  the  ability  to  generate  elec¬ 
tron  impact  rate  constants  as  a  function  of  position  and  time 
in  rf  discharges. 

A  reoccurring  conclusion  from  previous  modeling  stud¬ 
ies  is  that  there  are  interdependencies  between  the  subpro¬ 
cesses  (electron  kinetics,  plasma  chemistry,  and  surface  ki¬ 
netics)  which  put  a  premium  on  performing  as 
self-consistent  an  analysis  as  possible.  Recognizing  this  fact, 
previous  investigators  augmented  their  models  by  using  ex¬ 
perimental  data  for  variables  such  as  electron  density  or  dis¬ 
sociation  rates.  A  totally  self-consistent  analysis  has  thus  far 
been  lacking. 

In  this  paper  a  model  for  the  discharge  kinetics  and  plas¬ 
ma  chemistry  of  rf  and  dc  discharges  sustained  in  silane  gas 
mixtures  for  the  PECVD  of  a-Si;H  is  presented.  This  model 
differs  from  previous  treatments  in  that  no  assumptions  are 
made  concerning  parameters  such  as  the  electron  distribu¬ 
tion  function,  dissociation  rate,  reactor  geometry,  or  reac¬ 
tion  chain.  The  goal  of  the  model  is  to  simulate  the  PECVD 
of  a-Si:H  from  discharges  sustained  in  silane  while  making  a 
minimum  of  a  priori  assumptions.  The  simulation  consists  of 
a  one-and-one-half-dimensional,  time-dependent  plasma 
chemistry  model  whose  electron  impact  rate  constants  are 
obtained  from  a  Monte  Carlo  simulation  for  the  electron 
distribution  function  in  parallel  plate  rf  discharges.'"''*  The 
plasma  chemistry  model  includes  a  sufficient  number  and 
variety  of  species  and  reactions  that  gas  mixtures  of  SiH^  or 
3i,H„  diluted  by  Ar  or  H:  can  be  simulated.  Although  some 
of  the  reaction  rates  are  estimates,  the  model  is  so  construct¬ 
ed  that  as  new  measurements  of  these  rates  become  available 


they  can  be  easily  incorporated  into  the  simulation. 

The  model  uses  as  input  gas  mixture,  reactor  geometry, 
and  power  deposition:  and  produces  as  output  the  density  of 
radicals  and  ions  as  a  function  of  time  and  position  between 
the  electrodes,  the  flux  of  radicals  and  ions  incident  onto  the 
substrates  or  electrodes,  and  the  deposition  rate  of  a-Si:H. 
The  flux  of  radicals  incident  onto  the  electrodes  can  also  be 
used  as  input  to  a  separate  model  for  surface  deposition  ki¬ 
netics.  This  second  model  simulates  characteristics  of  the  o- 
Si:H  film  such  as  hydrogen  fraction  and  dihydride  density. 
The  three  submodels  for  electron  kinetics,  plasma  chemis¬ 
try,  and  surface  deposition  kinetics,  represent  the  first  at¬ 
tempt  at  a  totally  self-consistent  model  for  the  PECVD  of  a- 
Si:H.-' 

In  Sec.  II,  the  plasma  chemistry  model  is  described  with 
an  accompanying  discussion  on  the  source  and  selection  of 
rate  constants,  cross  sections,  and  transport  coefficients.  The 
intent  of  the  discussion  and  listing  of  reactions  is  to  summa¬ 
rize  the  current  availability  of  rate  constants  and  cross  sec¬ 
tions  for  silane  plasma  chemistry.  In  Sec.  Ill,  a  sampling  of 
results  from  the  model  are  discussed  for  various  aspects  of 
the  PECVD  of  fl-Si:H.  Time  and  spatially  dependent  densi¬ 
ties  of  radicals  and  saturated  molecules  are  presented  for  rf 
discharges,  followed  by  a  discussion  of  the  rate  of  survival  of 
radicals.  Discharge  properties  as  a  function  of  electrode  sep¬ 
aration  and  radical  contributions  to  film  growth  are  dis¬ 
cussed,  followed  by  a  discussion  on  the  formation  of  higher 
silane  ions. 


II.  DESCRIPTION  OF  THE  MODEL 


The  plasma  chemistry  model  consists  of  a  set  of  partial 
diflerential  equations  ( PDEs)  describing  the  time  evolution 
of  the  density  of  individual  radical  and  saturated  species  as 
the  input  gases  flow  downstream  in  an  rf  glow  discharge. 
The  model  principally  simulates  a  capacitively  coupled  par¬ 
allel  plate  (diode)  rf  glow  discharge,  and  by  modifying  the 
boundary  conditions,  dc  discharges  are  simulated.  We  as¬ 
sumed  plug  gas  flow  in  the  downstream  direction,  thus  inte¬ 
gration  in  time  is  equivalent  to  integrating  in  space  along  the 
flow  direction  by  specifying  the  flowrate  and  dimensions  of 
the  reactor. 

The  PDE  for  each  species  has  the  general  form  of 


dN,(x,t) 

Jt 


-N, 


+  '7  D,VN, -V  p,ENi,  (1) 

where  all  quantities  are  a  function  of  position  and  time.  In 
Eq.d),*,  is  the  electron  impact  rate  coefficient  for  produc¬ 
tion  of  species  i  by  electron  collisions  with  species  j,  k,j,  is  the 
analogous  coefficient  for  reactions  between  heavy  particles  j 
and  /.  D,  is  the  diffusion  coefficient  for  species  i,  p,  is  the 
mobility  of  charged  species  /,  and  E  is  the  local  electric  field. 
The  first  and  third  terms  in  Eq.  ( 1 )  are  for  electron  impact 
processes  (excitation,  dissociation,  ionization,  attachment) 
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resulting  in  generation  and  depletion  of  species  ;  the  sec¬ 
ond  and  fourth  terms  are  for  generation  and  depletion  by 
neutral  and  ion  heavy  particle  reactions;  the  fifth  term  is  for 
transport  by  diffusion  for  neutrals  and  ions;  and  the  sixth 
term  is  for  transport  of  ions  by  drift  in  the  electric  field. 
Using  a  finite  difference  formulation  for  the  spatial  deriva¬ 
tives,  the  set  of  PDEs  was  reduced  to  a  set  of  ordinary  differ¬ 
ential  equations.  Initial  conditions  were  assigned,  and  the 
system  of  equations  was  integrated  in  time  using  a  third- 
order  Runga-Kutta  technique. 

The  diffusion  approximation  for  gas  transport  is  valid 
only  for  conditions  where  the  gas  flow  is  laminar,  and  the 
molecular  mean  free  path  is  small  compared  to  other  spatial 
scales  of  interest,  such  as  the  electrode  separation.  For  a  gas 
density  of  1.5  X  10”  cm~^  (F  =  100  mTorr,  F,  =  650  K), 
the  molecular  mean  free  path  is  s;0.2  cm.  The  diffusion  ap¬ 
proximation  is  therefore  valid  for  electrode  separation  d, 
obeying  £f>0.2  (0.1/F(Torr)  ]  [r^(K)/650]  cm. 

A.  Species 

The  species  included  in  the  model  are  listed  in  Table  I. 
The  species  and  reactions  used  here  are  sufficient  to  simulate 


TABLE  I.  Species  in  the  plasma  chemistry  model. 


Argon  species 

Hydrogen  species 

Ar 

H, 

Ar*(3i) 

H? 

Ar*»(4p) 

H,* 

Ar* 

H,‘ 

ArH* 

H 

H« 

H* 

H- 

Silane  species 

SiH, 

SiH,(vl.3) 

SiH,(t>2.4) 

SiH, 

SiH; 

SiH 

Si 

SiH,* 

SiH,* 

SiH* 

Si* 

SiH,- 

SiH,- 

Disilane  species 

Si.H. 

Si,H,(i;l.3) 

Si2H,(w2.4) 

Si,H, 

SijH, 

Si,H, 

Si.H,* 

Si,H,* 

Si.H,* 

Si,H* 

Si,H,* 

Si,H,* 

Si.H* 

Higher  silane  species 

Si.H, 

Si,H, 

Si,H, 

Si.H,o 

Si,H,, 

Si,H„ 

DUST 

Si.H,* 

Si.H,* 

Si,H,; 

discharges  sustained  in  gas  mixtures  initially  consisting  of 
SiH4/Si2H6/H2/Ar.  The  species  include  lower  silane  radi¬ 
cals  and  ions,  higher  silane  radicals,  saturated  molecules, 
and  ions,  and  vibrationally  excited  silane  and  disilane.  Vi- 
brationally  excited  silane  and  disilane  are  explicitly  included 
in  the  model  because  a  significant  fraction  of  the  power  can 
be  coupled  into  the  discharge  by  vibrational  excitation.  The 
fraction  of  power  deposited  in  this  fashion  can  be  >  0.2  re¬ 
sulting  in  a  fractional  density  of  vibrationally  excited  silane 
of  <0.05. 

Mass  spectroscopy,^^  laser  scattering  measurements,^^ 
and  qualitative  observations  of  particulate  (or  dust)  in  si¬ 
lane  glow  discharges  are  evidence  for  the  existence  of  higher 
silane  radicals,  molecules,  and  ions.  These  species  range 
from  higher  silane  molecules  containing  <  10  silicon  atoms 
to  clusters  exceeding  many  micrometers  in  size,  thereby  con¬ 
taining  >  lO'”  silicon  stoms.  In  order  to  investigate  the  clus¬ 
tering  process,  we  included  higher  radical  and  saturated  si¬ 
lane  species,  Si,H2,.Hi  Si.Hj,,.,.  d,  respectively,  for 
A  <5.  We  also  included  a  category  labeled  DUST  represent¬ 
ing  molecules  or  clusters  having  greater  than  five  silicon 
atoms.  The  polymerization  chain  linking  the  higher  silanes  is 
based  on  the  insertion  of  lower  silane  radicaU  into  higher 
saturated  molecules,  and  the  saturation  of  higher  radicals  by 
hydrogen.  For  example, 

SiH-l-SuH,o-Si,H,„ 

SiH2 -)- Si4H|o-»Si5H,2,  . 

SiH,-|-Si4H,o-H,H-Si,H,„  ’ 

Si,H„-l-H-Si,H,j. 

Reactions  in  the  polymerization  chain  which  result  in  mole¬ 
cules  having  more  than  five  silicon  atoms  contribute  to  the 
density  of  DUST.  The  polymerization  chain  is  propagated 
with  a  uniform  rate  constant,  which  is  k,  =  10“ ' '  cm’  s~ ' 
unless  otherwise  noted.  The  formation  of  dust  is  sensitive  to 
the  choice  of  the  propagation  rate  constant  since  the  density 
of  the  nth  member  in  the  polymerization  is  proportional  to 

)".  The  polymerization  model  and  its  consequences  will 
be  discussed  in  more  detail  in  an  upcoming  publication. 

Higher  silane  ions  ( n  >  2 )  can  have  large  fractional  den¬ 
sities,  particularly  at  high  silane  pressures  (F>0.3Torr)  and 
large  electrode  separations  {d>2  cm).’^  Therefore,  we  in¬ 
cluded  higher  silane  ions  having  n<6.  An  analogous  poly¬ 
merization  chain  to  that  proposed  for  the  higher  neutral  si¬ 
lanes  was  included  for  ions  and  will  be  discussed  in  Sec. 
Ill  E. 

Although  the  dissociative  attachment  cross  sections  for 
silane  and  disilane  producing  SiH2~  and  SiHf  are  not  par¬ 
ticularly  large,  the  residence  time  for  negative  ions  in  the 
discharge  can  be  long.  These  ions  were  therefore  included  in 
the  model.  For  conditions  typical  of  rf  discharges,  the  elec¬ 
tronegativity  [N,/{n,  +  N,),  IS  the  electron  density,  N, 
is  the  negative  ion  density]  is  0. 1-0.3.  For  dc  positive  col¬ 
umn  discharges,  where  the  electron  distribution  function 
more  closely  matches  the  dissociative  attachment  cross  sec¬ 
tions,  the  electronegativity  is  >0.9. 

B.  Reactions,  rate  constants,  and  cross  sections 

A  comprehensive  set  of  rate  constants  and  cross  sections 
was  compiled  for  the  electron  and  heavy  particle  reactions 
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occurring  in  glow  discharges  sustained  in  gas  mixtures  of 
SiH4,  Si2H6,  Hj,  and  Ar.  The  reactions  are  separated  into 
electron-neutral,  neutral-neutral,  ion-molecule,  electron- 
ion,  and  ion-ion.  Each  class  of  reactions  will  be  briefly  de¬ 
scribed  below. 

1.  Electron  impact  cross  sections  and  rate  constants 

The  electron-neutral  reactions  for  excitation,  dissocia¬ 
tion,  and  ionization  included  in  the  plasma  chemistry  model 
are  listed  in  Table  II.  In  addition  to  these  reactions. 


TABLE  11.  Electron  impact  excitation,  dissociation,  and  ionization. 


Process 

Rale  constant* 

Cross-section 

reference 

Collisions  with  silane 

e  +  SiH,  — SiH4(i;2,4)  +  e 

8.34(  -  9) 

16 

e-H  SiH4  — SiHjCvl.S)  +e 

8.74(  -9) 

16 

e-(-SiH4-SiH, -l-H-i-e 

l.59(  -  10) 

24-28 

e  +  SiH4-SiHj  +  H-l-H-)-e 

l.87(  -  II) 

24-28 

e  +  SiH4-SiH-l-H, +  H  +  e 

9.34(  -  1,2) 

24-28 

e-l-SiH4-SiHr  -l-  ‘H-)-2e 

3.64(  -  12) 

29 

e  +  SiH4-SiH/  +  H,  +  2< 

7.I9{  -  12) 

29 

e-)-SiH4-SiH*  4-H, +  H  +  2< 

l.39(  -  12) 

29 

e  +  SiH4  — Si"  +  2H,  +  2e 

I.19(  -  12) 

29 

e  +  SiH4-SiH,-  +  H 

2.50(  -  13) 

30,31 

e4-SiH4-SiHj  +Hj 

1.50(  -  13) 

30,31 

Collisions  with  disilane 

t  +  Si.H*  — Si;H,(v2,4)  +  e 

1.29(  -8) 

16 

t  +  Si]H4  — Si,H4(vl,3)  +  e 

1.06(  -8) 

16 

«  +  SijHt  — SiH4  +  SiHj  +  e 

2.86(  -  10) 

24 

t  +■  Hj  +  € 

1.23(  -  10) 

24 

e  +  Si,H4-Si.H,-  +  2e 

l.44(  -  II) 

29 

#-t-Si,H.-Si;H4"  +H, +  2e 

2.4S(  -  11) 

29 

e  +  Si,H,-SiH,'  +SiH,  +  2# 

l.75(  -  II) 

29 

e  +  Si,H,-SijH,"  -kH-kZe 

1.73(  -  ID 

29 

e4-Si,H,-Si:H,"  +  2Hj  +  2e 

1.15(  -  11) 

29 

e  +  Si,H,-SiH,-  +  SiHj 

9.4S(  -  14) 

16.32 

e-t-Si-H^-SiH,-  +  SiH4 

9.45(  -  14) 

16.32 

Collisions  with  hydrogen 

e  +  H,  — Hjtvl )  -t- 1 

8.8I(  -  10) 

33 

e  -i-  Hj  — H,(i»2)  -t-  e 

3.88(  -  11) 

33 

e  +  H,-H?-l-e 

1.28(  -  11) 

33 

e-l-Hj-H,*  -l-2e 

2.32(  -  ID 

34 

e+H,-H  +  H  +  e 

4.49(  -  12) 

35 

e+H-H*  +  e 

9.90(  -  12) 

36 

e  +  H-H"  +t 

2.62(  -  13) 

36 

Collisions  with  argon 

e  +  Ar— Ar*  +  e 

1.21(  -  12) 

37,38 

e  +  Ar— Ar**  +  e 

1.27(  -  12) 

39 

e  +  Ar— Ar*  +  2e 

10S(  -  ID 

34 

e  +  Ar*-Ar**  +  e 

8.07(  -  8) 

40 

e  +  Ar*  — Ar"  +  2e 

2.40(  -  9) 

41 

e  +  Ar**-Ar‘  +  2e 

7.02(  -  9) 

42 

*I.S(  —  10)  a  I.S  X  I0~  Rate  constants  have  units  of  cm^  '  and  were 
obtained  by  convolving  cross  sections  from  the  indicated  references  with 
the  electron  distribution  function,  as  described  in  the  text.  Rale  constants 
in  this  table  are  averaged  over  the  volume  of  the  reactor  and  over  the  rf 
cycle.  Rate  constants  for  electron  collisions  with  SiHt,  SijH,,  and  H,  are 
for  the  following  discharge  conditions:  SiH^Hi  »  I/I,  0.2S  mW  cm  ~ 

I  }.S6  MHz,  electrode  separation  2.S  cm,  0.2S  Torr,  and  substrate  tempera¬ 
ture  SOO  K.  The  rate  constants  for  collisions  with  Ar  are  for  an  Ar/SiH^  gas 
mixture  having  similar  discharge  conditions. 


elastic  (momentum  transfer)  collisions  were  included  in  the 
solution  for  the  electron  distribution  function  (see  below). 
Cross  sections  for  electron  impact  excitation,  ionization,  and 
dissociation  for  Ar  and  are  well  known  and  the  sources 
for  these  cross  sections  are  listed  in  Table  II.  Electron  impact 
cross  sections  for  silane  and  disilane,  though,  are  less  well 
known.  Only  a  subset  have  been  directly  measured;  total 
dissociation,’^  dissociative  ionization,’''  and  dissociative  at¬ 
tachment.’*''”  Electron  swarm  parameters  (drift  velocity 
and  net  ionization  coefficient),  though,  have  been  measured 
as  a  function  of  E  /N  (electric  field/number  density ) .  ”  With 
the  measured  cross  sections,  swarm  parameters  and  knowl¬ 
edge  of  the  approximate  shapes  of  cross  sections  for  other 
collisional  processes,  a  self-consistent,  though  not  unique, 
set  of  cross  sections  can  be  generated.  In  doing  so,  one  hy¬ 
pothesizes  a  set  of  cross  sections  and  uses  those  cross  sections 
to  compute  the  electron  distribution  function  and  swarm 
parameters.  The  computed  swarm  parameters  are  compared 
to  experiment.  The  process  is  iterated  until  agreement  is  ob¬ 
tained. 

This  method  for  generating  cross  sections  was  used  by 
Garscadden  et  a/.'’  and  Ohmori  ei  ai.''  for  silane,  and  by 
Hayashi"’  for  silane  and  disilane.  The  cross  sections  so  gen¬ 
erated  are  not  unique  and  depend  upon  the  method  used  to 
solve  for  the  electron  distribution  function.  The  cross  sec¬ 
tions  obtained  by  Garscadden  and  Hayashi,  for  example,  are 
quantitatively  different.  Both  sets,  though,  reproduce  swarm 
data  and  are  therefore  acceptable  solutions.  We  repeated  the 
process,  generating  self-consistent  sets  of  cross  sections  for 
silane  and  disilane  which  are  accurate  for  our  method  of 
solution.  We  used  as  a  basis  the  measured  dissociation,  disso¬ 
ciative  ionization,  and  attachment  cross  sections,  and  the 
momentum  transfer  and  vibrational  excitation  cross  sections 
suggested  by  Hayashi.  We  used  a  modified  form  of  the 
Monte  Carlo  simulation  described  in  Ref.  18  for  solving  the 
electron  distribution  function  in  a  dc  field.  It  was  necessary 
to  use  an  energy-dependent  anisotropic  elastic  scattering 
cross  section  to  obtain  agreement  with  swarm  parameters. 
The  cross  sections  for  silane  thus  derived  are  shown  in  Fig.  2. 


FIG.  2.  Electron  impact  cross  sections  for  SiH,  used  in  the  model  (see  Sec. 
II  B  I  and  Table  11  for  references). 


2535 


J  Appl.  Phys..  Vol.  63.  No.  8, 15  April  1968 


Mark  J.  Kushner 


2535 


The  agreement  with  swarm  measurements  is  shown  in  Fig. 

Electron  impact  rate  coefficients  were  obtained  as  a 
function  of  phase  during  the  rf  cycle  and  position  between 
the  electrodes  by  parameterizing  the  electron  kinetics  simu¬ 
lation  program  rfdis.  This  model  is  a  multidimensional, 
time-dependent  Monte  Carlo  particle  simulation  for  the 
electron  distribution  function  in  a  parallel  plate  capacitively 
coupled  rf  discharge,  and  has  been  described  in  detail  in  pre¬ 
vious  publications."*  Electron  impact  rate  coefficients  were 
obtained  by  convolving  the  cross  sections  with  the  electron 
distribution  generated  by  rfdis.  A  sample  of  rate  constants 
so  calculated  are  shown  in  Fig.  4.  High  threshold  energy 
processes,  such  as  ionization,  have  rate  constants  which  are 
significantly  modulated  during  the  rf  cycle,  a  consequence  of 
a  modulation  of  the  high-energy  tail  of  the  electron  distribu¬ 
tion  function.  The  high-energy  tail  is  “pumped”  at  the  peaks 
of  the  rf  cycle  near  the  electrodes,  a  result  of  the  oscillating 
sheaths.  Processes  with  lower  threshold  energies,  such  as 
vibrational  excitation,  have  rate  constants  which  are  more 
uniform  as  a  function  of  position  and  time,  as  lower  energy 
“bulk”  electrons  participate  in  these  processes.  Electron  im¬ 
pact  rate  constants  averaged  over  position  and  phase  of  the  rf 
cycle  appear  in  Table  II  for  typical  discharge  conditions,  as 
noted  in  the  table.  These  rate  constants  are  only  representa¬ 
tive;  the  manner  of  selecting  the  appropriate  coefficients 
during  execution  of  the  model  are  described  below. 

For  each  initial  gas  mixture  and  plate  separation,  rfdis 
was  parameterized  while  using  different  values  of  the  applied 
rf  voltage  and  electric  field  in  the  bulk  plasma.  For  each  run 
ofRFDis,labeled7,  we  computed  the  quantity/}^  =  (F/n,  )^, 
where  P  is  the  spatially  averaged  power  deposition  per  unit 
volume  and  n,  is  the  spatially  averaged  electron  density.  A 
set  of  time  and  spatially  dependent  electron  impact  rate  con¬ 
stants  was  generated  and  labeled,  k,  {pj ) ,  where  i  denotes  the 
electron  impact  collision.  A  look-up  table  was  constructed  of 
these  rate  constants  using  pj  as  the  parametric  variable. 
When  executing  the  plasma  chemistry  model,  we  used  power 
deposition  per  unit  volume,  F,  as  a  scaling  parameter.  Know- 


FIG.  }.  Net  ionization  coefficient  (ionization  attachment)  Tor  dc  dis¬ 
charge!  in  silane  as  a  function  of  electric  fleld/number  density.  Theoretical 
results  are  as  generated  with  the  cross  sections  from  Fig.  2.  The  experimen¬ 
tal  results  are  from  Shimozuma  and  Tagashira  (Ref.  43). 


(oi  —  Ar'*  2e 


=  4  5  X 


(Di  e-SiH,  — neutral  fragments  t-e 
(k„„,  =  2.0xi0AmV‘l 


(cle*SiM4  -S'M,i,|*e 
(•■moi*  4  5  *  I0'*cm^s'') 


FIG.  4.  Typical  electron  im¬ 
pact  rale  constants  as  used  in 
the  plasma  chemistry  model 
plotted  as  a  function  of  posi¬ 
tion  between  parallel  plate 
electrodes  and  time,  ex¬ 
pressed  as  fraction  of  the  rf 
cycle;  (a)  electron  impact 
ionization  of  argon,  (b)  neu¬ 
tral  dissociation  of  silane, 
and  (c)  ribrationa)  excita¬ 
tion  of  silane.  As  the  thresh¬ 
old  energy  increases  from  vi¬ 
brational  excitation  to 
ionization,  the  uniformity  of 
excitation  (position  and 
time)  decreases. 


ing  the  instantaneous  electron  density,  n,(f),  the  electron 
impact  rate  constant  k,  was  obtained  by  interpolating  the 
table k/[(P/n,)j]  using  the  parameter P /n,  (I).  The  power 
includes  contributions  from  ions  which  are  accelerated 
through  the  sheaths  and  strike  the  electrodes.  For  the  refer¬ 
ence  conditions  presented  in  Sec.  Ill  (SiH4  2S0  mTorr,  25 
mW  cm“’),  <  5%-10%  of  the  power  is  dissipated  in  this 
fashion. 

2.  Neutral-neutral  reactions 

The  neutral  chemistry  occurring  in  silane  plasmas  has 
been  the  topic  of  numerous  previous  investigations  and  is 
summarized  in  Refs.  2  and  7.  We  have  benefited  by  those 
works  and  others,  and  have  collated  the  suggested  reaction 
pathways  into  a  global  set  of  reactions.  The  neutral-neutral 
reactions  included  in  the  model  are  listed  in  Table 
III  i2.2s.26.28.M-62  i^^  ajjjfjQn  {(,  yeactions  suggested  by 
previous  investigators,  we  included  polymerization  reac¬ 
tions  leading  to  the  formation  of  higher  silane  radicals  and 
molecules,  and  higher  silane  ions.  These  reactions  are  dis¬ 
cussed  in  Sec.  II  A  and  III  E.  We  also  included  homogen¬ 
eous  pyrolysis  reactions  of  silane  and  disilane.  The  substrate 
temperatures  typically  used  in  PECVD  are  <300  K,  and 
based  on  these  temperatures,  pyrolysis  should  not  be  impor¬ 
tant.  For  moderate  power  deposition  {  >  25  mW  cm^)  and 
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TABLE  III.  Neutral  heavy  particle  reactions. 


Reactions  with  silane  and  disilane 
SiH4  +  H-SiH,  +  H, 

SiH4  +  Si-SiH,+SiH, 

SiH4  +  SiH^SiiH, 

SiH4  +  SiH-Si,H,  +  H, 

SiH4  -t-  SiH2-<Si;H; 

SiH4  +  SiH,-Si,H,  +  Hj 
SiH4  +  SijH4-.Si,H, 

SiH4  +  SijHj-SijH,  +  SiH, 

SijH»  +  SiH,-SiH4  +  SijH, 

SijH.  +  SiHj-.Si,H, 

SijH»  +  SiH-Si,H, 

SijH*  +  H^SiH,  +  SiH4 
SijH,  +  H-Hj  +  Si,H, 

Si,Hr-Si,H4  +  Hj 
SijH; +  M-SijH»  +  M 
SijHr-.Hj  +  Si,H4 
SijHr*-.SiH4  +  SiHj 
SijH;* +  M -SijH* +  M 


Reactions  of  radicals  with  radicals,  higher  silanes,  and  hydrogen 
SiH,  +  H-SiHj  +  Hj 
SiH, +  SiH2-Si,H, 

SiH, +  SiH,-SiH, +  SiH4 
SiH,  +  SiH, -SijH;* 

SiH,  +  Si2Hj-Si,H, 

SiH2  +  Si-Si,H, 

SiH2  +  SiH-Si,H, 

SiHj  +  H2-SiH4 
SiH, +  Si,H,-.Si4H,« 

SiH2  +  H-SiH  +  H2 
SiH, +  H -SiH, 

SiH  +  Hj-SiH, 

Si  +  H,-SiH2 
SijH4  +  H2-Si,H* 

Si:H4  +  H2-SiH4  4-SiH, 

SijH,  +  Hj-Si.H, 

SijHj  +  H,-Si,H4 
SijHj  +  H-SijH, 

SijH,  +  Si,H,-Si4H,„ 

SijH4  +  SijH*-Si4H,o 

SiH  +SijH,-Si,H, 

Si,H,  +  H-Si,H, 

SiH  +  SijH.-SuH, 

SiH, +  Si,H,-Si.H,  +  H. 

Si4H,  +  H-Si4H,„ 

SiH  +  Si4H,„-Si,H„ 

SiHj  +  Si4H,„-Si,H„ 

SiH, +  Si4H  ,0-Si, H„  +  H, 

Si,H,,  +  H-Si,H,j 


Excitation 
lAr’,Ar**l  + 
'Ar*,Ar**]  + 
(Ar*.Ar**)  + 
(Ar*.Ar**l  + 
[Ar*.Ar**l  + 
(Ar*,Ar**)  + 
(Ar*,Ar**l  + 
(Ar*,Ar**)  + 
Hf +  H*-H,* 
Ar**-Ar* 
H?-H, 
H*-H 


■  transfer  and  radiative  decay 
Hj-H  +  H  +  Ar 
SiH4-SiH,  +  H  +  Ar 
SiH4-SiH,  +  H  +  H  +  Ar 
SiH.-SiH, +  H  + Ar 
SiH,-SiH  +  H  +  Ar 
SiH-Si  +  H  + Ar 
SijH.-SijH4  +  H  +  H  +  Ar 
SijHj-SijHj  +  H  +  H  +  Ar 
'  +  Hj  +  e 


Rate  constant* 

Reference* 

2.6>(  -  12) 

26,44-46 

S.33(  -  13) 

47 

2.S0(  -  12) 

25,46,48 

l.70(  -  12) 

25,46,48 

1.00(  -  II) 

12.25,26,47,49-51 

l.7«(  -  15) 

47 

l.00(  -  II) 

est.  (52) 

5.00{  -  13) 

est. 

I.QO(  -  12) 

est.  (53) 

l.20(  -  10) 

12.50-52 

l.00(  -  II) 

PC 

).1I(-I2) 

53 

2.I6(  -  12) 

49,53 

5.00(6)s-' 

est.  (28) 

l.00(  -  10) 

est.  (54) 

2.30(7)  s-' 

est.  (28) 

2.30(7)  s-' 

esL  (28) 

2.00(  -  10) 

est.  (54) 

l.00(  -  10) 

est.  (44) 

3.77(-I3) 

46 

7.00(  -  12) 

12,26-28,44 

l,00(  -  II) 

12,26,44 

l.00(  -  II) 

est.  (52) 

4.53(  -  13) 

47 

7.22(  -  13) 

47 

2.00(  -  13) 

50-52,55 

1.00(  -  II) 

PC 

7.96(  -  13) 

46 

Ml(  -  12) 

46 

l.98(  -  12) 

47 

6.S9(  -  12) 

47 

5,33(  -  13) 

47 

3.56(  -  9) 

47 

1.70(  -  12) 

47 

l.•W(  -  11) 

47 

4.94(  -11) 

47 

1.00(  -  II) 

PC 

l.00(  -  ID 

PC 

l.00(  -  ID 

PC 

l.00(  -  ID 

PC 

1.00(  -  ID 

PC 

l.00(  -  ID 

PC 

1.00(  -  ID 

PC 

1.00(  -  ID 

PC 

l.00{  -  11) 

PC 

l  .00(  -  11) 

PC 

1.00(  -  11) 

PC 

7.00(  -  1 1) 

56 

1.40{  -  10) 

c,  56-58 

2.60(  -  10) 

c,  56-58 

l.00(  -  10) 

c,  56-58 

1  00(  -  10) 

c,  56-58 

1.00(  -  10) 

c,  56-58 

().00(  -  10) 

c,  56-58 

6.60(  -  10) 

c,  56-58 

1.00(  -  11) 

est. 

1.00(7)  s-' 

est. 

1.00(7)  s-' 

est. 

2  00(7)  s  ■ 

est. 
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TABLE  III.  {continued). 


Process 

Rate  constant* 

Reference* 

Vibmtioiul  relaxation  (see  text) 

SiH4(i;2.4)  +  SiH4-SiH4  +  SiH^ 

1.92(  -  13) 

39,60 

SiHiCvM)  +  Ar— At  -I-  SiH4 

1.86(  -  14) 

60,61 

SiH4(i;2.4)  +  H,-Hi  -k  SiH, 

3.05(  -  12) 

60-62 

SiH4(«2.4)  +  SijHt-SijH,  -1-  SiH* 

S.7l(-n) 

60-62 

SiH4(t/l.3)  +  SiH4-SiH«  +  SiH4(i»2.4) 

6.08(  -  11) 

60-62 

SiH4(Bl.3)  +  Ar-Ar  +  SiH4(B2.4) 

5.37(  -  12) 

60-62 

SiH4(i>1.3)  +  Hj-Hj  -(-  SiH4(  1-2.4) 

l.53(  -  10) 

60-62 

SiH4(»l,3)  -1-  Si.H*-Si,H*  -|-  SiH,(B2.4) 

3.00(  -  10) 

60-62 

Si,H»(i-2,4)  -1-  SiH4-SiH4  -t-  Si.H* 

6.40(  -  12) 

60-62 

SijH»(i-2.4)  +  Ar-Ar  +  SijH* 

I.09(  -  12) 

60-62 

Si,H,(ti2,4)  -1-  Hj-Hi  -1-  SijH. 

1.80(  -  10) 

60-62 

Si,H»(i;2.4)  +  SijH^-Si.H.  +  SijH» 

1.13(  -11) 

60-62 

SijH,(i>l,3)  -1-  SiH4-SiH4  +  SijH*(u2.4) 

3.00(  -  10) 

60-62 

SiiH^Cul.S)  +  Ar-Ar  +  Si2H4(v2.4) 

3.00(  -  10) 

60-62 

SijH»(i-l.3)  -1-  H,-Hj  ^  SijH,(B2.4) 

3.00(  -  10) 

60-62 

SijH,(i>l,3)  +  SijH^-SijH,  +  SijH*(B2.4) 

3.00(  -  10) 

60-62 

SiH4(B2.4)-SiH4 

2.30(00) 

est.  (62) 

SijH»(i(2.4)-Si,H» 

2.30(00) 

est.  (62) 

Homogeneous  pyrolysis 

SiH4-SiHj  +  Hj 

3.0(12)exp(  -  32.2  kcal/Zir)  s"* 

47 

SiH4-SiH,  +  H 

3.7(  13)exp(  -  93.0 kcal/AD  s" ' 

47 

Si,H4-SijH4  +  Hj 

2.3(  14)exp(  -  52.2  Itcal/kD  s"' 

47 

Reactions  leading  to  formation  of  DUST  {X  +  T— DUST). 

(Uniform  rate  constant  1.0(  —  11)  cm*s*' 

.  See  text) 

SiH.  +  (Si,H„3i,H,:l 

Si4H,-l-lSi4H*Si4H,o] 

SiH,  +  DUST 

SuH,-!- (Si,H.,.Si,H,jJ 

Si,H.  -1-  (Si4H«Si4H,ol 

Si4H,  +  DUST 

Si,H.  +  (Si,H„.Si,H„l 

Si,H„  +  lSi,H,„Si,H,:l 

SijH,  -1-  DUST 

Si,H„  +  DUST 

SiaH^  +■  [SijH7,Si5H4l 

Si,H,  +  (Si4H*Si4H„l 

Si,H,+  (Si,H„.Si,H„l 

Si,H,  +  DUST 

*I.SO(  —  10)  si. SOX  10''°.  Rateconsunu  have  units  of  cm^  s~ '  unless  otherwise  noted. 

‘PC  denotes  that  the  reaction  is  hypothesized  as  being  part  of  the  polymerization  chain  leading  to  formation  of  DUST  and  has  an  estimated  rate  constant  ( see 
text). 

‘Rate  constant  estimated  by  analogy  to  CH^,  and  their  radicals. 


electrode  separations  (d  >  3  cm),  though,  the  centerline  gas 
temperature  can  exceed  the  substrate  temperature  by  many 
hundreds  of  degrees  (see  Sec.  II  C).  For  these  conditions, 
pyrolysis  will  occur  in  the  bulk  plasma.  Rates  for  homogen* 
eous  pyrolysis  of  silane  and  disilane  were  obtained  from  Col- 
trin,  Kee,  and  Miller.'*’ 

A  subset  of  the  rate  constants  listed  in  Table  III  was 
estimated  and  they  are  so  noted  in  the  table.  The  major 
classes  of  rate  constants  that  were  estimated  are  those  for  the 
polymerization  chain  (see  above).  Penning  reactions  of  elec¬ 
tronically  excitated  argon  with  silane  and  disilane,  and  colli- 
sional  deactivation  of  SiH4(i»l,3;  b2,4)  and  Si,Hft(pl,3; 
p2,4).  For  the  Penning  reactions,  (e.g.,  Ar*  +  SiH4— SiHj 
H  Ar),  we  used  the  analogous  rates  for  Penning  reac¬ 
tions  of  argon  with  methane  (CH4  Rate  constants  for 
vibrational  relaxation  of  silane  and  disilane  were  also  ob¬ 
tained  by  analogy  to  methane.  In  methane,  rapid  vibrational 
relaxation  proceeds  by 
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CH4(i;1,3)-CH«(i;2,4)-CH4(p0,0).  (3) 

Rate  constants  for  the  vibrational  relaxation  of  CH4  by  Ar, 
Hj,  CH4,  and  other  molecules  are  available,  as  well  as  for  the 
vibrational  relaxation  of  silane  by  silane. By  comparing 
the  rate  of  self-relaxation  of  silane  to  that  of  methane,  one 
obtains  a  scaling  parameter  representing  the  ratio  of  the 
average  mutual  interaction  potential,  approximately  8.9. 
This  value  was  used  to  scale  the  vibrational  relaxation  rates 
of  methane  with  various  collision  partners  to  silane.  The  re¬ 
sults  are  shown  in  Table  Ill. 

Rate  constants  were  also  estimated  when  quoted  or 
measured  values  varied  over  a  large  range.  In  this  class  of 
reactions  are  the  important  processes: 

H-(-SiH4-SiHj-(-Hj,  (4) 

-  .  Mark  J  Kushner  2538 


SiH2  +  SiH4-Si2H?*  +  M-SijHft  (5a) 

i 

Si2H4  +  H, 

r 

SiHj  +  SiHj  -  SijH?*  +  M  -  SijH^  ( 5b ) 

-.SiH2  +  SiH4.  (5c) 

Reactions  (4)  and  (5a)  are  important  because  they  are  first- 
order  reactions  of  the  dissociation  products  of  silane  with 
silane.  The  former  hydrogen  abstraction  reaction  is  effec¬ 
tively  a  radical  multiplication  reaction.  A  radical  which  does 
not  directly  contribute  to  film  growth,  H,  is  exchanged  for  a 
radical  that  will  be  incorporated  into  the  film,  SiHj.  The 
latter  reaction  is  important  for  the  opposite  effect;  a  radical 
which  can  be  directly  incorporated  in  the  film,  SiH,,  is  ex¬ 
changed  for  either  a  saturated  molecule  ( Si2H(,)  which  is  not 
incorporated  into  the  film,  or  a  molecule  ( Si2H4 )  which  like¬ 
ly  has  a  lower  sticking  coefficient.  Although  reactions  (5c) 
and  (5d)  are  second  order,  they  are  important  because  the 
most  probable  branches  reduce  by  at  least  half  the  number  of 
radicals  which  can  contribute  to  film  growth.  Obser\'a(ions 
that  film  growth  increases  less  than  linearly  with  increasing 
discharge  power  could  be  attributed  to  this  reaction,  whose 
rate  scales  as  the  square  of  the  rate  of  dissociation  of  silane.*^^ 
The  rate  constants  for  these  reactions  were  studied  parame¬ 
trically,  as  reported  in  Ref.  26.  The  values  for  the  rate  con¬ 
stants  listed  inTable  III  are  those  resulting  from  this  study. 


3.  Ion-neutral 

The  ion-molecule  reactions  included  in  the  model  are 
listed  in  Table  iv.''*  '’*"’’  Rate  constants  or  cross  sections 
have  been  measured  for  Si„ H,;  H/,'”'  and  Ar""^’® 

collisions  with  silane  and  disilane  as  a  function  of  ion  energy. 
Ion-molecule  reactions  were  placed  in  endothermic  and  ex¬ 
othermic  categories  based  on  the  average  gas  temperature. 
The  endothermic  reactions  typically  have  energy  defects  of  a 
few  eV,  Those  reactions  classed  as  e.xothermic  were  included 
in  the  model  at  all  locations  in  the  discharge.  Those  reactions 
classed  as  endothermic  were  included  only  in  the  sheaths. 
Positive  ions  in  the  sheaths  are  accelerated  to  energies  in 
excess  of  the  energy  defect,  thereby  allowing  reactions  which 
are  endothermic  in  the  bulk  plasma  to  occur  in  the  sheaths. 
The  endothermic  ion-molecule  reactions  for  which  this  ap¬ 
plies  are  labeled  by  in  Table  IV. 

Ion  molecule  reactions  between  Si,  !!„;■  ions  and  silane 
or  disilane  generally  add  silicon  atoms  to  the  ions  for  n  <  5-6. 
For  example,  the  major  branches  for  reactions  of  SijHj^ 
with  silane  and  disilane  are''* 

SijHj"  -HSiH4^Si,H,*.  (6a) 

SijH,*  -f  (6b) 

Such  addition  reactions  have  nearly  gas  kinetic  rate  con¬ 
stants  (  =:  10”  "'cm^  s” ').  Asa  result,  higher  silane  ions  can 
be  generated  quickly  in  silane  plasmas  from  the  products  of 
the  dissociative  ionization  of  silane,  primarily  SiH,"  and 
SiH,‘ . 

In  analogy  to  the  neutral  polymerization  model  de¬ 


scribed  in  Sec.  II  A,  we  hypothesized  a  “clustering”  chain 
for  higher  silane  ions.  The  clustering  reaction  of  lower  silane 
ions  with  higher  silane  saturated  molecules  has  the  form  of 

"F  Si„H2(,„4. 1)  -»Si,„4.  |H2'^,„4.  I)  +  H,  (7a) 

SiH2*  +Si„H2<,,.^.,-Si„^,H2t„4..  +H2,  (7b) 

SiH,^  +Si„H2,„,„-SL^,H2t„,n  +  H,  +  H. 

(7c) 

Rate  constants  for  reactions  in  the  clustering  chain  were  esti¬ 
mated  to  be  10“ '  ‘  cm’  s“ ',  if  not  otherwise  available. 

Recent  work  by  Mandich,  Reents,  and  Jarrold’'*  has 
shown  that  ion-molecule  reactions  of  t  ie  form  of 

Si,D,:  -l-SiD4^Si,^,D,:^2 +D2  (8) 

form  a  chain  which  may  terminate  with  n  =  4  due  to  an 
activation  energy  barrier.  In  a  plasma,  a  significant  activa¬ 
tion  energy  is  available  in  the  form  of  excited  states  of  the 
reactants  and  which  may  aid  in  breaching  any  activation 
energy  barrier  in  the  chain.  Since  silane,  as  the  feedstock,  has 
the  highest  number  density  of  any  reactant,  sequential  ion- 
molecule  reactions  with  silane  could  potentially  and  rapidly 
lead  to  higher  silane  ions  if  the  chain  propagates.  This  topic 
will  be  discussed  further  in  Sec.  Ill  £. 

Ion-molecule  reactions  of  H2*  and  Ar"^  with  silane  (and 
disilane)  usually  result  in  charge  exchange  to  silane  and  re¬ 
duce  the  number  of  H  atoms  in  the  product  ions.  (The  SiH4"^ 
ion  has  only  been  observed  under  near  collisionless  condi¬ 
tions  and  is  believed  to  be  unstable  in  typical  discharges. ) 
The  differences  in  the  ionization  potentials  of  hydrogen 
(15.4  eV)  and  argon  ( 15.8  eV)  with  the  appearance  poten¬ 
tials  of  ions  from  silane  ( 1 1 .6  eV )  and  disilane  (9.9  eV )  are  a 
reservoir  of  energy  for  dissociative  branches.  Typical  reac¬ 
tions  are'^*’’*-’” 

H2*  +  SiH4-.SiH3*  -t-  H2  -(-  H,  (9a) 

Ar" -l-ZHj-fH.  (9b) 

Ion-molecule  reactions  of  Ar"^  with  silane  have  unusually 
small  rate  constants  ( l-4x  10“  ”  cm*  s“ ')  when  compared 
to  the  same  reactions  for  ,  Ne'^ ,  and  He  *  The  frac¬ 
tion  of  Ar’"  ions  in  Ar/SiH4  gas  mixtures  will  therefore  be 
larger  than  the  analogous  fractions  in  gas  mixtures  contain¬ 
ing  hydrogen,  neon,  or  helium. 

4.  Ion-ion  and  electron-ion  reactions 

The  ion-ion  and  electron-ion  reactions  and  rate  con¬ 
stants  included  in  the  model  are  listed  in  Table  V.’’"”*  Nega¬ 
tive  ions  are  formed  in  silane  plasmas  by  the  dissociative 
attachment  reactions 

e-f-SiH4^SiH3-  H  (10a) 

-SiHj+H,.  (10b) 

The  cross  sections  for  dissociative  attachment  to  silane  and 
disilane  are  small  (<2x  10“  "*  cm’),'*’”  as  are  the  rate  con¬ 
stants  (10“'’-10“'’  cm’s”').  Negative  ions,  though,  are 
not  lost  rapidly  to  the  walls  due  to  the  negative  space  charge 
there.  They  are  lost  dominantly  by  electron  impact  detach¬ 
ment  and  ion-ion  neutralization  reactions.  The  electronega¬ 
tivity  ( fraction  of  negative  ions  compared  to  the  total  density 


2539 


J  Apol  Phys..  Vol  63.  No  8.  15  April  1988 


Mark  J  Kushner 


2539 


TABLE  IV.  Ion-molecule  reactions. 


Process 

Rate  constant* 

Reference* 

SiH,*  +  SiH^-SijH"  -1-  H,  -h  H,  -l-  H, 

0.13(  -  10)  E 

14,65 

SiH,*-  +  SiH^-SijH/  -(-  Hj  -t-  Hj  -1-  H 

0.25(  -  10)  E 

14,65 

SiH,*  -1-  SiH4-Si,H/  -f-  Hj  Hj 

0.07(  -  10) 

66 

SiH,*  -1-  SiH«  -1-  SijHr  +  H, 

0.25(  -  10) 

14,66 

SiH,*  -t-  SiH.  -1-  M-SijH/  +  M 

I.g0(  -26)  cm*s-' 

67 

SiHj'  -(-Si,H«-SiiH,*  -l-SiH. 

5.20(  -  10) 

68 

SiH,^  -1-  SijH.-Si,H/  -(-  H,  -l-  H 

3.S0(  -  10) 

68 

SiH,*  -1-  Si,H,-Si.H,*  -1-  H,  -|-  H 

1.00(  -  11) 

PC 

SiH,*  -(-  Si.H,o-Si,H,;  -1-  H,  H 

I.00(  -  11) 

PC 

SiH,*  +  SiH. -SiH,*  -l-SiH, 

1.07(  -9) 

66 

SiH,*  +  SiH.-Si,H*  -l-  H,  -l-  H,  -i-  H 

0.!3(  -  10)  E 

65 

SiH,*  SiH.-Si,H,*  -l-  H,  -t-  H, 

0.55(  -  10) 

66 

SiH,*  -h  SiH.-Si,H,*  -t-  H,  -t-  H  H 

0.70(  -  10)  E 

65 

SiH,*  +  SiH.-Si,H,*  -H  H,  -l-  H 

0.12(  -  10)  E 

14,65 

SiH,*  -(-SiH.-Si,H,*  -(- H, 

2.50(  -  10) 

14,65 

SiH,*  -i-SiH.-Si,H,*  -l- H 

0.08(  -  10)  E 

65 

SiH,*  -1-  H,-SiH,*  -1-  H 

I.OK  -  10)  E 

69 

SiH,*  -1-  Si,H,-Si,H,*  -f-  SiH.  +  H, 

1.40(  -  10) 

68 

SiH,*  -1-  Si,H,-Si,H,*  -1-  SiH.  +  H, 

3.10(  -  10) 

68 

SiH,*  -hSi,H*-Si,H,*  -(-SiH, 

3.60(  -  10) 

68 

SiH,*  -i-Si,H*-Si,H4*  -1- H, 

5.62(  -  10) 

68 

SiH,*  -1-  Si,H,-.Si.H,*  +  H, 

l.00(  -  11) 

PC 

SiH,*  +  Si.H,o-Si,H,J  -1-  H, 

l.00(  -  11) 

PC 

SiH*  -(-SiH.-SiH,*  -l-SiH, 

0.60(  -  10)  E 

14,65 

SiH*  SiH.-Si,H*  -(-  H,  -(-  H  -1-  H 

1.20(  -  10)  E 

65 

SiH*  -I-  SiH.-Si,H*  -1-  H,  -(-  H, 

0.70(  -  10) 

14,66 

SiH(i.)*  -(-Si,H*  -(-  H,-t-H, 

1.20(  -  9) 

90d,e 

SiH*  +  SiH. -Si, H,*  -1-  H,  -1-  H 

0.52(  -  10)  E 

65 

SiH*  +  SiH. -Si, H,*  -i-H, 

2.80(  -  10) 

66 

3.20(  -  10) 

90d 

SiH*  -(-  SiH.  -I-  M-Si,H,*  +  M 

l.30(  -  26) 

67 

SiH*  -t-  H,-SiH,*  -1-  H 

l.84(  -  11)  E 

69 

SiH*  -1-  Si,H»-Si,H*  +  SiH.  -1-  H, 

1.10(  -  10) 

68 

SiH*  -(•  Si,H»-Si,H,*  +  SiH. 

3.30(  -  10) 

68 

SiH*  -t-Si,H»-Si,H,*  -h  H 

3.50(  -  10) 

68 

SiH*  -t-Si,H,-Si.H,*  +H 

1.00(  -  11) 

PC 

SiH*  +Si.H,o-Si,H,;  -l- H 

1.00(  -11) 

PC 

Si*  -(-SiH.-Si,H,*  -(- H, 

4.g0(  -  10) 

66 

6.90(  -  10) 

79d 

Si*  +SiH.-Si,H,*  -l- H 

0.15(  -  10)  E 

14,65 

Si* -(- SiH.-Si,H* -1- H, -h  H 

0.1J(  -  10)  E 

14,65 

Si*  -(-SiH. M-Si,H.*  -(-M 

1.40(  -26)cm‘s-' 

67 

Si*  +  Si,H»-Si,H,*  +  SiH. 

7.80(  -  10) 

68 

Si*  +  SijH*— Si,H,* 

2.40(  -  10) 

68 

Si,H,*  -(-  SiH.-Si,H**  -(-  H,  -(•  H,  -(-  H 

5.00(  -  11) 

14 

Si,H,*  +  SiH.-Si,H,*  -(■  H,  -(-  H, 

5.001  -  11) 

14 

Si,H,*  -(-  SiH.-Si,H  *  -(-  H,  -(-  H 

5.00(  -  11) 

14 

Si,H.*  +SiH.-Si,H*  -(- H, 

5.20(  -  11) 

14,67 

Si,H,*  -f  SiH.-Si,H.*  +  H 

5.00(  -  ID 

14,67 

SijH,* -(-SiH.-Si,H.* 

3.76(  -  11) 

14,67 

Si,H,*  -(-SiH.-Si,H.*  -|- H, 

3.20(  -11) 

79d,e 

Si,H*  -(-  SiH.-Si,H,*  -(-  H,  +  H, 

3.10(  -  10) 

90d,e 

SijH*  -(-SiH.-Si,H.* 

l.08(  -  10) 

14,67 

Si,H.*  -(-Si,H»-Si,H,*  -(■  Si,H, 

5.80(  -  10) 

68 

Si,H,*  -(-Si,H,-Si,H4*  -(- SiH. -(- H 

0.86(  -  10) 

68 

Si;H.*  -i-Si,H«-Si,H.*  -(-SiH. 

0.65(  -  10) 

68 

Si,H.*  +Si,H4-Si.H,*  -(-H, 

1.40(  -  10) 

68 

Si,H,*  -f-Si,H»-Si,H.*  -(-SiH, 

1.87(  -  10) 

68 

Si,H,*  -(-  Si,H»-Si.H,*  -(-  H 

3.34(  -  10) 

68 

Si,H,*  -(-Si,H4-Si,H.*  -(-SiH, 

I.10(  -  10) 

68 

Si,H,*  -(-  Si,H»-Si.H,* 

2.95(  -  10) 

68 

Si,H*  -(-Si,H»-Si,H,*  -(-SiH 

0.30(  -  10) 

68 
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TABLE  IV.  (continued). 


Process 

Rate  constant* 

Reference" 

SijH*  +  SijH.-Si.H,’^ 

4.90(  -  10) 

68 

Si.H,"  -hH, 

1.00{  -  11) 

PC,c 

SljH,*^  -l-SiH.-Si.H.’^  -I-Hj 

1.70(  -  10) 

79d,e 

Si,H.*  -t-SiH.-SijHiS 

1.00(  -  13) 

79d,e 

-)-SiH,-Si,H,*  -l-H, 

1.30(  -  10) 

90d.e,f 

H”  -t-SiH^-SiH,"^  -t-Hj 

3.00(  -  10) 

est. 

H,’  -1-  SiH.-Si*  -1-  Hj  -k  H,  -)-  H, 

3.66(  -  11) 

69 

+  SiH,-SiH*  -)-  H,  -)-  Hj  -(-  H 

3.66(  -  11) 

69 

H,*  +  SiH,-SiH,*  -(-  H,  -1-  Hj 

6.39(  -  11) 

69 

Hj*  -(-SiH^-SiH,*  -l-Hj-hH 

6.23(  -  10) 

69 

-t-  SiH^-SiHj  +  H,* 

1.83(  -  ID 

69 

H,’  -)-SiH,-SiH,*  -l-Hj-l-Hj 

5.16(  -  10) 

69 

Ar  4*  SiH4  -•  Si  +  Hj  -+■  -f-  Ar 

l.lg(  -  11) 

64,70 

Ar*  -(-  SiH,-SiH*  -1-  H  -1-  Hj  -1-  Ar 

4.20(  -  12) 

64,70 

Ar*  -1-  SiH,-SiH,*  -t-  Hj  -1-  Ar 

1.66(  -  12) 

64,70 

Ar*  +  SiH,-SiH,*  +  H  +  Ac 

2.40(  -  12) 

64,70 

Ar*  +  Si,H.-Si,H*  -l-  Hj  -)-  H,  -)-  H  -l-  Ar 

1.80(  -  10) 

64,70 

Ar*  -1-  SijH,— SijHj*  -H  Hj  -1-  H,  -f  Ar 

6.30(  -  U) 

64,70 

Ar*  -1-  Si,H„-Si:H,*  -1-  H,  -f  H  -1-  Ar 

3.30(  -  11) 

64,70 

H,*  -(-Hj-H,"  -t-H 

2.10(  -9) 

71 

Ar*  -|-H,-ArH*  -l-H 

7.40(  -  10) 

72 

H,*  -1-  Ar  — ArH*  -i-  Hj 

l.00(  -  11) 

73 

ArH*  -1-  H,-H,*  -1- Ar 

I.30(  -9) 

73 

'I.S(  —  10)  =  l.Sx  10*'°.  Rate  constants  have  units  of  cm’ s  ' '  unless  otherwise  noted.  The  notation  "E”  by  the  rate  constant  specifies  that  this  reaction  is 
endothermic  in  the  bulk  plasma  and  occurs  only  in  the  sheaths  (see  Sec.  II  B  3). 

"PC  denotes  that  the  reaction  is  hypothesized  as  being  part  of  an  Ion  clustering  chain  and  has  an  estimated  rate  constant  (see  Sec.  II  B  3). 

‘See  Sec.  Ill  E. 

’This  most  recent  value  of  reaction  rate  constant  for  analogous  deuterated  species  was  not  available  at  the  time  of  formulation  of  the  model.  The  rate  constant 
is  listed  for  reference. 

This  reaction  was  not  explicitly  included  in  the  model. 

'Structural  isomeric  reactant  ion,  estimated  abundance  0.73. 


TABLE  V.  Charged  particle-charge  particle  reactions. 


Reaction 

Rate  constant' 

Reference 

Electron-ion  recombination" 

e-i-SiH,'  -SiH-l-H 

l.b9(  -7) 

73 

e-i-SiH,*  -SiHj -l-H 

1.69(  -7) 

73 

e-l-SiH*-Si-l-H 

1.69(  -7) 

75 

e-i-Si,H*-Si-i-H-i-H 

1.69(  -7) 

75 

e  +  SijHj'  —Si  -i-  Si  -1-  Hj 

1.69(  -7) 

73 

e  +  SijH,;  (m>3)— SijHi^.  ,,  -1-  H 

1.69(  -7) 

73 

e-i-Si,H.-  -SijH, -fSiH, 

1.69(  -7) 

75 

e  -t-  SijH,*  —  SijH,  -i-  SijH, 

1.69(  -7) 

75 

e  +  Si,H|J  — Si,H,  +  SijHj 

1.69(  -7) 

75 

e  +  Hj*  -H*  +  H 

3.66(  -  8) 

76 

e-l-H,-  -HJ-i-H 

9.73(  -8) 

76 

e-l-  H*-H* 

2.62(  -  13) 

76 

e  +  Ar*  —  Ar** 

8.13(  -  13) 

77 

Ion-ion  neutralization 

SiH;  -i-Si,H„-SiH.  -l-Si,H„ 

3.00(  -  7) 

est.  (78,79) 

SiH.-  +M*-M-i-SiH, 

1.00(  -7) 

est.  (78,79) 

'I.3(  —  10)  a  1. 3  X  10* '°.  Rate  constants  have  units  of  cm’ s  *' unless  oth¬ 
erwise  noted. 

"Rate  constants  for  dissociative  recombination  of  Si,  H.,'  were  estimated  by 
analogy  to  methane. 


of  negative  ions  and  electrons)  of  silane  plasmas  can  there¬ 
fore  be  many  tenths  to  near  unity.  We  estimated  that  ion-ion 
neutralization  reactions  between  silane  ions  (e.g.,  SiHf 
+  SiHj*  —2  SiHj)  have  uniform  rate  constants  of  5  X  10"^ 
cm'  s“',  whereas  neutralization  between  unlike  ions  (e.g., 
SiHj*  -I-  Ar’"  —  SiHj  -i-  Ar)  have  rate  constants  of  1 X  lO*^ 
cm’s-'.’'*’'' 

Electron-ion  recombination  in  the  bulk  plasma  of  glow 
discharges  typically  occurs  by  radiative  recombination  (e.g., 
e -I- A"*"  —  A*),  collisional  radiative  recombination  (e.g., 
e -I- e -I- A*"  —  A* -I- e),  or  dissociative  recombination  of 
molecular  ions  [e.g.,  e  +  Hj’"  — H  -|-  H(i;)].  Typical  rate 
constants  for  these  processes  are  10“ cm’s"', 
IQ-^’/r^^  cm'"  s"',  and  cm’s”',  respectively, 

where  the  electron  temperature  T,  has  units  of  eV.  For  typi¬ 
cal  bulk  electron  temperatures  of  1-2  eV  and  ion  densities  of 
lO’-lO'”  cm"',  only  dissociative  recombination  is  a  signifi¬ 
cant  electron/ion  loss  mechanism  when  compared  to  diflTu- 
sion  to  the  walls,  attachment,  and  ion-ion  neutralization. 

The  products  for  dissociative  recombination  of  silane 
ions,  as  well  as  the  rate  constants,  have  not  been  measured. 
The  rate  of  dissociative  recombination  of  CH/  ,  though,  has 
been  obtained  by  Kline  ”  By  analogy,  the  cross  section  cor- 
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responding  to  this  rate  was  used  for  the  dissociative  recom¬ 
bination  of  all  silane  ions  in  the  model.  These  reactions  have 
the  general  form 

-t-H.  (11) 

The  actual  product  distribution  depends  upon  the  details  of 
the  curve  crossing  of  the  potential  surfaces  for  the  intermedi¬ 
ate  excited  state  and  products.  There  would  though,  be  suffi¬ 
cient  energy  in  the  intermediate  state  to  break  any  multiple 
bond  and  therefore  result  in  a  variety  of  products.  In  the 
absence  of  knowledge  of  the  actual  product  distribution,  we 
chose  the  indicated  products  for  all  dissociative  recombina¬ 
tions. 


C.  Gas  heating 

Heating  of  the  gas  by  power  deposition  is  important  in  at 
least  two  respects.  First,  since  the  deposition  reactor  is  iso- 
baric,  the  reactant  gas  density  is  inversely  proportional  to 
the  local  gas  temperature.  Second,  for  sufficiently  high  tem¬ 
peratures  ( r,  >800  K)  gas  phase  pyrolysis  will  occur. 

Heat  transfer  within  the  plasma  can  be  described  by 
Picks  law  provided  the  molecular  mean  free  path  is  much 
less  than  the  electrode  separation  (see  Sec.  II  A);  and  AV 
a  4  At,  where  a  is  the  thermal  diifusivity.  A  is  the  reactor 
diffusion  length,  and  At  is  the  gas  residence  time  between  the 
electrodes.  The  thermal  diffusivity  is  a  =  v/(Prp),  where  v 
is  the  viscosity,  Pr  is  the  Prandtl  number,  and  p  is  the  gas 
mass  density  (see  Sec.  II  D).  For  our  conditions  in  pure 
silane,  asrlx  10^  cm^  s"‘.  Thermal  diffusion  is  therefore  a 
good  description  for  heat  transfer  provided  that  A/>  0.05  </  * 
ms,  where  the  electrode  spacing  d  is  in  cm.  This  condition  is 
met  for  all  electrode  spacings  and  residence  times  of  interest. 
The  gas  temperature  is  therefore  in  quasiequilibrium  with 
discharge  power  deposition  and  the  gas  temperature  in  a 
symmetric  parallel  plate  geometry  is  approximately 

T,  is  the  substrate  temperature,  k  is  the  thermal  conductivity 
( sSx  10“’ W  cm“‘  K“'  for  pure  silane),  Pis  power  depo¬ 
sition  excluding  that  dissipated  by  ions  striking  the  elec¬ 
trodes,  and  X  is  measured  from  the  midplane  of  the  reactor. 
In  this  expression,  we  have  assumed  that  power  deposition 
by  electron  collisions  is  uniform  throughout  the  discharge. 
The  average  gas  temperature  is  T  =  T,  -t-  {2d^P)/{it^K). 
For  example,  with  an  electrode  separation  of  2.5  cm  and 
power  deposition  of  25  mW  cm“’,  the  centerline  gas  tem¬ 
perature  T,  -f-  200  K  and  the  average  gas  temperature  is  T, 
+  130  K. 

D.  Transport  coefficients  and  thermodynamic 
properties 

Diffusion  coefficients  for  neutral  species  were  obtained 
by  assuming  diffusion  through  the  feedstock  gases.  Coeffi¬ 
cients  were  computed  using  Lennard-Jones  parameters  in 
the  manner  described  in  Ref.  80.  Lennard-Jones  parameters 
for  silane  species  were  assumed  to  be  the  same  as  their  meth¬ 
ane  counterparts,  when  available,  and  estimated  when  not. 


The  parameters  as  used  in  the  model  are  listed  in  Table  VI. 
These  parameters  were  also  estimated,  with  qualitatively  dif¬ 
ferent  results  in  Ref.  46. 

The  heat  capacities  of  the  input  gases  ( As,  H2,  SiH^,  and 
Si2H«)  were  obtained  as  a  function  of  temperature  from  a 
various  sources  during  model  development.  They  are  now 
summarized  in  a  thermodynamic  data  base  published  as  Ref. 
81.  In  that  compilation,  the  heat  capacity  is  expressed  as  the 
sum 

c,/R=  X 

i  -  OS 

The  coefficients  for  silane  (/  =  0,5)  are:  2.475  1668, 
9.003  721  18X10“’,  2.185  394  35x10“®,  -2.681423  08 
X10“’,  -6.621  080  64X10“”. 

The  thermal  conductivity,  k,  of  the  gas  mixture  was  ob¬ 
tained  from  K  =  wc^/Pr*’  where  v  is  the  viscosity  and  Pr  is 
the  Prandtl  number.  The  viscosity  of  the  mixture  was  com¬ 
puted  from  Lennard-Jones  parameters  as  described  in  Ref. 
80.  The  Prandtl  number  is  Pr  =  4y/(9j'  —  5 )  where  the  spe¬ 
cific  heat  ratio  y  =  (c^/R)/{c^/R  —  1 ). 

The  mobility  of  ions  in  the  gas  mixture  were  calculated 
from  p  =  1/1,  ( 1//X, ),  where  p,  is  the  mobility  of  the  ion  in 
gas  /.  Mobilities  of  silane  ions,  Si.H,;!',  in  silane  were  ob¬ 
tained  from  Chatham  and  Gallagher,  '*  having  low-field  lim¬ 
its  of  1 180  and  675  cm’  Torr/V  s  forn  =  1  and  n  =  2.  The 
mobilities  of  At*  in  At,  H*  in  H2,  and  H,"^  in  H2  were 
obtained  from  the  compilations  of  Ellis  and  co-workers.®’ 
The  mobilities  of  other  ions  were  calculated  from  po 
=  2.46X  10*/[  (M /m°*) T°*a]  cm’/V  s,  where po  is  the 
mobility  at  standard  temperature  and  pressure,  M  is  the  mo¬ 
lecular  weight  (amu),  m  is  the  reduced  mass  with  the  colli¬ 
sion  partner  (amu),  Tg  is  the  gas  temperature  (K),  and  cr  is 
the  sum  of  the  elastic  and  charge  exchange  cross  sections 
(A’).  The  charge  exchange  cross  sections  were  obtained 
from  the  charge  exchange  rate  constants  listed  in  Table  IV. 
The  elastic  cross  section  was  assumed  to  be  10  A’  for  all 
species. 

E.  Boundary  and  Initial  conditions,  and  deposition  rate 

The  simulation  begins  with  the  density  of  all  radicals 
being  zero,  and  the  density  of  the  feedstocks  being  deter¬ 
mined  by  the  absolute  pressure,  gas  temperature,  and  their 


TABLE  VI.  Estimated  Lennard-Jones  parameters  for  neutral  species. 


Species 

<r(A) 

e(K) 

Ar 

3.44 

120 

H. 

2.94 

36 

H 

2.50 

30 

Si 

3.84 

100 

SiH, 

3.84 

140 

Si.H, 

4.22 

185 

Si.H, 

4.23 

195 

Si.H, 

4  23 

205 

Si.H, 

4.33 

217 

Si.H„ 

4.42 

230 

Si.H, 

5.08 

254 

Si,H. 

5.80 

254 

SUH, 

6.50 

254 
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mole  fractions.  The  electron  and  ion  densities  are  set  to  small 
initial  values  ( =:  10^  cm  ~  ^ .  All  positive  ions  are  assumed  to 
be  neutralized  at  the  boundaries  by  either  electrons  or  nega¬ 
tive  ions,  with  the  products  returning  to  the  plasma.  The 
fraction  of  negative  ions  which  are  neutralized  at  the  wall  is 
determined  by  the  appropriate  ratio  of  fluxes  between  elec¬ 
trons  and  negative  ions. 

Molecules  which  do  not  stick  to  the  surface  have  the 
boundary  condition  dN/dx{x=  ±<//2)=0.  Radicals 
with  unity  sticking  coefficient  have  the  boundary  condition 
ofiV(x  =  ±  d  /2 )  =  0.  Those  radicals  with  a  nonunity  stick¬ 
ing  coefficient  s  are  treated  as  having  a  net  flux  “sink”  at  the 
wall  of  magnitude  The  sticking  coefficients  of 

saturated  molecules  is  zero,  whereas  those  for  Si„  H,„  ,m<ln 
are  unity.  The  sticking  coefficient  of  radicals  with 
m  =  2n  -h  1  (e.g.,  SiH])  is  O.IS,  in  accordance  with  recent 
measurements.^^  The  initial  transient  in  gas  temperature, 
which  lasts  for  a  ms  or  less,  is  usually  ignored  and  the  initial 
gas  temperature  is  set  equal  to  the  values  as  given  in  Sec. 
II  C. 

The  deposition  rates  of  amorphous  silicon  on  the  sub¬ 
strates  at  ac  =  ±d/2  are  given  by 

R,=2,b\s,D,VN,^x= 

where  b  is  the  average  bond  length  (2.8  A)  between  silicon 
atoms  in  the  fllm,  n,  is  the  number  of  silicon  atoms  in  radical 
/,  s,  is  the  sticking  coefficient,  Dj  is  the  diffusion  coefficient, 
and  number  density  of  the  radical  is  Af,  . 

III.  RESULTS 

In  this  section,  results  from  the  plasma  deposition  mod¬ 
el  will  be  discussed  for  discharges  in  mixtures  of  silane  and 
argon  or  hydrogen.  In  each  case,  we  have  assumed  that  the 
parallel  plate  reactor  is  symmetric;  therefore  results  are  plot¬ 
ted  as  a  function  of  position  only  from  the  left  electrode  to 
the  center  line.  Power  deposition  is  cited  as  W  cm~’  aver¬ 
aged  over  the  volume  of  the  reactor.  The  substrate  tempera¬ 
ture  is  5(X)  K.  Species  densities  presented  as  a  function  of 
position  between  the  electrodes  are  for  a  time  after  the  input 
species  have  flowed  downstream  and  reached  their  quasi- 
steady-state  conditions.  These  conditions  are  defined  as 
when  mono-  and  di-silane  radicals  and  ions  have  reached  an 
equilibrium,  and  the  silane  has  been  depleted  by  <0.1.  This 
value  corresponds  roughly  to  20  mJ/Torrcm“’  of  silane. 
Higher  silane  radicals,  saturated  molecules,  and  DUST  may 
continue  to  increase  in  density  beyond  the  equilibration  time 
for  lower  radicals. 

A.  Typical  results 

Species  densities  for  a  discharge  sustained  in  250  mTorr 
of  silane  between  electrodes  with  2.S-cm  separation  with  a 
power  deposition  of  25  mW  cm“’  are  shown  in  Figs.  5-8. 
The  computed  deposition  rate  of  <i-Si:H  for  these  conditions 
is  s  250  A/min,  which  agrees  well  with  experimental  values 
for  similar  conditions.^  The  time  evolution  of  the  density  of 
saturated  molecules  is  shown  in  Fig.  5(a).  At  the  end  point, 
silane  has  been  depleted  by  8%.  The  stable  product  with  the 


highest  density  is  Hj  representing  a  net  deposition  reaction 
of  SiH4(g)— Si:H(5)  Hjfg),  where  (g)  and  (s)  denote 
gas  phase  and  solid  species,  respectively.  Based  only  on  volu¬ 
metric  processes,  the  next  most  abundant  saturated  product 
is  initially  SijHg,  but  after  a  short  flow  downstream  is  re¬ 
placed  by  SijH^.  This  condition  is  likely  a  result  of  the  rapid 
reaction  sequence 

SiH. 

SiH  j  -t-  SiH«  -  Si  jH?*  -  Si,H4  -I-  Hi  -  Si,H»  -f-  H^  (13a) 
L-t-M-SiiHft-SijHg  (13b) 

SiH, 

which  depletes  disilane.  Only  after  a  sufficient  density  of  Hj 
builds  up  does  the  reverse  reaction  Hi  -I-  SiiH4— SiiH^  be¬ 
come  frequent  enough  to  reduce  the  relative  importance  of 
reaction  (13a).  Operating  at  a  higher  pressure  also  reduces 
the  relative  importance  of  this  reaction,  thereby  increasing 
the  density  of  disilane  by  increasing  the  branching  ratio  of 
reaction  ( 13b).  This  effect  that  can  also  be  obtained  by  dilut¬ 
ing  with  Hi.  In  fact,  in  a  SiH4/Hi  =  1/1  gas  mixture  for  the 
same  discharge  conditions,  the  density  of  disilane  is  always 
at  l«ist  twice  that  for  trisilane.  Higher  silanes  having  more 


FIG.  5.  Volume  averaged  densities  of  saturated  molecules  in  an  rf  discharge 
sustained  in  silane  as  a  function  of  Row  time;  (a)  without  the  wall  associ- 
ation  reaction  2-SiH,(g)  —2'SiH,(a)  —SijHttg),  (b)  with  the  wall  associ¬ 
ation  reaction.  The  gas  pressure  is  230  mTorr  with  power  deposition  of  23 
m  W  cm  ~ '  and  substrate  temperature  of  300  K.  The  electrode  separation  is 
2.3  cm.  Hydrogen  is  the  most  abundant  saturated  product,  but  only  with  the 
wall  association  reaction  is  disilane  the  most  abundant  higher  silane. 
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(a) 


FIG.  6.  Volume  avenged  n> 
dical  densities  for  the  condi¬ 
tions  of  Fig.  S;  (a)  the  sums 
of  monomer,  dimer,  and 
trimer  radical  densities,  (b) 
the  individual  monomer  n- 
dical  deiisities.The  monomer 
radicals  are  the  most  abun¬ 
dant,  however  only  the  dimer 
radicals  come  quickly  into 
equilibrium. 


than  four  silicon  atoms  obtain  only  modest  densities  for 
these  conditions.  The  density  of  the  higher  silanes,  though, 
depends  critically  on  the  choice  of  the  chain  propagating  rate 
constant,  which  is  1.0 X  10“  "  cm“^  s“'  in  this  example. 

It  appears  that  at  tow  pressures  (<200  mTorr)  gas 
phase  reactions  do  not  favor  formation  of  disilane.  Large 
mole  fractions  of  disilane,  though,  have  been  observed  for 
similar  discharge  conditions  by  Chatham  and  Gallagher.'^ 
They  attributed  the  high  density  to  the  heterogeneous  wall 
association  reaction 


FIG.  7.  The  density  of  monomer  ((•)  top)  end  dimer  radicals  ((b)  bot¬ 
tom  ]  as  a  function  of  position  between  the  electrodes  for  the  discharge  con¬ 
ditions  of  Fig.  S.  The  discharge  is  symmetric.  The  electrode  is  at  the  left 
boundary,  and  the  discharge  centerline  is  at  the  right.  SiH,  and  SiH  are 
peaked  near  the  electrode  as  electron  impact  dissociation  of  silane  are  their 
dominant  source  and  they  rapidly  react  with  silane.  SiHj  is  nearly  uniform 
as  a  Ametion  of  position  due  to  its  mote  distributed  production  and  low 
sticking  coefficient  SijH,  is  peaked  near  the  electrode  due  to  its  source  be¬ 
ing  insertion  of  SiH,  into  silane.  [Note;  The  label  on  the  lowest  curve  in 
7(b)  is  incorrect.  It  should  be  SijHj.) 


2SiH,(g)  -2SiH,(fl)  -SUHJg),  ( 14) 

where  (a)  denotes  adsorbed  species.  This  reaction  was  in¬ 
cluded  in  the  model  for  the  nonsticking  silyl  flux  striking  the 
electrodes.  The  density  of  SiiH^  increased  as  shown  in  Fig. 
S(b).  For  these  results,  half  of  the  nonsticking  silyl  flux 
striking  the  electrodes  participates  in  reaction  14.  The  wall 
association  reaction  results  in  disilane  becoming  the  most 
abundant  higher  silane  molecule,  consistent  with  Chatham 
and  Gallagher’s  observations. 

The  time  evolution  of  SiH„, ,  Si^H,, ,  and  SijH,,,  radicals 
are  shown  in  Fig.  6(a) .  The  dimer  radicals  quickly  come  into 
equilibrium  and  the  trimer  radicals  do  so  at  a  later  time.  The 
monomer  radicals  display  an  extremum,  due  primarily  to 
the  density  of  SiH3.  This  behavior  is  shown  in  more  detail  in 
Fig.  6(b),  where  the  densities  of  the  individual  SiH,,  radi¬ 
cals  are  plotted.  The  density  of  SiH  quickly  comes  into  equi¬ 
librium  between  its  production  by  electron  impact  dissocia¬ 
tion  of  silane  and  its  subsequent  insertion  into  silane. 
Although  the  production  rate  of  SiHj  is  relatively  constant, 
the  density  of  SiHj  initially  peaks  then  settles  to  a  near  equi¬ 
librium  value.  (The  production  of  SiH,  by  hydrogen  ab¬ 
straction  from  SiH^  exceeds  that  by  electron  impact  disso¬ 


ciation  of  silane  in  the  ratio  of  1.5/1.)  The  decrease  in 
SiH}  density  is  due  to  the  increase  in  density  of  other  radicals 
and  saturated  molecules  with  which  SiHj  reacts.  These  spe¬ 
cies  are  principally  H,  SiHj,  SijH,,  and  higher  silanes.  The 


FIG.  8.  Silane  ion  denaitiei  aa  a  function  of  position  between  the  electrodes 
for  the  conditions  of  Fig.  7.  The  progressively  more  uniform  densities  going 
from  SiH,'  toSi,H,'  indicate  that  the  source  of  the  ion  is  further  displaced 
from  the  originating  step  in  the  chain,  electron  impact  ionization  of  silane. 
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density  of  SiH2  increases  monotonically  towards  an  equilib¬ 
rium  value  as  its  dominant  destruction  mechanism,  insertion 
into  silane,  does  not  appreciably  change  its  rate. 

The  spatially  dependent  densities  of  SiH„  radicals  at  the 
terminal  time  shown  in  the  previous  figures  are  in  Fig.  7(a). 
The  spatial  distribution  of  SiHj  is  fairly  uniform  whereas 
those  for  SiHj  and  SiH  are  maximum  near  the  electrode  at 
the  location  of  the  peak  in  the  production  rate  by  electron 
impact  dissociation.  The  density  of  SiH3  near  the  electrode  is 
about  two-thirds  that  of  the  maximum,  a  consequence  of  its 
nonunity  sticking  coefficient  (0.15);  the  densities  of  SiH2 
and  SiH  are  negligible  at  the  electrode  due  to  their  unity 
sticking  coefficient.  The  spatially  dependent  densities  of  the 
dimer  radicals  are  shown  in  Fig.  7(b).  The  density  of  SijHj 
is  largest  near  the  electrode  because  neutralization  of  Si2H2'^ 
at  the  surface  is  an  important  source  for  the  radical.  The 
density  of  Si2H4  has  an  extremum  near  the  electrode  at  the 
site  of  production  of  SiH2,  as  the  insertion  of  that  radical  into 
silane  is  the  major  source  of  Si2H4.  The  densities  of  Si2H]  and 
Si2H,  have  the  characteristic  profiles  of  radicals  having  volu¬ 
metric  sources  and  edge  sinks.  Initially,  the  major  source  of 
SijH]  is  insertion  of  SiH  into  silane;  downstream  after  the 
density  of  disilane  increases,  the  major  source  is  hydrogen 
abstraction  of  disilane. 

Ihe  spatially  dependent  densities  of  Sinj**,  SiHj^, 
SijHT** ,  and  Si^H^^  are  shown  in  Fig.  8  for  the  same  condi¬ 
tions.  This  sequence  of  ions  represents  species  which  result 
from  1, 1  or  2, 2,  and  three  collisions,  respectively,  where  the 
first  collision  is  electron  impact  ionization  of  silane.  The  den¬ 
sity  of  SiHj*  closely  mirrors  the  source  function.  The  SiH,* 
ion  is  produced  either  directly  by  electron  impact  dissocia¬ 
tion  or  by  the  ion-molecule  reaction  of  SiH2*  with  silane. 
The  Si2H7*  ion  results  from  reaction  of  SiH,*  with  silane; 
and  the  trimer  ion  results  from  collisions  of  dimer  ions  with 
silane.  As  more  collisions  are  required  to  form  the  species  by 
ion-molecule  reactions,  the  ion  appears  to  have  a  more  uni¬ 
form  rate  of  production.  This  issue  is  further  discussed  in  the 
following  section. 

B.  Rate  coefficient  versus  production  rate  versus 
survival  rate 

The  local  electron  impact  rate  coefficient  A;  is  a  function 
of  phase  during  the  rf  cycle  and  position  between  the  elec¬ 
trodes.  The  rate  of  production  R  of  a  radical  by  electron 
impact,  though,  is  the  product  R  =  n,kN,  where  N  is  the 
heavy  particle  reactant  density.  The  rate  constant  is  maxi¬ 
mum  near  the  electrodes,  whereas  the  electron  density  is 
maximum  near  the  center  of  the  discharge  due  to  the  loss  by 
diffusion  to  the  walls.  The  heavy  particle  density  is  a  mini¬ 
mum  at  the  center  line  because  of  gas  heating.  The  result  is 
that  the  production  rate  has  a  maximum  displaced  from  the 
electrodes.  The  displacement  of  the  rate  of  production  away 
from  the  electrodes  has  important  implications  with  respect 
to  the  flux  of  radicals  incident  onto  the  substrates. 

Film  properties  are  a  function  of  the  composition  of  the 
flux  of  radicals  incident  on,  and  sticking  to,  the  substrate. 
The  composition  of  the  flux  which  strikes  the  substrate  is 
different  from  that  which  results  from  electron  impact  disso¬ 


ciation  of  silane.  This  discrepancy  results  from  reactions 
which  occur  between  the  site  of  production  of  radicals  by 
electron  impact  and  their  arrival  at  the  substrate.  In  this 
respect,  we  define  the  survival  ratio  as 

^=5Z)/A'Ax,A^„  (15) 

/  I 

where  D  is  the  diffusion  coefficient,  s  is  the  sticking  coeffi¬ 
cient.  A  is  the  diffusion  length  and  x,  is  the  rate  of  reaction 
with  heavy  particle  i  having  density  IV,  .  If  1,  the  rate  of 
diffusion  is  large  compared  to  the  rate  of  chemical  reaction 
and  the  radical  survives  to  reach  the  surface  as  produced  by 
electron  impact  dissociation.  The  spatially  dependent  den¬ 
sity  of  a  sticking  radical  then  resembles  that  having  a  uni¬ 
form  volume  source  and  edge  sink.  If^<  1 ,  the  radical  reacts 
soon  after  production  and  its  spatially  dependent  density 
mirrors  that  of  the  production  rate.  A  comparison  of  the 
spatial  distribution  of  production  and  density  of  a  radical  is 
therefore  an  indication  of  its  survivability  and  contribution 
to  growth  of  the  film. 

To  illustrate  these  points,  we  examine  the  spatial  distri¬ 
butions  of  SiHj  and  SiH,  as  a  function  of  silane  mole  fraction 
in  hydrogen.  These  radicals  are  the  major  products  of  elec¬ 
tron  impact  dissociation  of  silane.  SiH,  is  not  particularly 
reactive  with  silane  (/r=;  10~  cm^  s~ ' )  and  has  a  low  stick¬ 
ing  coefficient,  while  SiH,  is  reactive  with  silane  (kx  10~ 
cm^  s~ ' )  and  has  a  unity  sticking  coefficient.  The  densities 
of  SiH,  and  SiH,  as  a  function  of  position  and  dilution  of 
silane  by  hydrogen  are  shown  in  Fig.  9.  The  distributions  of 
SiH,  do  not  appreciably  change  as  silane  is  diluted  with  hy¬ 
drogen.  The  distribution  for  the  SiHyH,  =  1/1  mixture, 
though,  is  displaced  from  the  electrode.  This  motion  reflects 
the  fact  that  the  production  of  SiH,  is  more  dominantly  by 
hydrogen  abstraction  from  silane  in  the  1/1  mixture  as  op¬ 
posed  to  electron  impact  dissociation  of  silane  (a  ratio  of 
10/1  compared  to  =;  1.5/1  fora  pure  silane).  The  produc¬ 
tion  rate  by  electron  impact  dissociation  is  maximum  near 
the  electrodes,  whereas  that  by  hydrogen  abstraction  is  max¬ 
imum  in  the  bulk,  as  shown  in  the  figure.  The  distributions  of 
SiH,  are  highly  peaked  near  the  electrode  in  a  pure  silane  and 
maximum  in  the  bulk  for  a  SiH^H,  «=  1/1  mixture.  This 
trend  indicates  production  occurs  near  the  electrodes  and 
the  primary  reactant  is  silane.  The  survival  ratio  for  SiH, 
analogously  increases  with  dilution  with  H,,  as  shown  in 
Fig.  9(c). 

Similar  trends  will  occur  if  the  location  of  production  of 
radicals  is  displaced  from  the  electrodes.  This  is  best  illus¬ 
trated  in  a  negative  glow  dc  discharge  where  the  production 
rate  of  radicals  by  electron  impact  is  maximum  at  the  Klge  of 
the  cathode  fall.  The  cathode  fall  thickness  is  a  function  of 
many  parameters,  but  generally  increases  with  increasing 
cathode  fall  voltage.  Therefore,  operating  at  larger  cathode 
fall  voltages  moves  the  production  of  silane  radicals  by  elec¬ 
tron  impact  away  from  the  cathode,  thereby  changing  the 
probability  that  a  particular  radical  survives  to  reach  the 
substrate  (cathode).  We  modified  the  plasma  chemistry 
model  to  simulate  a  dc  negative  glow  discharge  and  comput¬ 
ed  the  ratio  of  the  fluxes  of  SiH,  to  SiH,,  ^(SiHjl/^lSiH,), 
incident  on  the  cathode.  Typical  SiH,  and  SiH,  densities  are 
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FIG.  9.  Radical  densities  for  rf  discharges  in  mixtures  of  silane  and  hydro¬ 
gen;  (a)  SiH,  and  (b)  SiHj.  The  densities  of  the  radicals  for  each  mixture 
are  normalized  to  their  maximum  value  to  emphasize  the  difference  in  their 
spatial  distributions.  The  density  of  SiH,  is  always  fairly  uniform,  reflecting 
its  low  reaction  rate  and  low  sticking  coefficient.  The  density  of  SiH.  is 
uniform  only  with  low  mole  fractions  of  silane.  The  relative  probability  of 
survival  of  SiH,  radicals  to  reach  the  substrate  increases  with  decreasing 
silane  mole  fraction,  as  shown  in  (c). 


shown  in  Fig.  10(a)  for  a  cathode  fall  length  of  2.0  cm.  The 
survival  probability  for  SiH^  to  reach  the  cathode  is  low  due 
to  the  long  residence  time  of  SiH,  radicals  to  diffuse  to  the 
cathode,  and  the  high  probability  of  reaction  before  reaching 
it.  As  a  result,  the  flux  ratio  ^(SiH])/^(SiH.)  is  large, 
>  1(X).  As  the  cathode  fall  distance  decreases,  this  ratio  in¬ 
creases,  a  result  of  the  increased  probability  of  survival  for 
SiH;  [see  Fig.  10(b)  ]  as  its  residence  time  decreases.  Since 
the  quality  of  a-Si:H  films  depends  on  this  ratio  (see  Sec. 
Ill  D),  the  thickness  of  the  cathode  fall  directly  effects  the 
properties  of  the  film. 


(b) 

FIG.  10.  Radical  densities  and  flux  ratios  fora  dc  discharge  in  ISOmTorrof 
silane  sustained  between  two  parallel  electrodes  (cathode  at  left,  anode  at 
right ) :  ( a )  SiH ,  and  SiH,  density  for  a  cathode  fall  distance  of  2  cm  with  the 
distribution  of  production  of  radicals  by  electron  impact  dissociation.  (b> 
The  ratio  of  the  flux  of  SiH,  to  SiH,  incident  on  the  cathode  as  a  function  of 
cathode  fall  distance.  The  ratio  increases  with  increasing  cathode  fall  dis¬ 
tance  due  to  the  increased  residence  time  for  SiH,. 


C.  Electrode  separation 

A  common  measure  of  power  deposition  in  parallel 
plate  discharges  is  W  cm~^  based  on  electrode  area.  This 
scaling  parameter  implies  that  the  deposition  rate  of  a-Si:H 
should  be  independent  of  electrode  separation.  This  depo¬ 
sition  rate,  though,  has  generally  been  observed  to  decrease 
with  increasing  electrode  separation.^'^  One  would  predict 
that  if  specific  power  deposition  ( W  cm~^)  is  held  constant, 
the  deposition  rate  should  increase  with  increasing  electrode 
separation  since  the  area  to  volume  ratio  scales  as  d.  This 
trend  is  also  not  always  observed.  There  must  therefore  be 
other  processes  which  effect  deposition  rate  when  the  elec¬ 
trode  separation  is  changed.  In  this  section,  we  will  discuss 
some  of  these  processes. 

Changes  in  the  geometry  of  the  reactor  (that  is,  elec¬ 
trode  separation )  can  affect  the  manner  of  power  deposition, 
residence  time  for  radicals  and  ions,  and,  as  a  result,  the 
composition  of  the  flux  incident  onto  the  substrate.  Power 
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deposition  by  electrons  occurs  in  two  distinct  regions;  by 
higher  energy  electrons  in  and  near  the  sheaths  and  by  lower 
energy  electrons  in  the  bulk  plasma  which  resembles  a  posi¬ 
tive  column.  The  power  deposition  associated  with  the 
sheaths  can  more  easily  access  processes  with  higher  thresh¬ 
old  energies.  The  power  deposition  in  the  bulk  typically  ac¬ 
cesses  lower  threshold  processes.  Since  the  sheaths  occupy  a 
fixed  volume  near  the  electrodes,  the  fraction  of  power  de¬ 
posited  by  processes  associated  with  the  sheaths  decreases  as 
the  electrode  separation  increases,  and  that  deposited  in  the 
bulk  by  lower  energy  electrons  fractionally  increases.  The 
fractional  amount  of  power  dissipated  by  ions  striking  the 
electrodes  also  changes,  usually  decreasing  as  the  electrode 
separation  increases.  The  contribution  for  our  conditions 
therefore  remains  small,  <  10%. 

Consequences  of  this  change  in  the  manner  of  power 
deposition  are  shown  in  Fig.  1 1  where  the  densities  of  SiHj 
radicals  and  vibrationally  excited  silane  ((SiH4(<;i,3)) 
+  [SiH4(u2,4)  ])  are  plotted  as  function  of  position  in  the 
discharge  for  three  different  electrode  separations.  The  dis¬ 
charge  conditions  are  =  1/1,  250  mTorr,  with 

power  deposition  constant  at  25  m W  cm  ~  As  the  electrode 


b) 

FIG.  II.  The  density  of  SiH,  ((■),  lop)  and  vibrationally  excited  silane. 
SiHt(v),  I  (b),  bottom)  as  a  function  of  position  between  the  electrodes  in 
an  rf  discharge  for  dilTerent  electrode  separations.  The  discharge  conditions 
are  otherwise  the  same  as  for  Fig.  S.  Relative  rates  of  high  threshold  pro¬ 
cesses  such  as  dissociation  decrease  with  increasing  electrode  separation  ax 
the  discharge  becomes  dominated  by  the  positive  column. 


separation  t/  increases,  the  density  of  SiHj  decreases  while 
the  density  of  SiH4(t;)  increases.  This  trend  indicates  that  as 
d  increases,  a  larger  fraction  of  power  deposition  occurs 
through  lower  energy  processes,  such  as  vibrational  excita¬ 
tion,  resulting  in  less  dissociation.  The  discharge  appears 
more  like  a  positive  column.  Less  dissociation  leads  to  a 
smaller  radical  density  and  hence  lower  deposition  rate  of  a- 
Si:H.  The  predicted  deposition  rate  of  n-Si;H  as  a  function  of 
electrode  spacing  appears  in  Fig.  12  and  decreases  with  in¬ 
creasing  d.  For  comparison,  experimental  results  from  Ross 
and  Jaklik^  and  Ishihara  et  al.^  are  also  shown. 

A  second  consequence  of  increasing  d  is  that  the  average 
residence  time  for  radicals  also  increases,  thereby  increasing 
the  probability  of  chemical  reactions  after  the  initial  electron 
impact  dissociation  of  silane.  The  change  in  residence  will 
result  in  a  change  in  the  composition  of  the  flux  incident  onto 
the  electrodes.  An  example  is  shown  in  Fig.  13  where  the 
ratio  of  fluxes  incident  on  the  electrodes  for  SiHy  and  SiHj, 
^(SiHjl/^fSiHj),  is  plotted  as  a  function  of  d.  This  ratio 
increases  with  increasing  electrode  separation  as  the  prob¬ 
ability  for  reaction  of,  and  consumption  of,  SiH;  by  insertion 
into  silane  increases  with  increasing  residence  time.  The  de¬ 
pletion  of  SiH)  increases  at  a  lower  rate.  Also  shown  in 
Ftg.  1 3  is  the  density  of  DUST,  which  decreases  with  increas¬ 
ing  separation.  Both  of  these  trends  coincide  with  film  quali¬ 
ty  increasing  with  large  electrode  separation,  trading  off 
against  lower  deposition  rates.  Further  parametric  tradeoffs 
for  film  quality  are  discussed  in  the  following  section. 

D.  Contributions  to  deposition  ' 

Conditions  for  which  the  a-Si:H  film  grown  from  silane 
discharges  is  judged  “device  quality”  are  empirically  known. 
Generally,  one  operates  where  gas  phase  polymerization  and 
particulate  formation  is  not  likely  to  occur  (i.e.,  low  power 
and  low  partial  pressures  of  silane),  and  conditions  where 
the  ion  flux  is  a  small  fraction  of  the  radical  flux.  The  latter 
criteria  insure  that  the  process  is  dominantly  chemical  vapor 
deposition,  as  opposed  to  physical  vapor  deposition*^  and 
minimizes  the  formation  of  columnar  structures.**  High- 
quality  films  have  also  been  correlated  with  conditions 


FIO.  1 2.  Computed  depotition  rate  of  a-Si:H  at  a  funclioa  of  electrode  sepa¬ 
ration  compared  to  the  results  of  Ross  and  Jaklik  (Ref.  3)  and  Ishihara ero/. 
(Ref  6)  for  similar  conditions.  The  results  have  been  normalized  to  Ross 
and  Jaklik  at  d  •  2  cm. 
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FIG.  13.  The  ratio  of  the  fluxofSiH]  toSiH.  incident  on  the  substrate  and 
the  density  of  DUST  as  a  function  of  electrode  separation.  The  increase  in 
this  ratio,  and  decrease  in  the  density  of  DUST,  with  increasing  electrode 
separation,  are  indications  of  increasing  quality  of  the  a-Si:H  films  deposit¬ 
ed. 

where  the  ratio  of  SiHj  to  SiH^  incident  on  the  surface  is 
high.  SiH,  has  a  high  surface  mobility  compared  to  SiH;."** 
Films  grown  from  primarily  SiH3  radicals  therefore  will  re¬ 
semble  the  preferred  chemical  vapor  deposition,  as  opposed 
to  physical  vapor  deposition.  In  this  respect,  we  quantify 
three  measures  of  film  quality;  a  large  fraction  of  deposition 
resulting  from  lower  silanes  (Si,H„,  /t<2),  a  large  ratio  of 
deposition  from  SiH,  as  compared  to  SiH^,  and  large  ratio  of 
the  flux  of  radicals  containing  Si  atoms  compared  to  the  flux 
of  ions  of  any  type  incident  onto  the  substrate. 

Discharges  were  simulated  for  gas  mixtures  of  silane 
diluted  in  argon  and  hydrogen.  The  argon  diluted  gas  mix¬ 
tures  are  at  a  total  pressure  ofO.S  Torr  and  power  deposition 
of  50  m  W  cm  the  hydrogen  diluted  mixtures  are  at  a  total 
pressure  of  0.25  Torr  and  power  deposition  of  25  mW  cm  ~  \ 
each  case  being  the  same  specific  power  loading.  The  results 
are  shown  in  Fig.  14  as  a  function  of  the  ratio  of  [SiH^J/ 
([SiH^l  -t-  [Aror  H;l). 

For  both  dilutions,  the  fraction  of  deposition  due  to  low¬ 
er  silanes  decreases  with  increasing  silane  mole  fraction  as 
the  probability  of  radicals  colliding  with  other  silane  species 
before  being  incorporated  into  the  film  increases  with  in¬ 
creasing  silane  mole  fraction.  The  fraction  of  deposition 
from  lower  silanes  decreases  from  =:0.9  to  ~0.6  (see  Fig. 
14),  and  by  this  measure,  film  quality  should  decrease  with 
increasing  silane  mole  fraction.  By  other  measures,  though, 
film  quality  should  increase  with  increasing  silane  mole  frac¬ 
tion.  For  example,  the  ratio  of  SiHj/SiH,  incorporated  into 
the  film  increase  from  about  2  to  20  as  the  mole  fraction  of 
silane  in  hydrogen  increases,  as  shown  in  Fig.  14(a).  This 
trend  is  a  consequence  of  the  higher  probability  of  SiH,  in¬ 
serting  into  silane  before  the  radical  reaches  the  surface. 

Our  other  measure  of  film  quality,  the  ratio  of  the  silicon 
radical  flux  to  the  total  ion  flux,  also  increases  with  increas¬ 
ing  silane  mole  fraction,  as  shown  in  Fig.  14(b).  This  ratio 
increases  by  a  factor  of  4  from  a  5%  mixture  of  silane  in 
argon  to  pure  silane.  The  increase  is  a  result  of  a  larger  frac¬ 
tion  of  discharge  power  being  coupled  through  the  neutral 
dissociative  modes  of  silane,  and  a  decrease  in  the  average 
electron  energy  with  increasing  silane  fraction,  which  favors 
neutral  dissociation  over  dissociative  ionization. 


[SiH4] 


[SiH4],[H2] 

(0) 


(b) 

FIG.  14.  Measures  of  film  quality  for  discharges  sustained  in  mixtures  of 
silane  diluted  by  hydrogen  ((a),  top)  and  silane  dilu'-'d  by  argon  [(b), 
bottom).  Quality  increases  with  an  increasing  fraction  of  the  deposition 
from  SiH)  compared  to  SiH-,  increasing  Si  radical  flux  compared  to  ion 
flux,  and  increasing  deposition  by  lower  silanes.  These  measures  of  quality 
have  different  dependencies  on  silane  mole  fraction. 


These  measures  of  film  quality  have  different  dependen¬ 
cies  on  silane  mole  fraction  and  partial  pressure.  A  tradeoff 
must  be  performed  to  optimize  the  quality  of  a-Si:H  films 
when  diluting  silane  or  changing  its  panial  pressure. 

Other  parameters  are  also  important  to  film  deposition, 
particularly  the  rate  of  deposition.  The  deposition  rate  does 
not  simply  scale  with  silane  mole  fraction  due  to  the  changes 
in  the  electron  distribution,  and  hence  electron  impact  rate 
constants  when  the  mixture  is  diluted.  For  example,  when 
diluting  with  argon,  the  average  electron  energy  and  rate 
constant  for  dissociation  increases.  Therefore,  in  spite  of  re¬ 
ducing  the  mole  fraction  of  silane,  the  total  rate  of  dissocia¬ 
tion,  and  hence  deposition  rate,  can  increase  as  more  net 
power  is  coupled  into  the  dissociation  of  silane.  This  trend  is 
illustrated  in  Fig.  1 5,  where  deposition  rate  is  plotted  as  a 
function  of  Ar/SiH4  gas  mixture.  Normalized  experimental 
results  of  Street,  Knights,  and  Biegelsen*^  are  shown  for 
comparison.  The  maximum  in  deposition  rate  at  SiHy 
Ar  s;  1  / 1  is  a  result  of  the  tradeoffbetween  an  increasing  rate 
constant  for  dissociation  and  decreasing  mole  fraction  of  si¬ 
lane. 
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GAS  MIXTURE 

FIG.  1  i.  Deposition  rate  of  a-Si:H  for  silane/argon  gas  mixtures,  compared 
to  the  normalized  experimental  results  of  Street  era/.  (Ref.  87). The  depo¬ 
sition  rate  optimizes  as  a  tradeoff  between  increasing  electron  impact  disso¬ 
ciation  rate  constant  and  decreasing  silane  mole  fraction. 


E.  Formation  of  higher  allane  Ions 

The  formation  of  particulate  matter  (synonomous  with 
clusters  and  dust)  is  relevant  with  respect  to  the  detrimental 
effect  particulate  matter  is  believed  to  have  on  the  quality  of 
films.  Particulates  are  often  observed  and  their  spatial  distri* 
butions  have  been  measured  by  laser  light  scattering.*'  The 
formation  of  particulates  is  generally  associated  with  high- 
power  deposition  and  a  high  partial  pressure  of  silane.  Low 
him  quality  under  these  conditions  is  attributed  to  silicon 
“flakes”  which  are  incorporated  into  the  him  as  polymeric 
groups  (SiHi),.*’ 

When  quantitatively  observed,  particulates  usually 
have  a  higher  densities  adjacent  to  the  sheaths,  leading  one  to 
assume  that  they  are  charged.  One  then  requires  a  mecha¬ 
nism  to  explain  the  growth  of  higher  silane  ions,  presumably 
the  precursor  to  charged  particulates.  Silane  ions  having 
greater  than  6-8  silicon  ions  have  been  observed  in  dc  dis¬ 
charges  by  Weakliem,  Estes,  and  Longeway.^^  Chatham  and 
Gallagher  have  calculated  that  Si„H„|'  ions  with  n  >  2  are 
the  most  abundant  ions  in  dc  negative  glow  discharges  for  j/ 
p  <0.05  A/cm^  Torr."  We  studied  the  formation  of  higher 
silane  ions  by  comparing  results  from  our  model  to  the  ex¬ 
periments  of  Weakliem  et  al.,  an  example  which  is  repro¬ 
duced  in  Fig.  16(a)  for  a  silane  pressure  of  0.5  Torr  and 
current  density  of  13  mA  cm~^.  In  this  figure,  the  fractional 
abundances  of  Si„H,*  ions  having  m<6  arriving  at  the  an¬ 
ode  of  a  dc  discharge  are  plotted  as  a  function  of  electrode 
separation.  Note  that  the  dominant  silane  ions  have  >4  sili¬ 
con  atoms  and  that  the  least  abundant  ion  is  mono-silicon. 

Pathways  for  the  formation  of  higher  silane  ions  are  il¬ 
lustrated  in  Fig.  1 7.  The  rate  of  loss  of  Si,  ions  by  ion-ion 

neutralization,  dissociative  recombination,  and  diffusion 
have  typical  timescales  of  2-8  ms.  The  formation  of  higher 
silane  ions  by  ionization  of  higher  silane  radicals  has  a  typi¬ 
cal  timescale  of  50  ms.  Clustering  of  ions  with  radicals  of  all 
description  has  a  time  scale  of  approximately  2  ms.  For  high¬ 
er  silane  ions  to  be  the  most  abundant  ion,  one  must  invoke 


FIG.  16.  Fractional  densities  oflower  and  higher  silane  ions  incident  on  the 
anode  of  a  dc  discharge  in  silane  as  a  function  of  distance  between  cathode 
and  anode  ( 13  mA/cm^,  O.S  Torr);  (a)  experimental  results  from  Weak¬ 
liem  etal.  (Ref.  22),  (b)  results  from  the  model. 


FIG.  1 7.  Pathways  and  reaction  times  for  creating  higher  silane  ions  for 
typical  discharge  conditions.  The  values  in  parenthesis  are  the  characteris¬ 
tic  times  for  the  indicated  loss  process  ( ion-ion  neutralization,  diffusion,  or 
recombination)  or  for  the  ion-molecule  reaction  leading  to  clustering.  The 
entry  for  clustering  by  ion-molecule  reactions  with  silane  has  a  characteris¬ 
tic  time  or0.02  ms/k,  where  the  rate  constant  is  expressed  in  10  ~ '  cm '  s  ' . 
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ion-tnolecule  reactions  with  silane,  as  only  this  mechanism 
has  a  typical  time  scale  less  than  that  for  losses  when  a  rea^ 
sonable  rate  constant  (/c<  10“  cm’  s“ ' )  is  used. 

The  electron  kinetics  and  plasma  chemistry  models 
were  modified  to  simulate  a  dc  discharge  and  were  used  to 
model  the  experiments  of  Weakliem  et  al.  We  hypothesized 
that  higher  silane  ions  are  formed  by  consecutive  ion  mole* 
cule  reactions  of  the  form 

Si.H.:  +SiH,^Si,,,H„%,  +Hj  (16) 

and  parametrized  the  rate  constant  for  this  class  of  reactions. 
The  results  of  the  study  are  shown  in  Fig.  16(b),  which 
should  be  compared  to  the  experimental  results  of  Fig. 
16(a).  The  systematic  agreement  is  quite  good,  showing  that 
for  similar  conditions  (silane  partial  pressure  >0.3-0.5 
Torr,  electrode  separation  2-4  cm),  the  majority  of  ions  are 
higher  silanes.  To  obtain  the  results  in  Fig.  14(b),  the  chain 
propagating  rate  constant  for  reaction  16  was  1x10“" 
cm’  s“ ',  with  the  exception  of  the  reaction  with  n  =  3  where 
the  rate  constant  was  2  X  10“  ‘ '  cm’  s“ '.  We  also  found  it 
necessary  to  reduce  the  three-body  rate  constant  for  the  reac¬ 
tion  SiHj^ -t-SiH^H-M-SUH,*  -l-M  from  1.8x10“’* 
cm^s"'  to  1.0x10“’*  cm*s“'.  Note  that  if  the  reaction 
sequence  of  Eq.  ( 16)  holds  beyond  n  =  5,  the  detection  limit 
of  Weakliem's  experiments,  the  absolute  density  of  all  ions 
having  greater  than  6  silicon  atoms  exceeds  that  of  the  lesser 
ions  for  large  electrode  separations. 

These  results  only  qualitatively  agree  with  the  recent 
measurements  by  Mandich  et  alJ*  who  obtained  rate  con¬ 
stants  for  the  clustering  chain  of  cold  ions  (300  k)  given  by 


Si-"  -t-SiD^-SiD4Si2Di"  +D2 

SiD^I 


L-SijD/  -I-  D, 
SiD4  I 


+  Dj Si5Dg^  +  D2. 

(17) 


plasma  and  on  the  initiating  step  of  the  chain  (e.g., 
SiH,-*-  -f  SiH4— products,  0</i <3 ). 

There  are  many  implications  of  the  observation  that 
ions  in  silane  discharges  of  moderate  electrode  separation 
and  pressure  are  dominantly  higher  silanes.  Since  the  reac¬ 
tion  chain  is  first  order,  the  formation  of  higher  silane  ions 
depends  only  upon  the  pressure  x  electrode  separation  (pd) 
product  and  not  upon  having  higher  power  deposition. 
Therefore,  the  formation  of  large  silane  ions  can  be  prevent¬ 
ed  only  by  operating  at  small  pd.  Since  large  ions  are  less 
efficient  at  transferring  momentum  to  individual  species 
when  they  collide  with  the  substrate,  the  rate  of  sputtering 
and  activation  of  surface  sites  will  be  lower  with  higher  silane 
ions.  To  the  extent  that  site  activation  may  be  required  for 
incorporation  of  radicals  such  as  SiH3  into  passivated  sur¬ 
face  bonds  (i.e.,  sSi-H),  this  change  in  the  rate  of  surface 
activation  will  change  both  film  growth  rates  and  surface 
morphology. 

iV.  CONCLUDING  REMARKS 

A  plasma  chemistry  model  for  the  deposition  of  amor¬ 
phous  silicon  from  rf  and  dc  discharges  sustained  in  gas  mix¬ 
tures  containing  silane  has  been  described,  and  results  from 
the  model  discussed.  A  compendium  of  rate  constants  for 
silane  plasma  chemistry  has  also  been  presented  and  impor¬ 
tant  reactions  discussed.  We  conclude  that  the  spatial  distri¬ 
bution  of  the  production  of  radicals  and  the  subsequent 
probability  of  survival  to  reach  the  substrate  of  those  radi¬ 
cals  greatly  impact  the  characteristics  of  the  a-Si:H  film. 
Geometry  will  also  impact  the  quality  of  the  film  through 
changes  in  the  manner  of  power  deposition  and  the  forma¬ 
tion  of  higher  silane  species.  Diluting  silane  with  either  inert 
or  active  gases  changes  both  the  manner  of  power  deposition 
and  the  reaction  chain,  thereby  affecting  film  properties.  The 
formation  of  higher  silane  ions  requires  first-order  reactions 
with  silane  and  also  may  affect  film  properties.  Clearly,  even 
small  changes  in  discharge  parameters  may  impact  the  qual¬ 
ity  of  a-Si:H  films  produced. 


Their  measurements  indicate  a  bottleneck  in  the  step 
between  Si4D^"^  and  Si,D,^ ,  with  the  rate  constant  decreas¬ 
ing  by  a  factor  of  >  2(X)  from  the  previous  step.  To  reproduce 
experimental  results,  we  must  sidestep  this  bottleneck.  If  the 
bottleneck  at  n  =  5  is  due  to  an  activation  energy  barrier, 
electronically  and  vibrationally  excited  species  in  the  plasma 
may  provide  the  energy  reservoir  required  to  relieve  the  bott¬ 
leneck,  a  reservoir  not  available  during  Mandich’s  measure¬ 
ments.  In  this  respect,  subsequent  measurements  by  Reents 
and  Mandich^'  have  shown  that  the  rate  constants  for  ion- 
molecule  clustering  reactions  of  Si,  D,"  with  SiD4  are  sensi¬ 
tive  to  the  internal  energy  of  the  reactant  ion,  as  well  as  on  its 
structural  isomeric  state.  For  example,  structural  isomeric 
pairs  of  Si4D2"  either  react  rapidly  with  SiD4  to  form  8^04^ 

( rate  constant  1.3x  10“"’ cm"’ s“')  or  have  no  measura¬ 
ble  reaction.  Reaction  rates  for  sets  of  Si,D,;|^  ions  (n  fixed, 
m  variable)  with  SiD4  were  also  found  to  have  significantly 
different  rate  constants.^  These  results  imply  that  the  prop¬ 
agation  length  of  an  ion-molecule  clustering  chain  may  criti¬ 
cally  depend  on  the  availability  of  activation  energy  in  the 
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SIMULATION  OF  THE  BULK  AND  SURFACE 
PROPERTIES  OF  AMORPIKXJS  HYDROGENATED  SILICON 
DEPOSITED  FROM  SILANE  PLASMAS 


Michael  J.  McCaughey  and  Mark  J.  Kushner 
University  of  Illinois 

Department  of  Electrical  emd  Computer  Engineering 
Gaseous  Electronics  Laboratory 
607  E-  Healey 
Champaign.  Illinois  61820 


ABSTRACT 

A  Monte  Carlo  simulation  for  the  growth  of  amorphous  hydrogenated  silicon 
(a-Si:H)  thin  films  by  plasma  enhanced  chemical  vapor  deposition  (PECVD)  Is 
presented.  The  goal  of  the  model  is  to  predict  the  bulk  and  surface 
properties  of  films  (eg.,  hydrogen  content,  deposition  rate,  buried 
hydride/di hydride  ratios,  porosity,  and  surface  roughness)  having  thicknesses 
of  £2000  A.  The  effects  on  the  film  properties  of  the  composition  of  the 
radical  flux  incident  on  the  surface  are  examined.  Film  properties  were  found 
to  be  critically  dependant  on  the  ratio  of  SiH.^/SiH2  in  the  radical  flux. 

High  values  for  this  ratio  results  in  film  properties  resembling  chemical 
vapor  deposition.  Film  properties  obtained  with  small  values  resemble 
physical  vapor  deposition.  Rough  films  (RMS  roughness  >10’s  of  A)  accordingly 
result  from  radical  fluxes  having  high  SiHg  fractions.  We  find  that  surface 
roughness  and  hydrogen  fraction  increase  with  increasing  growth  rate  and 
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Increasing  film  thickness,  though  thin  films  (<10's  layers)  have  large 
hydrogen  fractions  due  to  there  being  a  hydrogen  rich  surface  layer  with  a 
thickness  approximately  equal  to  the  surface  roughness.  We  also  find  a 
correlation  between  porosity  (subsurface  voids)  emd  hydrogen  fraction, 
implying  that  the  inner  surfaces  of  voids  are  lined  with  ^i-H  configurations. 
In  simulating  films  grown  from  Ar/SiH^  gas  mixtures,  we  find  that  a  decrease 
in  the  SiH2/SiH2  ratio  in  the  radical  flux,  and  Increase  in  the  ion/radical 
ratio  incident  on  the  surface  are  largely  responsible  for  the  degradation  of 
film  properties  observed  with  decreasing  sileme  fraction. 
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I.  Introduction 

Thin  films  of  hydrogenated  amorphous  silicon  (a-Si^H)  fabricated  by 

plasma  enhanced  chemical  vapor  deposition  (PECVD)  are  valuable  materials  for 

the  manufacture  of  thin  film  transistors,  direct-line  contact  sensors, 
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xerographic  materials,  and  solar  cells.  The  efforts  of  many  researchers 

have  resulted  in  a  generally  accepted  working  theory  that  device  quality  films 

-5  -1 

(photoconductivity  >  10  (fl-cm)  .  low  surface  roughness,  hydrlde/dlhydride 

16  ““3 

ratio  >1.  spin  density  <  10  cm  )  are  obtained  from  pure  silane  plasmas 
having  radical  fluxes  Incident  on  the  substrate  composed  dominantly  of  silyl 
(SiH^).  The  mechanisms  responsible  for  this  behavior  were  proposed  by 
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Gallagher.  He  observed  that  silyl  cannot  directly  Insert  into  hydrogen 
passivated  (^1-H)  silicon  bonds  on  the  surface  and  as  a  result  they  are 
highly  mobile  in  the  adsorbed  state.  The  residence  time  and  sticking 
coefficient  of  these  radicals  then  depends  on  the  generation  of  dangling  bonds 
on  the  surface  by  ion  bombardment  or  etching  reactions.  The  diffusive  motion 
of  these  adsorbed  radicals,  filling  in  surface  roughness  results  in  smooth 
films.  Plasma  conditions  resulting  in  large  fractions  of  silylene  (SiHg)  in 
the  radical  flux,  which  can  directly  insert  into  passivated  surface  sites,  or 
conditions  which  result  in  a  high  rate  of  activating  surface  sites  by  removing 
passivating  H  atoms,  result  in  rougher  poorer  quality  films. 

In  simulating  the  deposition  of  thin  films  of  a-Si:H  by  PECVD  one  desires 
to  predict  physical  and  processing  parameters  such  as  growth  rate,  hydrogen 
fraction,  hydrlde/dihydrlde  ratio,  local  bond  structure  (Si-Si  and  Si-H  bond 
angles),  surface  roughness  and  porosity.  An  extension  of  this  predictive 
capability  would  be  to  relate  these  physical  pareuneters  to 
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electronic  properties  of  the  material,  either  on  an  ab  initio  or 

semi -empirical  basis.  Although  the  qualitative  description  of  film  growth 

described  above  is  logical,  there  are  few.  if  any,  quantitative  models  for 

film  growth  which  can  predict  these  properties  based  on  the  cited  mechanisms. 

In  a  similar  fashion,  many  models  and  descriptions  of  the  plasma  kinetics  of 
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silane  discharges  have  been  presented  however  few  have  correlated  plasma 
properties  to  film  characteristics.  To  date.  Gleason  et  al . has  presented 
the  most  complete  physical  model  of  film  growth  based  on  first  principles 
simulation  ’echnlques.  Their  model  is  capable  of  predlcitlng  bond  angles, 
hydride  ratios,  hydrogen  fraction  and  local  order.  Their  work,  though,  did 
not  address  thick  films  having  many  hundreds  of  layers.  Therefore,  bulk 
properties  such  as  surface  roughness  and  porosity  were  not  predicted. 

In  this  paper,  we  present  a  Monte  Carlo  simulation  for  the  deposition 
of  a-Si:H  from  silane  glow  discharges  and  we  relate  plasma  properties  (in  the 
form  of  radical  and  ion  fluxes  to  the  substrate)  to  bulk  and  surface 
characteristics  of  the  film.  As  described  below,  we  do  not  address  the  issue 
of  short  range  order  (eg.,  bond  2mgles}  as  discussed  by  Gleason  et  al , 
though,  we  do  address  bulk  film  properties  such  as  porosity  and  surface 
roughness.  By  using  computed  radical  fluxes  from  a  companion  model  for  the 
plasma  chemistry  of  silane  discharges,  a  complete  first  order  model  for  the 
growth  of  a-Sl:H  films  is  presented. 

The  model  presented  here  emphasizes  the  bulk  properties  of  the  film  as 
opposed  to  local  bond  structure.  The  modeling  techniques  we  use  are  capable 
of  simulating  the  deposition  of  thick  films  exceeding  many  hundreds  of  layers 
while  de-emphaslzing  the  details  of  the  short  range  order,  particularly  bond 
angles.  This  tradeoff  is  accomplished  by  assuming  that  the  film  is  composed 
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of  a  cubic  lattice,  where  each  lattice  point  may  be  occupied  by  either  a 
silicon  atom  or  a  void.  The  use  of  the  term  lattice  does  not  imply  that  there 
is  any  long  range  order  In  the  occupation  probability  of  a  site  by  a  silicon 
atom.  No  assumptions  are  made  concerning  the  coordination  partners  of  each 
silicon  atom.  The  amorphous  character  of  the  film  results  from  the  random 
distribution  of  coordination  partners  for  each  silicon  atom  in  the  lattice, 
and  from  the  finite  but  random  probability  that  a  lattice  point  is  occupied  by 
either  a  silicon  atom  or  a  void. 

In  Section  II,  the  Monte  Carlo  model  for  film  growth  is  described  in 
detail.  Simulated  properties  of  a-Sl;H  films  (roughness,  hydrogen  content, 
porosity)  deposited  from  silane  plasmas  are  discussed  in  Section  III  using 
parametric  radical  fluxes.  The  simulation  of  PECVD  a-Si-'H  films  grown  from 
Ar/SIH^  gas  mixtures  is  discussed  in  Section  IV  followed  by  concluding  remarks 
in  Section  V. 

II.  Description  of  the  Model 

Our  model  for  the  growth  of  a-Si:H  films  uses  Monte  Carlo  techniques  to 
simulate  the  impingment  and  adsorption  of  gas  phase  radicals  emd  ions  on  the 
surface  of  a  growing  film,  as  well  as  for  the  processes  which  occur  after 
adsorption.  These  processes  Include  the  sputtering  of  surface  species  by 
energetic  ions,  diffusion  of  adsorbed  radicals  on  the  surface,  the 
incorporation  of  radicals  into  the  film,  the  desorption  of  radicals  from  the 
surface,  and  the  elimination  of  hydrogen  from  the  film  by  cross  linking 
between  Sl-H  bonds  in  the  near  surface  layers  of  the  film.  Gas  phase-surface 
reactions  are  also  simulated  by  Including  reactions  between  the  impinging  flux 
and  the  surface  lattice  species  when  non-sticking  radicals  strike  the  surface. 
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For  the  results  discussed  here,  we  typically  modeled  the  growth  of  a  square 
patch  of  a-Si:H  approximately  100-250  A  on  a  side  having  up  to  1000  layers  or 
a  thickness  of  2000-3000  A.  Periodic  boundary  conditions  are  used  for  the 
lateral  dimensions  of  the  film.  In  the  discussion  that  follows,  the  term 
"site"  denotes  the  x-y  coordinate  of  a  position  on  the  film  as  viewed  normally 
from  the  plasma.  The  term  "lattice  point"  denotes  the  position  of  a 
particular  atom  or  void  in  the  film  below  a  given  site.  There  are  nany 
lattice  points  stacked  vertically  at  each  site. 

The  simulation  proceeds  schematically  as  follows.  A  flux  of  radicals  is 
directed  towards  the  surface,  isotropically  for  neutrals  and  vertically  for 
ions.  The  identity  of  the  individual  radical  and  the  local  morphology  of  the 
impingement  site  determines  whether  the  radical  adsorbs,  sputters,  bonds,  or 
chemically  reacts  with  the  surface.  An  adsorbed,  but  non-bonded,  radical  may 
diffuse  on  the  surface  until  encountering  an  activated  site  (ie.,  a  dangling 
bond,  sSi-*)  at  which  time  it  bonds,  a  process  which  constitutes  film  growth. 
SiH^  radicals  which  bond  to  the  film,  may  later  Interconnect  with  adjacent 
surface  species  to  form  the  amorphous  network.  Cross  linking  may  occur 
between  both  nearest  and  next  nearest  neighbors,  the  latter  possibly  resulting 
in  the  formation  of  voids.  Adsorbed  species  exceeding  a  specified  residence 
time  without  bonding  to  the  film  may  desorb,  thereby  enabling  calculation  of  a 
sticking  coefficient.  By  specifying  characteristic  frequencies  for  each 
process,  the  rate  of  occurance  of  a  particular  event  is  determined  relative  to 
the  rate  of  arrival  of  radicals  on  the  surface.  The  model  then  proceeds  as  an 
effective  integration  in  time  of  the  equation  of  motions  of  each  particle  in 
the  simulation. 

An  important  feature  of  the  model  is  the  categorization  of  lattice  points 
as  being  either  in  the  near  surface  layer  or  in  the  bulk  film.  Since  we 
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desire  to  simulate  volumes  of  a-Si:H  films  having  10^  -  10^  individual  sites, 
it  is  not  computationally  practical  to  retain  statistical  information  on  each 
site.  The  model  instead  retains  statistics  for  each  lattice  point  in  the  near 
surface  layer  at  each  site.  The  near  surface  layer  is  typically  5-20  atoms 
deep.  When  a  radical  bonds  to  the  surface  at  the  top  of  the  near  layer,  the 
lattice  point  at  the  bottom  of  the  near  layer  is  removed  tind  its  statistical 
information  added  to  the  average  values  for  the  bulk  film.  Although  the  near 
surface  layers  at  each  site  are  only  5-20  layers  thick,  the  height  of  film  at 
adjacent  sites  may  differ  by  an  arbitrary  amount.  Therefore,  the  rms 
difference  in  height  (or  roughness)  of  the  film  may  exceed  the  near  surface 
layer  thickness. 

A.  Lattice  Points  and  Site  Heights 

Lattice  points  in  the  near  surface  layers  are  classified  as  being  either 
unoccupied  or  occupied  by  a  silicon  atom.  An  occupied  site  is  identified  by 
its  coordination  partners.  For  example  a  point  identified  as  SSSS  denotes 
occupation  by  a  silicon  atom  having  four  other  Si  atoms  as  coordination 
partners.  A  point  identified  as  SSHD  has  two  Si  atoms  and  an  H  atom  as 
coordination  partners,  and  a  dangling  bond.  Lattice  points  at  a  particular 
site  which  are  unoccupied  and  buried  beneath  an  occupied  point  are  denoted  as 
a  void.  The  fraction  of  voids  beneath  the  surface  layer  is  therefore  an 
Indication  of  the  porosity  of  the  film.  The  topmost  occupied  lattice  point  is 
the  surface  of  the  film  at  that  site.  The  average  thickness  of  the  film  is 
then  the  average  of  the  topmost  heights  for  each  site.  The  roughness  of  the 
film  is  the  rms  deviation  of  the  thickness  of  the  film  from  its  average  value. 


8 


B.  Radical  Fluxes 

The  flux  of  neutral  radicals  Incident  on  the  surface  may  consist  of  an 
arbitrary  mix  of  SiH^  (0  1  n  ^  3)  or  H.  We  do  not  presently  consider  disllane 
radicals.  The  flux  of  ions  is  assumed  to  consist  of  only  a  single  generic 
charged  species.  The  distribution  of  radicals  is  either  specified  as  a  part 
of  a  parametric  survey,  or  is  obtained  from  the  results  the  electron  kinetics 
and  plasma  chemistry  model  described  in  Ref.  6.  The  time  interval  between 
impingement  of  a  neutral  radical  from  the  plasma  onto  the  surface  is 


Atj^  =  (1) 

where  i  is  the  total  neutral  radical  flux  and  L  is  the  physical  length  of  the 
square  of  film  we  are  simulating.  The  length  is  based  on  the  width,  in  number 
of  sites,  of  the  lattice  assuming  an  average  lattice  constant  of  2.5  A.  Every 
At^  during  the  simulation  a  radical  is  directed  towards  the  surface  based  on 
selection  of  a  random  number  r  (0  <  r  $  1).  The  particular  radical  which  is 
directed  towards  the  surface  is  that  radical  satisfying 

F  <  r  ^  F  F  =  Z  f  (2) 

i  '  J  $  i  J 

where  f^  is  the  mole  fraction  of  radical  i  in  the  flux  of  radicals  Incident  on 
the  surface.  A  similarly  defined  timer  increment,  Atj,  for  the  ion  flux 
determines  the  frequency  with  which  ions  are  directed  towards  the  surface. 

The  site  on  the  surface  (x^.y^)  at  which  the  selected  radical  impacts  is 
randomly  selected  by  use  of  two  additional  rtindom  numbers: 
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Xj  =  ri*L.  Yj  =  r2*L.  The  particular  site  is  constrained  in  the  following 
fashion.  Assuming  that  the  neutral  radical  flux  is  incident  isotropically 
onto  the  surface,  the  probability  that  a  radical  impacts  in  a  deep  well 
(resulting  from  surface  roughness)  is  small.  The  local  morphology  of  the 
selected  site  is  therefore  Inspected  to  determine  whether  the  site  sits  at  the 
bottom  of  a  well.  A  well  is  defined  as  a  site  which  is  lower  than  the  height 
of  its  nearest  neighbor  sites  by  an  amount  greater  than  or  equal  to  the  near 
surface  layer.  If  this  is  the  case,  the  radical  which  would  have  otherwise 
Impacted  in  the  well  is  placed  on  the  surface  site  adjacent  to  the  well. 
Assuming  that  ions  are  incident  vertically,  they  may  impact  on  any  site 
including  those  in  wells. 

C.  Selection  Criteria  for  Plasma-Surface  Reactions 

The  surface  lattice  points  at  any  site  may  be  classified  as  inert, 
passivated,  or  activated.  An  inert  surface  site  is  occupied  by  a  silicon  atom 
fully  coordinated  to  other  silicon  atoms  (ie..  SSSS).  A  passivated  surface 
site  is  occupied  by  a  silicon  atom  having  at  least  one  Si-H  bond  and  no 
dangling  bonds  (eg.,  SSHH)  whereas  an  activated  surface  site  is  occupied  by  a 
silicon  atom  having  at  least  one  dangling  bond  (eg.,SSSD).  All  neutral  silane 
radicals  incident  on  an  inert  site  are  assumed  to  adsorb  onto  the  surface. 

All  neutral  silane  radicals  incident  on  an  activated  site  may  directly 
incorporate  into  the  lattice;  for  example. 


SiH2(g)  +  SSSD  -+  SSSS  +  SHHD. 


(3) 
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where  (g)  denotes  a  gas  phase  species. 

Neutral  silane  radicals  SlH  .  n  1  2.  may  also  incorporate  into  passivated 

n 

sites  by  displacing  a  hydrogen  atom  from  the  surface  species,  a  process  which 
is  nearly  thermo-neutral.  For  example, 

SlH(g)  +  SSSH  -*  SSSS  +  SHHD.  (4) 

SiH^  radicals  with  n  i  2  therefore  effectively  have  sticking  coefficients  of 
near  unity  unless  later  sputtered  off  the  surface.  In  the  model,  sllyl 
radicals  (SiH^)  may  only  Incorporate  into  an  active  surface  site  and  therefore 
most  often  move  directly  into  an  adsorbed  state  since  the  fractional  surface 
density  of  active  sites  is  small.  This  condition  is  important  to  determining 
the  surface  properties  since  the  mole  fraction  of  SiH^  radicals  in  the 
incident  flux  is  typically  ^  0.9  for  plasma  conditions  which  result  in  device 
quality  films. 

The  rate  of  generation  and  passivation  of  surface  sites  is  important  in 
two  respects.  First,  the  removal  of  hydrogen  from  the  surface  is  Important  in 
the  accounting  of  the  total  amount  of  hydrogen  in  the  film  and  in  determining 
the  characteristics  of  the  hydrogen  rich  surface  layers.  Second,  the  net  rate 
of  generation  of  dangling  bonds  on  the  surface  by  removing  passivating  H  atoms 
in  part  determine  the  sticking  coefficient  for  SiH^  radicals  since  this  specie 
requires  an  activated  surface  site  for  incorporation. 

New  passivated  and  active  surface  sites  are  continually  being  generated 
by  radicals  which  incorporate  into  the  film  and  constitute  the  new  surface 
specie  at  that  site.  Non-sticking  radicals,  molecules,  and  ions  from  the 
plasma  may  also  react  with  the  surface  species  to  generate  activated  sites. 
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Specifically,  In  the  model  we  assume  that  an  H  atom  (In  the  incident  flux) 
landing  on  2ui  active  site  passivates  the  dangling  bond,  whereas  H  atoms 
landing  on  a  passivated  site  extract  a  hydrogen  thereby  etching  the  surface. 
For  example. 


H(g)  +  SSSD  -»  SSSH 
H(g)  +  SSSH  SSSD  +  H2(g) 


The  latter  process  is  analogous  to  the  gas  phase  hydrogen  abstraction  reaction 


H  +  SiH^  SiH3  +  Hg.  (7) 

The  analogous  reactions  to  Eqs.  5  and  6  may  also  occur  with  SlH^  as  the 
reactant  Instead  of  H.  The  former  reaction 

SiH^(g)  +  SSSD  - »  SSSH  +  SiH^Cg)  (8) 

may  be  important  since  the  flux  of  SIH^  striking  the  surface  is  large  compared 
to  that  for  H.  The  etching  reaction 

SlH^(g)  +  SSSH  - »  SSSD  +  SiH^Cg)  +  HgCg)  (9) 

Is  likely  not  to  be  rapid  since  it  is  endothermic  by  approximately  1.9 
kcal/mole  [11].  Nevertheless,  as  will  be  discussed  in  Sec.  IV,  the  net  rate 
of  passivation  of  the  surface  by  SIH^  is  imjxjrtant  in  accounting  for  film 


properties . 
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Dangling  bonds  may  also  be  generated  by  sputtering  off  the  passivating  H 
atom  by  ion  impact  and  by  sputtering  of  loosely  adhered  silicon  containing 
surface  species,  for  example 

M*(g)  +  SSSH  -*  M*(g)  +  SSSD  +  H(g) .  (10) 

The  model  Ignores  the  direct  contribution  of  ion  impact  to  the  densif ication 
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of  the  film,  as  addressed  by  Mueller  and  by  Drevillon.  The  precise 
mechanism  for  hydrogen  removal  by  ion  impact  is  not  specified  in  the  model, 
though  it  is  assumed  to  be  a  combination  of  sputtering  and  thermal  spiking. 

In  the  simulation,  a  maximum  number  of  H  atoms  is  specified,  typically  10, 
that  may  be  removed  from  passivated  bonds  per  incident  ion  subject  to  the 
availability  of  =S1-H  bonds  in  the  vicinity  of  the  ion  impact.  In  practice, 
the  average  number  removed  per  Incident  ion  is  less  than  ten  due  to  this 
requirement . 

Ion  impact  may  also  sputter  off  loosely  adhered  SlH^  groups  on  the 
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surface,  a  process  which  has  been  described  as  "scouring"  of  the  surface. 

Typically,  these  groups  are  the  terminous  of  polymeric  chains  which  are  not 

presently  simulated  in  this  model.  To  investigate  the  importance  of  this 

process,  we  defined  loosely  adhered  SiH  groups  as  species  on  the  surface 

n 

which  have  only  1  silicon  atom  as  a  coordination  partner  (ie:  SHHH.  SHHD, 
SHDD,  SODD) .  We  then  included  the  process,  for  example, 

M'^Cg)  +  (SSSS-SHHH)  — ♦  +  SSSD  +  SiH^Cg)  (11) 

in  the  model.  The  importance  of  this  process  is  discussed  in  Section  IV. 
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D.  Surface  Diffusion  and  Reactions 

During  the  simulation,  adsorbed  radicals  on  the  surface,  principally 
SiH^.  diffuse  from  surface  site  to  surface  site.  Adsorbed  radicals  are  moved 
from  their  current  surface  site  with  each  increment  of  the  global  timer.  The 
direction  of  the  move  in  the  plane  of  the  surface  (one  of  eight  directions 
including  nearest  and  next  nearest  neighbors)  is  determined  by  choice  of  a 
random  number.  An  adsorbed  radical  is  allowed  to  make  a  specified  number  of 
such  moves  before  being  removed  from  the  surface;  that  is,  being  desorbed. 

The  specified  number  of  moves  is  an  input  parameter  chosen  to  give  a 
reasonable  sticking  coefficient.  Under  conditions  which  experimentally  yield 
high  quality  films,  allowing  10  such  moves  for  SiH^  results  in  a  sticking 

1C  1C 

coefficient  of  0.05  -  0.25,  in  fair  agreement  with  recent  measurements.  ’  ’ 

The  diffusive  motion  of  the  adsorbed  radical  is  constrained  in  the 
following  manner.  The  height  of  the  film  at  a  particular  site  is  h.  When  an 
adsorbed  radical  is  moved,  we  define  the  height  of  the  film  at  the  initial 
site  as  hj  and  the  height  of  the  film  at  the  site  to  which  the  adsorbed 
radical  is  chosen  to  move  as  h^..  Moves  with  h^  >  h^  are  allowed;  that  is  the 
adsorbed  radical  is  allowed  to  make  a  lateral  move  or  fall  onto  a  site  of 
lower  height.  Moves  with  h^  <  hj.  are  not  allowed;  that  is  the  adsorbed 
radical  is  not  allowed  to  climb  to  a  higher  surface  site.  An  adsorbed  species 
may  be  trapped  in  a  well  as  a  result  of  a  move.  A  trapped  radical  is  one  which 
does  not  move  from  a  surface  site  for  the  maximum  number  of  attempted  moves. 

In  this  case  the  model  assumes  that  the  radical  is  permanently  trapped  and 
will  eventually  bond  at  that  location.  It  is  therefore  bonded  in  place  and  If 
necessary  Hg  is  eliminated. 


When  an  adsorbed  radical  moves  over  a  new  surface  site,  it  may  bond  to 
the  occupying  Si  atom  and  incorporate  into  the  film  in  the  same  manner  as 
described  above  for  neutral  radicals  incident  from  the  plasma.  An  adsorbed 
radical  which  does  not  bond  to  the  surface  site  beneath  it  is  moved  eigain  on 
the  next  cycle.  When  an  adsorbed  radical  moves  to  a  surface  site  where  the 
height  of  the  film  at  an  adjacent  site  is  higher  tham  its  present  location 
these  are  additional  considerations.  Under  those  conditions,  as  well  as  when 
the  radical  attempts  a  move  with  h^  >  h^ ,  the  adsorbed  radical  may  bond  with 
the  lattice  atom  having  the  same  height  at  the  adjacent  site.  This  adjacent 
lattice  atom  is  buried  beneath  its  own  surface  site  cmd  constitutes  part  of 
the  near  surface  layer  for  that  site. 

E.  Cross  Linking,  Interconnection  eind  Hydrogen  Elimination 

Cross  linking  is  the  process  whereby  =S1-*  (eg.,  SSSD)  and  =Si-H  (eg., 
SSSH)  configurations,  bond,  and  form  the  amorphous  network  of  the  film.  Cross 
linking  is  also  an  important  process  for  removing  hydrogen  previously 
incorporated  into  the  lattice.  Cross  linking  may  occur  between  silicon  atoms 
in  adjacent  sites  having  dangling  bonds,  for  example, 

SSSD  +  SSSD  -»  SSSS  +  SSSS  (12) 

or  between  silicon  atoms  in  adjacent  sites  having  passivated  bonds,  for 


example. 
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The  former  process  is  exothermic.  The  latter  process,  though  being 

endothermic,  is  believed  be  largely  responsible  for  eliminating  hydrogen  from 

the  film.  An  activation  energy  of  1.35  kcal/mole  has  been  assigned  to  this 

process  by  observing  hydrogen  content  as  a  function  of  substrate 
17 

temperature.  Cross  linking  may  occur  between  ^i-H  configurations  topmost 

on  the  surface  and  those  burled  beneath  the  surface  but  exposed  on  a  sidewall. 

a  process  which  constitutes  film  densif ication. 

Cross  linking,  interconnection,  and  hydrogen  elimination  are  simulated  in 

the  following  fashion  in  the  model.  A  characteristic  site  examination 

frequency  for  nearest  neighbors,  u  ,  emd  for  next  nearest  neigbors,  u  are 

n  nn 

specified.  These  frequencies  yield  Intervals,  At^  =  which  specify  the 

average  time  between  Inspecting  adjacent  lattice  points.  After  At^  has 

elapsed  at  a  particular  lattice  point,  a  survey  of  nearest  neighbor  lattice 

points  is  made.  The  nearest  neighbor  lattice  points  are  the  ^1-*  and  =Si-H 

complexes  within  a  lattice  constant  of  the  central  point.  If  the  central 

lattice  point  is  a  ^l-«  complex,  the  survey  consists  of  examining  the  nearest 

neighbor  points  in  a  random  order  searching  for  another  =Si-*  complex.  A 

successful  search  results  in  an  Interconnection  as  shown  in  Eq.  12.  If  the 

central  lattice  point  is  a  =Si-H  complex,  the  survey  searches  for  another 

^i-H  complex.  If  the  search  is  successful,  the  choice  of  another  random 

number  determines  whether  there  is  enough  thermal  energy  available  for  the 

endothermic  process  to  take  place.  This  is  determined  by  having  -!ln(r)  ^ 

e  /T  ,  where  r  is  a  random  number  (0,1).  e  is  the  activation  energy,  and  T 

is  the  substrate  temperature.  If  these  conditions  are  met,  interconnection 

with  hydrogen  elimination  occurs  as  shown  in  Eq.  13.  The  procedure  described 

for  nearest  neighbor  interconnection  is  also  followed  for  next  nearest 

neighbors,  with  the  examination  interval  being  At  =  1/u 

®  nn  nn 
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F.  Void  Formation 

Cross  linking  and  interconnection  with  next  nearest  neighbors  cam  result 
in  the  formation  of  voids  (or  porosity).  The  algorithm  for  these  processes  is 
identical  to  that  described  above  for  interconnection  between  nearest 
neighbors.  The  frequency  of  interconnection  with  next  nearest  neighbors, 
however,  must  be  less  than  that  for  nearest  neighbors  as  evidenced  by  the  low 
porosity  of  films  (<  15%)  which  can  be  deposited.  To  account  for  this 
observation,  a  ratio  of  i  ^  typically  used.  As  cross  llnkii^  with 

at  most  next  nearest  neighbors  is  allowed,  the  voids  generated  in  the  model 
are  at  most  one  lattice  constant  wide.  They  may,  however,  have  either 
vertical  or  horizontal  extents  of  up  to  many  hundreds  of  lattice  constants. 

In  forming  voids  in  this  fashion,  low  porosity  films  may  occur  in  films  which 
have  an  extremely  rough  surface,  since  it  is  not  possible  (In  the  model)  to 
cross  link  at  greater  distances  than  next  nearest  neighbors  which  would 
otherwise  bridge  roughness  eind  form  large  voids.  The  large  surface  roughness, 
though,  gives  the  appearance  of  a  porous  near  layer,  and  in  measurements  of 
surface  roughness  a  void  fraction  Is  often  cited  based  on  the  near  layer 
roughness.  Our  definition  of  porosity  is  a  void  (vacant  lattice  points)  below 
the  topmost  occupied  lattice  point  at  a  site. 

Occupied  surface  lattice  points  on  both  vertical  and  horizontal  surfaces 
are  usually  passivated  by  hydrogen.  Therefore,  the  formation  of  voids  traps 
hydrogen  below  the  surface.  Our  algorithms  do  not  presently  allow  for 
evolution  of  hydrogen  from  sub-surface  sites  and  voids.  Therefore  the  bulk 
hydrogen  fractions  for  porous  material  (void  fraction  >10’s  of  percent)  are 
likely  to  be  overestimated.  This  effect  is  discussed  below. 


III.  Simulated  Properties  of  a-Si:H  Films 

In  this  section,  we  will  discuss  results  obtained  from  our  model  for  the 
simulation  of  the  growth  of  hydrogenated  amorphous  silicon  films  from  sileoie 
plasmas  when  the  radical  fluxes  are  specified  as  part  of  a  parametric  survey. 
In  Section  IV,  we  discuss  results  for  when  the  radical  fluxes  are  obtained 
from  a  companion  model  for  the  electron  kinetics  and  plasma  chemistry  of 
discharges  sustained  in  silane'^. 

Before  beginning  the  parametric  survey,  we  present  a  sample  of  a 
simulated  surface  of  a-Sl:H  for  typical  plasma  deposition  parameters  shown  in 
Fig.  1.  The  radical  fluxes  for  this  case  were  obtained  from  the  electron 
kinetics  and  plasma  chemistry  model.  The  discharge  conditions  are  a  silane 
pressure  of  250  mTorr.  substrate  temperature  of  500  K,  and  power  deposition  of 

3 

100  mW-cm  The  total  mono-silane  radical  flux  incident  on  the  substrate  is 

2.5  X  10^®  cm'^s"^  in  the  ratio  of  SlH^/SiHg/SiH/H  =  0.93/0.0064/1.8  x 

-4  14  -2-1 

10  /0.06.  The  total  ion  flux  incident  on  the  surface  is  1.3  x  10  cm  s 

and  at  most  10  surface  sites  may  be  activated  per  Incident  ion.  The 

frequencies  for  examination  of  nearest  and  next  nearest  neighbors 
-1  -2  -1 

are  u^=  3.9  s  and  =  8x10  s  .respectively.  The  calculated  deposition 
rate  from  the  model  using  these  parameters  is  4.5  A-s‘^.  The  hydrogen 
fraction  of  the  simulated  film  is  fj|  =  12%,  and  the  ratio  of  hydr Ide/dlhydr ide 
bonds  in  the  bulk  material  is  2.1,  which  does  not  include  the  hydrogen  rich 
surface  layers.  By  summing  the  SIH^  flux  incident  onto  the  surface  and 
comparing  the  fluence  to  the  amount  of  SlH^  desorbing,  we  obtain  a  sticking 
coefficient  for  SiHg  of  ~0.12.  The  calculated  surface  roughness  (the  rms 
deviation  from  the  average  film  height  of  750A)  is  approximately  20  A.  All  of 
these  parameters  are  in  good  agreement  with  experimental  results  for  device 
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quality  films  grown  by  PECVD. ^  The  qualitative  appearance  of  the  surface 

agrees  well  with  recent  scanning  tunneling  microscopy  measurements  of  the 

18 

surface  of  plasma  deposited  a-Si^H  films.  Our  predicted  spin  densities  are 

—4  18  -3 

fractionally  large  (a  few  times  10  .  or  ~  5-6  x  10  cm  )  for  films  which 

16 

otherwise  have  good  properties.  Desirable  spin  densities  are  typically  <  10 
-3  1 

cm  We  Interpret  this  discrepancy  as  there  being  thermal  structural 

relaxation  of  the  Si  network  af ter  burial  of  the  surface  species  which  reduces 
the  spin  density,^  a  process  we  do  not  include  in  the  model. 

A.  Surface  Roughness 

The  surface  roughness  of  a-Si:H  films  is  important  with  respect  to  their 

use  in  thin  film  devices  and  in  multilayer  structures  where  the  lattice  period 

19  20 

may  be  comparable  to  the  surface  roughness.  Works  by  Knights,  Tsai,  et 

2 1  22-24  25 

al . ,  Collins  and  Cavese,  and  Collins  and  Pawlowskl  have  resulted  in  a 

qualitative  understanding  of  the  origin  and  scaling  of  surface  roughness.  In 

summary,  plasma  conditions  which  have  a  predominance  of  silane  radicals  which 

have  high  sticking  coefficients  (SlH^,  n  ^  2)  generally  result  in  rougher 

films  euid  poor  sidewall  covere^e.  The  low  surface  mobility  for  these  radicals 

results  in  shadowing  of  precursor  or  random  surface  features,  and  magnify 

their  roughness.  These  conditions  are  found  in  low  pressure  silane  plasmas 

(<100  mTorr)  or  diluted  slleoie  plasmas  where  the  probability  of  radical 

scavanging  reactions,  exemplified  by 


S1H„  +  SIH.  -»  Sl„Hc. 
2  4  2  6 


(14) 
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is  low.  Deposition  under  these  conditions  resembles  physical  vapor  deposition 
(PVD).  It  is  generally  believed  that  plasma  conditions  which  result  in 
radical  fluxes  dominated  by  SiH^  result  in  smoother  films  and  higher  quality 
material,  and  resemble  chemical  vapor  deposition  (CVD) . 

To  investigate  the  effect  of  adsorbant  mobility  on  the  roughness  of 
a-Si:H  films,  we  parameterized  our  model  while  varying  the  fraction  of  SIHq  in 
an  impinging  flux  consisting  only  of  a  mixture  of  SiHg/SiH^.  Those  results 
are  shown  in  Fig.  2.  We  find  that  the  surface  roughness  is  relatively 
constant  until  the  silyl  fraction  is  less  than  0.4,  below  which  surface 
roughness  Increases  markedly.  As  the  growth  rate  is  a  function  of  SiH^ 
fraction  (see  below)  we  performed  our  calculations  for  both  constant  growth 
rate  and  constant  magnitude  of  the  radical  flux,  and  obtained  similar 
roughnesses.  The  ion  flux  was  kept  constant  throughout.  These  results  imply 
that  surface  adsorbant  mobility  is  most  influential  in  determining  the  surface 
roughness.  However,  only  moderate  fractions  of  a  highly  mobile  radical,  such 
as  SlH^.  is  sufficient  to  fill  in  and  smooth  the  surface.  The  sticking 
coefficient  for  SlH^  is  approximately  0.25  for  these  conditions.  It  appears, 
then,  that  the  surface  roughness  generated  by  SiH2  sticking  in  a  PVD-like 
manner  can  be  mitigated  by  having  approximately  >0.1  (0.5  flux  fraction  x  0.25 
sticking  coefficient)  of  the  deposition  result  from  the  higher  mobility 
specie,  SlH^. 

Surface  roughness  is  a  function  of  film  height,  as  well  as  deposition 

22 

specie.  Collins  and  Cavese  found  that  for  conditions  dominated  by  PVD-like 
deposition  surface  roughness  increased  with  increasing  film  height,  whereas 

this  dependence  on  film  height  was  not  as  strong  for  CVD-like  conditions. 

0  3-3  13 

They  found  that  roughness  scales  as  h  ’  '  for  0.2  pm  <  h  <  2  pm.  Drevlllon 
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has  also  found  a  correlation  between  growth  rate  and  surface  roughness.  In 
that  work.  Drevillon  viewed  surface  roughness  as  a  porous  overlayer  on  top  of 
homogeneous  bulk  material. 

We  simulated  a-Si:H  films  with  different  heights  for  CVD  and  PVD-llke 
conditions  by  using  incident  radical  fluxes  of  SlH^/SiH2  =  0.9/0. 1  and 
0.15/0.85,  respectively.  Our  results  are  shown  in  Fig.  3.  For  PVD-llke 
conditions  surface  roughness  scales  almost  linearly  with  thickness.  For 
CVD-llke  conditions  surface  roughness  has  only  weak  dependence  on  film 
thickness  for  thicknesses  greater  than  100  monolayers  (about  250A).  The  poor 
scaling  of  surface  roughness  for  PVD-llke  conditions  (le.,  large  ratios  of 
SiH2/SlH^)  results  primarily  from  the  shadowing  effect  of  precursor  or  random 
roughness  in  the  first  few  layers.  Roughness  increases  with  increasing 
thickness  as  the  shadowing  effect  becomes  more  severe.  This  effect  is  likely 
exaggerated  in  our  simulation  by  not  including  densif ication  by  ion  impact  and 
Interconnection  with  more  than  second  nearest  neighbors  (see  below)  though  the 
roughness  must  clearly  be  greater  than  for  the  CVD-llke  conditions.  For 
CVD-llke  conditions,  a  quasl-equl librium  is  reached  between  shadowing  and  the 
"filling”  effect  of  mobile  surface  species. 

For  CVD-llke  film  growth,  the  deposition  species  with  a  high  surface 
mobility,  SiH^,  requires  a  finite  residence  time  on  the  surface  to  fill  in 
roughness.  This  residence  time  must  be  sufficient  to  allow  the  mobile 
radicals  to  diffuse  a  dlsteince  approximately  equal  to  the  transverse  dimension 
of  the  surface  roughness.  One  would  then  expect  that  the  surface  roughness 
would  increase  with  increasing  deposition  rate  as  the  average  residence  time 
of  mobile  species  on  the  surface  decreases.  In  simulating  films  for  CVD-llke 
conditions,  though,  we  find  only  a  weak  dependence  of  surface  roughness  on 
growth  rate,  as  shown  in  Fig.  4. 
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B.  Hydride/Dihydride  Ratio 

The  ratio  of  burled  hydride  to  dihydrlde  configurations  In  a-Sl-'H  films 

Is  also  an  Indication  of  film  quality,  with  better  films  (e.g. ,  higher 

26 

conductivity)  having  large  (SSSH)/(SSHH)  bond  ratios.  In  our  simulations, 
this  ratio  decreases  with  Increasing  growth  rate  as  shown  In  Fig.  4.  Slower 
growth  rates  allow  longer  residence  times  of  SSHH  configurations  on  the 
surface,  and  therefore  there  Is  more  opportunity  for  interconnection  and 
hydrogen  elimination  which  Increases  the  hyd r 1 de/d 1 hydride  ratio.  Nemy  of  the 
dihydrlde  configurations  are  contained  in  voids  which  (see  below),  in  our 
model,  do  not  further  Interconnect.  Therefore  the  dihydrlde  density  may  be 
over-estimated.  Low  growth  rates  then  appear  to  produce  better  material  based 
on  both  roughness  and  hydrlde/dlhydrlde  ratio  criteria. 

C.  Hydrogen  Fraction  and  Porosity 

Before  discussing  the  relationship  between  hydrogen  fraction  and 
porosity,  we  are  careful  here  to  differentiate  between  surface  roughness  and 
voids.  The  former  is  an  rms  deviation  of  the  surface  height.  The  latter  is 
the  fraction  of  unoccupied  lattice  sites  below  the  surface  emd  Is  thus  a 
measure  of  film  porosity.  We  have  therefore  defined  the  term  porosity  to  mean 
a  true  volumetric  defect  density.  In  the  literature,  the  term  porosity  Is 
often  used  to  describe  a  defect  density  based  on  mean  film  thickness  which 
includes  surface  roughness.  For  example,  the  (volumetric)  void  fraction  of  a 
typical  film  is  10-20%.  A  defect  density  based  on  mean  film  thickness. 


def ined  as 
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R  =  (average  roughness)/ (mean  film  height)  (15) 

would  add  an  additional  defect  density  level  of  about  5%  for  a  film  having  a 
roughness  of  20A  and  height  of  '400A. 

The  atomic  fraction  of  hydrogen  in  a-Si:H  films,  f^j,  is  of  interest  with 

respect  to  its  effect  on  electronic  properties,  such  as  the  optical  bandgap 

and  photoconductivity.  It  is  generally  observed  that  fj^  increases  with 

14  27  28 

Increasing  deposition  rate.  '  ’  This  dependence  has  been  attributed  to  a 

"burial  effect"  where  the  shorter  residence  times  of  ^i-H  configurations  on 

the  surface  resulting  from  Increasing  growth  rates  reduce  the  probability  for 

cross  linking  and  hydrogen  elimination.  We  found  a  similar  dependence  of 

hydrogen  fraction  on  growth  rate,  as  shown  in  Fig.  5;  that  is,  fjjj  increases 

with  increasing  growth  rate.  When  we  do  not  allow  void  creation,  the  same 

general  dependence  of  fjj  on  growth  rate  is  obtained,  however  the  value  of  f^^ 

is  lower  at  slow  growth  rate.  The  increase  in  fj^  at  low  growth  rate 

correlates  directly  with  an  increase  in  the  porosity  of  the  film  with 

decreasing  growth  rate,  also  shown  in  Fig.  5.  This  correlation  between  voids 

and  hydrogen  fraction  has  been  made  experimentally  for  both  a-Si:H  films. 

30 

and  a-Ge:H  films.  We  find  that  large  tunounts  of  hydrogen  are  being  trapped 

in  voids  as  a  result  of  the  sidewalls  being  passivated  by  H  atoms. 

The  Increase  in  void  fraction  we  predict  at  low  deposition  rates  is 

typically  not  observed.  Void  fractions  are  either  constant  with  growth  rate 

28 

or  tend  to  Increase  with  increasing  growth  rates  under  some  conditions. 
Mathematically,  our  result  is  a  consequence  of  the  fact  that  interconnection 
with  next  nearest  neighbors  occurs  at  a  fixed  rate,  so  slower  growth  implies 
more  void  formation.  These  conditions  imply  that  ion  impact  densif Icatlon  and 
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heating  of  the  film,  both  of  which  will  evolve  hydrogen  and  close  voids  af ter 
burial,  must  be  Important. 

The  average  void  sizes  obtained  in  our  simulations  are  approximately 

20-40A  .  We  find  that  the  atomic  hydrogen  content  increases,  and  the  ratio  of 

hydride  to  dihydride  configurations  decreases,  as  the  avereige  void  size 

increases,  (see  Fig.  6).  This  correlation  suggests  that  dihydride 

configurations  are  more  likely  to  be  found  in  voids,  and  by  implication 

hydride  configurations  are  more  uniformly  distributed.  These  predictions 
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agree  with  the  results  of  Reimer,  et.al.  who,  using  NMR  methods,  were  able 

to  differentiate  between  Isolated  monohydride  and  grouped  dihydride 

configurations  in  a-Sl:H.  Since  films  that  have  high  monohydride  fractions 

typically  have  better  electrical  properties,  films  having  small  void  sizes 

should  similarly  have  better  properties. 

The  average  void  sizes  in  Fig.  6  are  small,  <40  A  ,  and  have  equivalent 

spherical  diameters  of  3, 5-4.0  A.  The  distribution  of  voids  is  compact,  with 

-4 

the  relative  density  falling  to  <10  for  voids  having  volumes  of  greater  thm 
a  few  hundred  A  (see  Fig.  7).  Beyond  this  size,  the  distribution  has  a  long 
tall  which  extends  to  larger  volumes.  There  are  infrequent  occurrences  of 
voids  having  volumes  greater  than  many  thousand  A  .  These  large  voids  tend  to 
be  long  fissures,  more  often  vertical  th2Ln  horizontal.  For  comijarison, 
though,  their  equivalent  spherical  diameters  are  <20-25A. 

Void  fractions,  or  porosity,  tend  to  be  constant  for  films  >750  A  thick 
euid  for  deposition  conditions  for  which  SiH^  is  the  predominant  depxjsition 
radical,  as  shown  in  Fig.  8.  For  the  cases  shown,  the  void  fraction  of 
PVD-like  depositions  is  less  th2m  that  for  CVD-like  conditions  when  only 
subsurface  porosity  is  considered.  The  porosity  of  the  PVD-like  material  may 
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be  artificially  low  since  the  model  limits  void  formation  to  those  created  by 
cross  linking  with  next  nearest  neighbors,  which  as  described  above,  is  not 
able  to  bridge  over  large  roug^hness  and  create  large  voids.  Cross  linking  or 
bridging  gaps  across  more  than  one  lattice  point  is  required  to  have  PVD-like 
deposition  have  a  higher  porosity  than  CVD-like  conditions.  When  including 
the  apparent  density  defect  due  to  surface  roughness,  as  given  by  R  in  Eq.15. 
the  defect  density  of  the  PVD-like  films  exceeds  that  of  PVD-like  films  as 
shown  in  Figure  8.  The  majority  of  density  defect  for  thin  PVD-like  films 
(<2000  A)  can  be  attributed  to  surface  roughness. 

D.  Hydrogen  Rich  Surface  Layers 

The  growth  mechanism  for  a-Si:H  films  suggested  by  Gallagher  implies  that 

both  latteral  and  vertical  surfaces  of  the  film  are  saturated  by  ^i-H 

configurations.^  This  mechanism  then  requires  that  there  be  a  hydrogen  rich 

layer  at  the  surface  whose  local  depth  is  Indicative  of  the  depth  of  active 

growth,  and  whose  average  depth  is  approximately  the  surface  roughness. 

Measurements  of  hydrogen  evolution  by  post  deposition  sputtering  of  a-Si:H 

films  have  shown  this  layer  to  be  5  layers  thick  at  room  temperature,  and  a 
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monolayer  at  >500K.  Simulations  of  film  growth  by  Gleason  et  al .  indicated 

that  the  active  layer  is  1.0  -  2.5  atoms  deep^®  not  Including  surface 

roughness.  To  account  for  the  effect  of  this  active  layer  on  fj^,  Gleason 

gives  the  expression  f„  =  N.  f.  +  N  f  ,  where  the  subscripts  denote  bulk 

n  DO  s  s 

material  (b)  or  surface  layer  material  (s).  N  is  the  fraction  of  atoms  in  a 

particular  region,  and  f^  is  the  hydrogen  fraction  in  region  i.  From  our 

simulations  we  find  that  f  is  0. 5-2.0  depending  on  conditions.  The  local 

s 

thickness  of  the  hydrogen  rich  active  region  in  our  model  is  <4-5  layers  as 
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our  results  computed  hydrogen  fractions  are  Insensitive  to  the  depth  of  our 
near  layer  for  values  >p. 

Extending  this  logic,  the  hydrogen  content  of  the  film  should  also  scale 
proportionally  to  Ah/h,  where  Ah  is  the  surface  roughness,  since  more 
roughness  provides  more  exposed  surfaces  which  are  passivated  by  H  atoms. 

Since  the  ratio  Ah/h  decreases  with  increasing  film  thickness,  one  would 
expect  that  the  hydrogen  content  of  thin  films  to  be  higher  than  that  of  thick 
films.  Simulated  hydrogen  content  as  a  function  of  film  thickness  is  shown  in 
Fig.  9.  In  very  thin  films  (less  than  a  few  hundred  A’s)  fj^  approaches  50% 
since  the  exposed  surface  areas  contain  a  sigpiif leant  fraction  of  all  atoms  in 
the  lattice.  For  moderately  rough  films,  we  have  found  that  up  to  30% 
(roughness  40  A,  height  750  A)  of  the  total  atomic  hydrogen  present  in  the 
films  may  be  in  the  near  layer  region  and  within  ~25  A  of  the  average  surface. 
PVD-llke  conditions  result  in  a  hlg^her  hydrogen  fraction  due  to  their  higher 
growth  rate  and  larger  surface  roughness. 

IV.  a-Si'H  Films  from  Ar/SiH.  Plasmas 

4 

It  has  generally  been  observed  that  the  properties  of  a-Si’-H  films 

degrade  upon  dilution  of  the  input  silane  gas  with  argon.  The  films  become 
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rougher,  begin  growing  with  a  columnar  structure,  emd  the  fraction  of 

14 

hydrogen  found  in  voids  or  poly-hydride  configurations  increases. 

29 

Concurrently,  deposition  rate  and  hydrogen  fraction  increase.  To  simulate 
a-Si:H  films  using  realistic  deposition  parameters,  we  computed  radical  fluxes 
using  the  plasma  chemistry  model  described  in  Ref.  6  for  a  parallel  plate  rf 
PECVD  reactor.  The  deposition  conditions  were  250  mTorr  gas  pressure  and  a 
substrate  temperature  of  500  K  with  an  Ar/SlH^  gas  mixture.  The  power 
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-3 

deposition  was  varied  about  50  mW-cm"  in  order  to  give  the  same  flux  of 
sileuie  radicals  onto  the  surface. 

Computed  flux  characteristics  are  shown  in  Fig.  10a  as  a  function  of 
silane  fraction  in  argon.  As  the  mixture  is  diluted  with  argon,  the  ratio  of 
SlH2/SiH2  in  the  incident  flux  decreases,  and  the  ratio  of  the  total  ion  flux 
to  silane  radical  flux  increases.  The  SlH3/SiH2  ratio  decreases  with  argon 
dilution  due  to  the  lower  net  rates  of  radical  scavenging  reactions  such  as 
that  in  Eq.  14  and 


H  +  SIH^  -*  SiH^  +  Hg.  (16) 

both  of  which  combine  to  Increase  the  SiH^/SiHg  ratio.  The  ion  to  silane 

radical  ratio  increases  due  to  there  being  a  larger  fraction  of  power  coupled 

into  argon  which  produces  dominantly  ions. 

Simulated  film  properties  for  PECVD  of  a-Sl:H  in  Ar/SlH^  mixtures  with  a 

constant  magnitude  for  the  radical  flux  are  shown  in  Fig.  10b.  We  find  that 

the  bulk  hydrlde/dihydride  ratio  decreases,  deposition  rate  Increases  and 

hydrogen  fraction  Increases  when  pure  silane  is  diluted  with  argon,  in 

14 

agreement  with  experiments.  The  as  buried  spin  density  also  increases  with 

argon  dilution,  in  spite  of  increased  hydrogenation  of  the  film.  This  effect 

is  likely  a  result  of  the  higher  rate  of  arrival  of  SiH2  and  higher  growth 

rate  which  reduces  the  time  for  Interconnection  on  the  surface.  Film 

roughness  remains  fairly  constant  for  these  conditions  at  320A.  This 

14 

systematic  agreement  with  experiment  was  obtained  by  limiting  the  role  of 
ions  to  sputtering  loosely  adhered  SiH^  groups,  euid  having  the  net  rate  of 
passivation  of  the  surface  be  a  constant.  These  results  imply  that 
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net  rate  of  passivation  by  gas  phase  species  is  dominated  by  SiH^  but  rate 
limited  by  an  activation  energy  barrier  for  the  reverse  process  which  produces 
dangling  bonds.  One  should,  therefore,  find  a  more  sensitive  dependency  of 
p>asslvation  on  substrate  temperature  than  silane  partial  pressure. 

V.  Concluding  Remarks 

A  model  for  the  growth  of  a-Si:H  thin  films  from  silane  plasmas  has  been 
presented.  From  the  results  of  the  model,  we  find  that  the  physical 
properties  of  the  film  depend  on  both  the  relative  fraction  of  the  radicals  in 
the  incident  flux  and  on  the  rate  at  which  the  radicals  impinge  on  the 
surface.  The  roughness  of  the  a-Si:H  films  depends  most  critically  on  the 
sticking  coefficients  of  the  impinging  radicals,  and  secondarily  on  the 
deposition  rate.  We  find  that  SlH^  fractions  of  >50%  in  the  incident  flux  and 
deposition  rates  of  <3-4  A/s  produce  smooth  films  (roughness  <20%).  We  also 
find  a  correlation  between  subsurface  voids,  hydrogen  content,  and  high 
dihydride  densities  implying  that  hydrogen  is  trapped  in  voids  when  their 
passivated  sidewalls  are  covered  over.  Subsurface  evolution  of  the  hydrogen 
trapped  in  these  voids  must  likely  take  place  in  order  to  reduce  the  hydrogen 
fraction  and  void  density  of  slowly  growing  films.  In  camparing  our  results 
to  experiment,  we  find  that  passivation  of  dangling  bonds  on  the  surface  is 
likely  dominated  by  SiH^,  and  not  radical  species.  The  observation  the  high 
quality  films  are  obtained  in  moderate  pressure  pure  silane  plasmas  having  low 
dissociation  fraction  is  therefore  jaartly  explained  by  the  following  criteria. 
The  silcuie  pressure  must  be  sufficiently  high  that  insertion  and  abstraction 
reactions  of  radicals  with  the  feedstock  maintain  a  high  ratio  of  SlH_/SiH„  in 
the  radical  flux.  Low  dissociation  fractions  are  required  to  insure  that 
passivation  of  the  surface  by  SIH^  is  complete. 
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FIGURE  CAPTIONS 


FIG.  1.  Simulated  surface  of  an  a-Sl:H  film,  being  relatively  smooth, 
deposited  from  a  pure  silane  discharge  (250  mTorr,  power=100  mW/cm"®).  The 
film  (area  125  x  125  A^)  has  an  average  thickness  of  750  A  with  an  RMS  surface 
roughness  of  20  A. 

FIG.  2.  RMS  surface  roughness  as  a  function  of  composition  of  the  radical 
flux  (SlH^,  SIH2) .  PVD-llke  deposition  resulting  from  a  large  SIH2  fraction 
in  the  radical  flux  (  >0.5)  has  a  high  surface  roughness. 

FIG.  3.  RMS  surface  roughness  as  a  function  of  film  height  for  PVD-like 
conditions  (SiH2/SiH2  =  0.15/0.85  in  incident  radical  flux)  and  CVD-like 
conditions  (SiH2/SiH2  =  0. 1/0.9).  The  RMS  roughness  for  PVD-like  conditions 
Increases  with  increasing  film  thickness  due  to  shadowing  effects.  CVD-like 
conditions  quickly  reach  an  equilibrium  when  mobile  adsorbed  SiH^  radicals  on 
the  surface  fills  in  surface  roughness  at  approximately  its  rate  of 
generation. 

FIG.  4.  Film  oropertles  as  a  function  of  film  deposition  rate;  RMS  surface 
roughness  and  the  ratio  of  burled  hydride/di hydride  configurations.  The 
composition  of  the  incident  flux  (SiHg/SlHa  =  0.9/0. 1)  was  kept  constant  for 
all  cases  while  the  magnitude  of  the  flux  was  varied.  The  average  film 
thickness  for  all  cases  was  750  A. 
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FIG.  5.  The  atomic  hydrogen  fraction  and  density  defect  due  to  porosity  as  a 
function  of  growth  rate  for  the  deposition  conditions  described  in  Fig.  4.  The 
Increase  in  hydrogen  content  at  low  growth  rate  for  simulations  where  voids 
are  created  results  from  hydrogen  being  trapped  in  subsurface  voids.  The 
effect  may  be  magnified  by  the  model,  which  does  not  allow  hydrogen  evolution 
from  voids  after  burial. 

FIG.  6.  Hydrogen  content  and  hydrlde/dihydride  ratio  of  a-Si:H  as  a  function 
of  average  void  size. 

FIG.  7.  Fractional  distribution  of  voids  as  a  function  of  void  size.  The 

3 

average  void  volume  is  40.9  A  and  the  total  void  fraction  is  16.3%.  The 
sample  size  for  this  case  is  125  x  125  x  750  A  . 

FIG.  8.  Density  defect  as  a  function  of  film  thickness.  The  simulated 
deposition  conditions  are  the  same  as  for  Fig.  3.  The  solid  lines  denote  bulk 
porosity.  The  predicted  values  for  PVD-like  conditions  may  be  artificially 
low  due  to  the  assumption  in  the  model  that  cross  linking  can  take  place  with 
at  most  next  nearest  neighbors  (see  text).  When  the  effect  of  surface 
roughness  (R  .  dashed  lines)  is  Included,  PVD  films  have  a  larger  density 
defect  than  CVD-like  films. 

FIG.  9.  Atomic  hydrogen  fraction,  f^j,  as  a  function  of  film  height  for  the 
same  deposition  conditions  as  used  in  Fig.  3.  Thin  films  have  higher  since 
larger  fraction  of  their  atoms  are  contained  in  the  hydrogen  rich  near  surface 
layer . 
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FIG.  10.  a)  Computed  parameters  for  the  radical  flux  in  SiH^/Ar  plasmas  as 
obtained  from  the  model  described  in  Ref.  6.  b)  Deposition  rate,  hydrogen 
fraction,  as  buried  spin  density,  and  hydride/d ihydrlde  ratio  for  simulated 
a-Si:H  films  having  a  constant  magnitude  of  silane  radical  flux  for  the 
conditions  of  Fig.  10  a).  These  results  imply  a  degradation  of  film 
properties  when  diluting  SIH^  with  argon. 
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